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Mobility  of  heavy  metals  in  soils  is  of  environmental  significance  due  to  their 
toxicity  to  both  humans  and  animals.  In  general,  heavy  metal  mobility  is  low  because  of 
its  low  solubility.  However,  enhanced  heavy  metal  mobility  has  been  reported  under  both 
laboratory' and  field  conditions.  It  has  been  attributed  to  enhanced  solubility  of  heavy' 
metals  and  colloid-facilitated  metal  transport. 

In  this  study,  the  solubility  and  mobility  of  heavy  metals  were  examined  in  a Pb- 
spiked  sandy  soil  and  two  Pb-contaminnied  soils.  For  the  Pb-spikcd  soil,  water-flooded 
and  non-water-flooded  incubations  were  used  to  alter  soil  solution  chemistry  in  soil 
columns,  which  were  then  leached  with  de-ionized  water.  It  was  found  that  Pb 
concentrations  in  leachates  were  related  to  the  ratios  of  Fe  concentration  in  the  leachates 


to  Fe  concentrations  extracted  with  HC1.  Enhanced  Pb  mobility  occurred  only  when  the 
Fe  ratios  were  lower  than  a threshold  value  for  a given  soil.  The  two  Pb-contaminated 
soils  were  incubated  for  different  times  under  water  flooded  condition  to  alter  their  redox 


status.  Metal  solubility  was  examined  by  analyzing  the  pore  water  of  the  incubated  soil 
columns  whereas  metal  mobility  was  examined  by  leaching  the  columns  with  0.01  M 
CaCli.  The  data  showed  that  metal  solubility  in  pore  water  and  metal  mobility  in  CaClj 
solution  were  not  always  directly  related.  There  has  been  sufficient  evidence  that  a 
reduction  in  cation  exchange  capacity  (CEC)  occurs  with  Fe  reduction  dissolution. 
However,  the  consequent  redistribution  of  metals  between  solution  and  exchange  phases 
with  incubation  was  determined  not  only  by  the  magnitude  of  CEC,  but  also  the 
characteristics  of  metal  ions,  the  competing  co-ions  and  counterions  (ligands). 

Colloid  mobility  in  the  two  Pb-contaminatcd  soils  was  examined  by  switching  the 
influent  from  0.01  M CaCf  to  de-ionized  water.  In  addition,  a batch  dispensability  test 
was  conducted  by  dispersing  the  incubated  soils  in  0.01  M CaCli  and  0.06  M CaCl?.  and 
then  a relative  stability  ratio  (W)  of  an  incubated  soil  was  estimated  from  the  absorbency- 
time curves  of  a dispersion  in  the  solutions.  The  results  showed  that  colloid  mobility  is 
greatly  influenced  by  incubation  and  that  colloidal  metal  mobility  is  enhanced  when 
exchangeble  Ca  is  easily  replaced. 


CHAPTER  1 
INTRODUCTION 

Undemanding  heavy  metal  migration  in  order  to  accurately  predict  heavy  metal 

effective  soil  remediation  (McCarthy  and  Zachara,  1989).  Although  much  effort'has 
been  spent  on  developing  models  to  predict  heavy  metal  mobility  in  soils  (Ccdcrbcrg  et 
al..  1 985;  Sposito,  1 984).  these  models  often  underestimate  metal  migration  in  soils.  In 
most  of  the  models,  the  interaction  between  dissolved  metals  and  soil  matrix  is  mainly 
described  using  a distribution  coefficient  (lcD),  which  is  generally  considered  as  a constant 
and  usually  measured  in  the  laboratory  under  a particular  condition.  In  a dynamic  soil 
system,  however,  ko  is  not  a constant.  In  other  words,  solubility  of  metals  and  related 
mobility  vary  with  soil  conditions. 

Besides  dissolved  species,  mobile  colloids  con  also  act  as  vectors  to  transport 
sorbed  contaminants  such  as  heavy  metals  in  soils  (Seaman  ct  al.,  1 995).  Understanding 
mobility  of  potential  hcavy*mctal-bearing  colloids  is  critical  for  predicting  the  fates  of 

and  Elimelech.  1 996)  since  colloidal  particles  arc  ubiquitous  in  those  systems.  In  concept, 
the  deposition  and  release  of  colloids  in  porous  media  can  be  described  with  the  DLVO 
(Detjaguin-Landau-Verwey-Overbcek)  theory.  However,  it  is  difficult  to  apply  to  soil 

colloid/dissolved  metal  transport  model,  COMET  (Mills  et  al.,  1991),  has  been  developed 


lo  predict  colloidal  metal  mobility.  It  is  difficult  to  apply  to  soil  environments  because  of 
the  difficulties  in  determining  the  parameters  needed  in  the  model.  It  has  been 
recognized  that  colloid  mobility  in  soil  is  mainly  influenced  by  the  physicochemical 
characteristics  of  colloids  and  soils  (Seaman  et  al..  1995).  Unfortunately,  little  is  known 
quantitatively  about  the  relationship  between  colloidal  metal  mobility  and  the 
characteristics  of  soil  colloids  (Amrhein  et  al.,  1993;  Seaman  et  al..  1995;  Kaplan  et  al., 
1993). 

In  this  study.  Ute  solubility  and  mobility  of  heavy  metals,  colloid  mobility  and  its 
association  with  heavy  metals  are  studied.  The  main  objectives  of  this  research  are  as 
follows 

• To  characterize  the  conditions  at  which  enhanced  solubility  and  mobility  of 
heavy  metals  occur  with  incubation  in  soils, 

• To  characlerize  the  conditions  at  which  enhanced  colloid  mobility  occurs  with 
incubation  in  contaminated  soils, 

• To  characterize  the  conditions  at  which  enhanced  mobility  of  colloidal  heavy 
metals  occurs  with  incubation  in  contaminated  soils. 

Therefore,  there  are  two  major  subtopics  in  this  dissertation:  colloidal  metal 
mobility  and  dissolved  nrctal  mobility  (Figure  1-1 ).  The  understanding  of  colloid 
deposition,  release  and  association  with  heavy  metals,  is  first  reviewed  and  discussed  in 
the  next  chapter.  In  Chapter  5,  changes  in  colloid  and  colloid-facilitated  heavy  metal 
mobility  with  incubation  are  examined  in  two  contaminated  soils.  In  Chapter  6.  a 
fundamental  approach  is  taken  to  relate  colloid  mobility  to  colloid  stability  ratio  (W).  A 
relative  colloid  stability  ratio  is  defined  and  estimated  by  the  curves  of  absorbency-time. 


nap  of  this  dissertation 


On  the  other  hand,  the  relation  of  Pb  solubility  in  soil  and  Fc  partitioning  behavior 
after  incubation  is  examined  in  Chapter  3.  In  Chapter  4,  an  extension  of  the  study  from 
an  uncontaminated  soil  to  two  Pb-contaminated  soils  is  made  (Figure  1-1). 


CHAPTER  2 

COLLOID  DEPOSITION,  RELEASE  AND  ASSOCIATION  WITH  HEAVY  METALS 
IN  SOILS 


Inlioi  an 


mobilized  and  transported  a significant  distance  (Ryan  and  Elimelech,  1996;  Swanton, 
1995)  in  soil  environments.  As  a result , the  mobility  of  heavy  metals  associated  with  the 
colloids  may  be  enhanced  (Kaplan  ct  al.  1995;  Amrhein  el  al„  1993:  McCarthy  and 
Zachara.  1989;  Newman  ctal..  1993;  Mills  el  al„  I991;0uyangetal-  1996).  Generally, 
colloid  movement  is  faster  than  an  inert  tracer  due  to  the  effect  of  size  exclusion  (Nocito- 
Gobcl  et  al..  1996).  Therefore,  colloid-facilitated  metal  transport  has  been  considered  as 
one  of  the  important  mechanisms  causing  heavy  metals  to  move  much  faster  than 
expected  (Newman  et  al.,  1993),  which  has  attracted  much  attention  in  the  scientific 


Figure  2-1.  Schematic  of  potential  energy  profile  of  the  interaction  of  surfaces  with 
inclusion  of  van  der  Walls  attraction,  electrical  repulsion,  and  Bom  repulsion,  which 
shows  both  primary  and  secondary  minima  and  an  energy  barrier  as  well  as  the  zones  in 
which  release  and  deposition  take  place. 


their  large  surface  areas  and  high  transportability  in  porous  media.  Therefore,  this 
dissertation  will  mainly  focus  on  submicron  colloids  excluding  microbes. 

Since  there  is  a collection  of  literature  published  on  colloid  release,  transport  and 
deposition  in  recent  years  (Ryan  and  Elimelech,  1996;  Swanton,  1995;  Mills  ctal.,  1991; 
Ouyanget  al„  1996;  Adamczyk  cl  al.,  1983;  Kallay  etal.,  1987;  Elimelech  el  a),  1995; 
McCarthy  and  Zachara,  1989;  McDowell-Boyer  et  al.,  1986),  no  attempt  is  made  to 
present  a complete  critique  on  the  topics.  Instead,  focus  will  be  pul  on  several  parameters 
affecting  the  mobility  of  colloids  and  colloid-facilitated  heavy  metal  transport  in  soils. 

At  the  present  time,  various  processes  responsible  for  colloid  deposition  and 
release  in  porous  media  have  been  well  established  (Ruckenstein  el  at.,  1 976;  Adamczyk 
etal.,  1983;  Kallay  ctal.,  1987;  Elimelech  et  al„  1995;  Ryan  and  Elimelech,  1996). 

These  processes  can  be  theoretically  described  by  the  DLVO  (Dcrjaguin-Landau- 
Verwey-Overbeek;  Israelachvili,  1992)  theory,  which  states  that  the  net  surface 
interaction  energy  is  the  sum  of  the  interactions  of  electrical  double  layer  (EDL)  and  van 
der  Waals-London  forces  (WL),  which  vary  with  separation  distance  between  particles. 
As  particles  approach  each  other,  the  net  interaction  energy  experiences  a secondary 
minimum  first,  then  a primary  minimum  (6«w)  after  a maximum  energy  barrier  (6m„) 
(Figure  2-1 ).  After  that  point,  further  movement  towards  each  other  causes  drastic 

(Israelachvili,  1992).  In  principle,  particles  may  be  attached  or  deposited  when  the  net 
attractive  forces  are  close  to  either  the  primary  or  secondary  energy  minimum.  Both  the 
magnitude  of  the  energy  barrier  and  the  depth  of  the  primary  and  secondary  minimum  are 


affected  by  solution  chemistry.  In  order  for  a deposited  particle  to  be  released,  repulsive 
forces  have  to  be  generated  between  the  surfaces  of  particles  and  stationary  grains,  as  a 
result  of  changes  in  solution  chemistry  or  fluid  flow.  However,  the  understanding  of 
colloid  release  and  deposition  is  for  from  complete.  Colloid  deposition  under  favorable 
conditions,  in  the  presence  of  attractive  interactions,  can  be  predicted  reasonably  well, 
whereas  colloid  deposition  under  unfavorable  conditions,  in  the  presence  of  repulsive 
interactions,  can  not  be  predicted  by  the  DLVO  theory.  In  fact,  colloid  deposition  rates 
observed  experimentally  are  many  orders  of  magnitudes  greater  than  those  predicted  and 
are  independent  of  the  size  of  colloid  particles  (Elimelech  et  al..  1995).  Elimclcch  ct  al. 
(1995)  have  recently  presented  on  extensive  discussion  on  possible  explanations  for  these 
discrepancies.  However,  compared  to  colloid  deposition,  colloid  release  is  poorly 
understood,  which  will  be  emphasized  in  this  review. 

Transport  of  mobile  colloids  is  a critical  process  for  evaluating  colloid  mobility  in 
soils.  However,  there  arc  several  extensive  reviews  on  the  topic  (Elimelech  et  al.,  1995; 
Tien,  1989;  McDowcll-Boycr  el  al.,  1986;Ouyangetal„  1996),  so  this  topic  is  not 
emphasized  in  this  review.  Association  of  soil  colloids  with  heavy  metals  has  been  an 
important  issue  for  many  years  (McBride,  1 994;  Jones  ct  al„  1 975);  however,  until 
recently,  the  attention  has  not  been  paid  to  the  effect  of  such  an  association  on  colloid 
mobility  (Kretzschmar  et  al.,  1997a).  In  fact,  in  flotation  processes  widely  used  in  the 
mining  industry,  it  is  well  documented  that  aqueous  heavy  metals  can  be  potential 
determining  ions  or  specific  adsorbing  ions  (o  charged  surfaces  of  minerals  (Fuerstennau 
and  Palmer,  1976),  therefore,  such  an  association  with  colloid  surfaces  may  alter  the 


surface  charge  significantly.  The  potential  impact  of  the  association  of  heavy  metals  with 
colloids  on  colloid  mobility  will  thus  be  emphasized  in  this  review. 

More  recently.  Wan  and  Wilson  (1994a.b)  have  demonstrated  that  colloid 
particles  tend  to  sorb  onto  the  water-air  interface  in  soil  and  that  the  process  is  almost 
irreversible.  In  fact,  these  findings  have  sound  theoretical  and  experimental  basis  in  the 
flotation  discipline  (Williams  and  Berg,  1991).  They  concluded  that  colloid  deposition 
onto  the  interface  in  unsaluraied  porous  media  is  one  of  the  important  mechanisms 
responsible  for  the  retardation  of  colloid  transport,  except  under  extremely  high  flow  rate 
that  may  drag  tile  bubbles  with  the  flow  (Wan  and  Wilson,  1994a).  Incorporating  those 
results  into  a convective-diffusion  equation,  Coiapcioglu  and  Choi  (1996)  have 
developed  a model  to  predict  colloid  transport  in  unsaturated  porous  media.  These  results 
have  undoubtedly  advanced  colloid  transport  in  soils  significantly,  however,  they  will  not 
be  further  discussed  in  this  review. 

Most  of  the  current  reviews  focus  on  either  well  defined  porous  systems 
(Elimelech  et  al.,  1995;  Tien,  1989)  or  subsurface  aquifers  (Ryan  and  Elimelech,  1996; 
Swanton,  1 995;  McCarthy  and  Zachora,  1 989).  There  is  apparently  much  similarity 
between  soil  and  subsurface  or  well-defined  porous  systems.  Different  from  the  later 
two,  however,  soil  is  more  dynamic  in  its  solution  chemistry.  In  surface  soil,  pH,  ionic 
strength,  and  redox  potential  of  soil  solution  are  influenced  by  rainfall,  bioactivity,  and 
changes  in  land  use;  whereas  in  vndose  zone  they  are  mainly  impacted  by  the  fluctuation 
in  water  table  (Boul  el  al.,  1 997).  In  addition,  as  mentioned  above,  the  gas  phase  is  not 
negligible  and  the  interfaces  between  colloids  and  air  need  to  be  taken  into  account  along 
with  that  between  colloids  and  water.  Therefore,  colloid  transport  is  more  dynamic  and 


complicated  in  soils  compared  to  that  in  subsurface.  In  this  review,  we  will  emphasize 
how  soil  environments  afreet  colloid  mobility  and  what  are  the  possible  mechanisms. 

The  DL VO  Theory 

The  classic  DLVO  theoiy  has  been  extensively  used  to  describe  colloid 
interactions  and  stability  (Verwcy  and  Overbeck,  1948).  It  has  been  extended  to  describe 
colloid  deposition  and  release  in  porous  media  by  including  short-range  Bora  repulsion 
forces  (Ruckcnstian  et  a!.,  1976).  According  to  the  DLVO  theory,  the  total  interaction 
energy  between  colloids  and  stationary  solid  phases  (fowl)  is  the  sum  of  electrostatic 
repulsion  energy,  $fi,  arising  fiom  the  overlap  of  electrical  double  layers  and  attractive 
energy,  4>lw,  due  to  van  der  Walls-London  dispersion,  i.e.,  4>mt«r  ♦ri.+  ♦w..  which  is  a 
function  of  particle  separation  distance.  At  close  separation  distance,  a third  repulsive 
force,  the  Bom  repulsion,  is  also  important.  Thus,  the  generalized  curve  of  interaction 
energy  vs.  particle  separation  distance  shows  one  energy  maximum  (6™,)  & two 
minimums  (primary  minimum  (4tan)  and  secondary  minimum;  Figure  2-1J.  Aggregation 
or  coagulation  occurs  when  particles  collide  with  sufficient  kinetic  energy  to  overcome 
the  energy  barrier  (itwO  so  that  they  reach  a distance  within  the  4>m»-  On  the  other  hand, 
if  coagulated  particles  gain  enough  energy  under  perturbation  to  overcome  the  energy 
barrier  Ms™  - <>mi„)  with  increasing  separation  distance,  they  tend  to  be  released.  The 

interactions  between  particle  surfaces,  which  change  with  surface  properties  of  particles 


observation  and  the  prediction  based  on  the  DLVO  theory  (Elimelech  and  OMelia, 

1 990b).  Therefore,  in  recent  years,  various  updated  theories  of  colloid  stability  have  been 
proposed  (Shulepov  et  al.,  1995;  Adamczyk  et  al.,  1986).  The  accuracy  of  prediction  has 

described  by  the  DLVO  theory  in  some  systems  (Verwey  and  Overbeek,  1948 );  colloid 

plus  the  Bom  force  (Ruckcnstcin  and  Prieve,  1976).  Although  measured  collision 
efficiency  between  colloids  and  porous  media  surfaces  is  much  greater  than  the 
theoretical  prediction  (Elimelech  and  O'Melia.  1990b),  colloid  deposition  and  release  can 
be  qualitatively  described  by  the  theory  (Roy  and  Dzorabk,  1 996a).  Therefore,  the 
concept  of  the  DLVO  theory  is  correct  even  though  it  is  incomplete  (Swanton,  1 995). 

other  forces  exist  between  particles  (Israclachvili.  1992).  The  interactions  between  polar 
electron  acceptors  (Lewis  acids)  in  solution  and  polar  electron  donors  (Lewis  bases)  from 

addition  to  electrostatic  and  van  der  Waals-London  force  (Oss  et  al.,  1988).  It  can 
account  for  up  to  85%  of  interactions  between  soil  mineral  particles  in  aqueous  solution, 
as  such  the  DLVO  theory  combined  with  the  AB  force  adequately  describes  some 
anomalous  colloid  stability  (Oss  etal,  1990;  Wuetal.,  I994a,b).  One  advantage  of  this 
approach  is  the  AB  force  can  be  estimated  simply  by  the  contact  angles  between  isolated 
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particles  and  liquids  (Oss  et  al.,  1988).  Several  contact  angles  are  required  to  calculate 
the  AB  force  using  the  formula  of  Oss  et  al.  (1988).  The  measurement  needs  to  be  made 
in  at  least  thtee  different  liquids,  two  of  which  must  be  polar.  In  their  approach,  the  WL 
force  (VWl)  can  be  estimated  from  surface  tension  instead  of  using  assumptions  for 
natural  colloids  as  done  by  Ryan  and  Gsehwend  (1994a)  and  Roy  and  Dzombk  (1996a). 
Here  the  classic  DLVO  theory  combined  with  the  AB  force  is  termed  the  extended 
DLVO  theory.  Besides  the  AB  force,  other  non-DLVO  forces  between  particles,  such  as 
osmotic  and  stcric  interactions,  should  also  be  included.  The  extended  DLVO  theory 
should  also  correct  for  colloid  surface  roughness  and  heterogeneity  of  colloid 
compositions,  charges  and  sizes  (Elimelech  and  O'Melia,  I990a,b).  Unfortunately,  no 
models  yet  successfully  take  all  these  factors  into  consideration.  Nevertheless,  the 
DLVO  theory  provides  a solid  foundation  and  a useful  tool  to  describe  colloid 
interactions  (Swanton,  1995;  Ryan  and  Elimelech,  1996).  Therefore,  at  the  present  time, 
it  is  practical  to  describe  the  interactions  between  colloids  using  the  extended  DLVO 
theory  (6r^=  .+  $LW  + foa)-  The  extended  DLVO  theory  has  been  more  successful 
than  the  original  one  in  describing  colloid  behavior  (Oss  el  al.,  1 990;  Wu  et  al.,  1 994a.b), 
which  has  been  discussed  to  quantify  the  terms  in  the  equation  above  by  Oss  et  al.  (1 988) 
As  shown  in  Figure  2- 1 , in  the  presence  of  an  energy  barrier,  the  deposition  of  a 
particle  from  solution  onto  a surface  is  mainly  affected  by  the  magnitude  and  shape  of 
this  barrier  extending  from  its  maximum  into  the  solution  phase.  These  are  primarily  due 
to  electrostatic,  van  der  Waals-London  and  AB  force  contributions.  In  contrast,  particle 

surface,  which  is  greatly  influenced  by  additional  repulsion  at  short  separation  (Bom 


repulsion).  Over  this  region,  ihe  energy  profile  is  very  sensitive  to  the  specific 
characteristics  of  the  interacting  surfaces  and  intervening  liquid  layer.  However,  the 
irregular  dimensions  of  natural  soil  colloids  make  the  interaction  even  more  complicated, 
which  will  be  discussed  later. 


Charge  development 

Solid  particles  present  in  water  are  often  charged.  The  mechanisms  of  charge 
development  have  been  extensively  investigated  (Parks.  1965;  Sumner.  1992).  In 
general,  based  on  the  nature  of  solid  particles,  charge  development  can  be  divided  into 
three  categories:  1 ) lattice  substitution  is  most  common  in  soil  clay  particles  and  referred 
to  as  permanent  charge;  2)  specific  chemical  interactions  between  surfaces  and  solution, 
including  hydrolysis  of  surface  functional  group  (e.g.,  hydroxyl  and  carbonyl)  and 
chemical  adsorption;  and  3)  preferential  dissolution  resulting  from  preferential  hydration 
of  surface  atoms  In  addition,  based  on  the  contribution  from  electrolytes  in  solution  to 
the  charges  of  particles,  ions  can  be  classified  into  three  categories:  potential  determining 
ions,  specific  adsorbed  ions  and  indifferent  ions.  Potential  determining  ions  arc  the 
constituent  ions  of  solid  panicles  (e.g>,  H'  and  OH'  for  Fc  oxides,  Ca2*  and  CO?  for 
CaCOj)  and  their  concentrations  determine  primarily  the  surface  potential  of  particles 
(Parks,  1 967).  Specific  adsorbed  ions  can  change  the  magnitude  and  sign  of  the  surface 
charge  (e.g„  Ca1',  Pb2‘)  by  adsorbing  onto  the  surfaces  and  forming  a stem  plane  in 
electric  double  layer  (EDL)  theory  (Parks,  1975).  Indifferent  ions  adsorb  physically  to 
surfaces  and  change  the  magnitude  of  the  surface  charge  (e.g.  Na*.  Cl*),  and  form  a 
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interactions  with  the  surface  and  bulk  solution  (Parks,  1975). 

structure  of  colloidal  particles  make  their  charge  development  a complicated  matter.  This 
may  be  further  complicated  by  the  dynamic  nature  of  soil  solution  chemistry,  which  may 

numerous  studies  have  shown  that  the  charge  developing  process  can  still  be  described 
using  the  HDL  theory,  at  least  in  principle. 

Electrical  dQublc.lavcrXEDLI  in  metal  oxides-walcr  interface 

Electrical  double  layer  of  a dispersed  particle  consists  of  a charged  surface  and  a 
diffusive  layer  of  counter  ions  next  to  the  surface  based  on  the  double  layer  theory.  Since 
metal  oxides  may  play  an  important  role  in  colloid  mobility  and  have  been  well 
characterized  among  various  surfaces  of  soil  colloid  particles  (Schindler  and  Stumm, 
1987;  Dzombak  and  Morel,  1990),  wc  use  them  os  an  example.  When  specific  adsorbing 
ions  other  than  H"  and  OH  are  absent  in  a system,  the  charge  developing  process  on  the 
surfaces  can  be  described  as  follows  (Schindler  and  Stumm,  1987): 

SOH2'  = SOH+H*  K,i““  (2-1) 

SOH  =SO*  + H*  K^“  (2-2) 

Where  SOH  denotes  a surface  site  and  K.,“  and  Ka“  are  the  intrinsic  surface  acidity 

(2-3) 


K,,“  = [SOH2*)/[SOH](H*) 
Kal“  = [(SO'][H']/[SOH] 


(2-4) 


Schindler  and  Siumm  (1987)  summarized  the  inlrinsic  acidity  constants  (K.i”  and  ) 
of  various  metal  oxides  and  found  that  the  constants  are  generally  correlated  to  those  in 
solution.  However,  the  affinity  of  protons  to  the  surfaces  of  metal  oxides  is  much  more 
complicated  and  influenced  by  changes  in  surface  composition  and  structure  (Yoon  ct 
al.,1979;  Svcijensky,  1994;  Bleam.  1993).  In  addition,  the  values  of  these  constants 
depend  on  the  particular  electrostatic  model  adopted,  such  as  constant  capacitance, 
diffuse  double  layer,  and  triple  layer . However,  the  sum  of  Eq  (2-1)  and  (2-2)  can  be 
expressed  as  follows; 

SO’  + 2H*  = SOHj'  (2-5) 

Eq.  2-5  is  simply  related  to  the  point  of  zero  charge  (PZC  or  pHo).  which  is 
experimentally  measured  and  is  independent  of  electrostatic  models  for  the  solid-water 
interface.  When  the  concentration  of  positively-charged  surface  species  is  equal  to  that 
of  the  negatively-charged,  that  is,  a zero-charge  surface,  the  pHo  is  related  to  the  intrinsic 
equilibrium  constant 

pHo  =0.5  log  K.""  = 0.5  log  K.  (2-6) 

Eq  (2-6)  states  that  intrinsic  acidity  constant  (K.”‘)  is  equal  to  the  apparent  acidity 
constant  (Ka)  at  a surface  potential  of  zero.  It  not  only  indicates  that  pHo  is  a measure  of 
K“,  but  also  raises  the  question  that  the  bonding  of  protons  at  metal  oxide  surfaces  is 
more  analogous  to  the  bonding  of  similar  complex  in  an  aqueous  phase  (Schindler  and 
Siumm,  1 987)  or  in  bulk  crystal  structure  (Sverjensky,  1994).  The  correlation  between 
the  surface  acidity  constants  of  metal  oxides  and  those  in  solution  support  the  former 
point.  However,  the  scatter  on  such  correlation  is  substantial  for  some  mineral  particles 


(Schindler  and  Stumm,  1987).  More  recently,  Blesa  et  al.  (1990)  reviewed  the 


I (hydrous)oxides.  They  found 


ihennodynamic  data  for  ihe  adsorption  of  protons  on  metal 
that  the  enthalpy  values  associated  with  K,i'“  and  are  very  similar  and  are  much 
lower  than  those  of  the  hydrolysis  of  aqueous  metal  ions.  This  is  consistent  with 
Lyklemas  work  ( 1 987a),  in  which  he  describes  the  adsorption  of  potential-determining 
ions  as  a process  characterized  by  a single  enthalpy  value  instead  of  two.  These  indicated 
that  it  is  unrealistic  to  represent  acidic,  neutral,  and  basic  surface  groups  using  Eqs.  (2-1) 

(OHWOH2)*,  (mo  and  no  being  the  stoichiometric  index  at  pH„)  offers  a better 
description  of  the  involved  processes;  for  pH<pHo,  m<mo,  and  n>no;  for  pH>pHo,  m>mo 
and  n<no.  In  fact,  the  existence  of  adsorbed  water  molecules  and  surface  hydroxyl  groups 
has  been  observed  experimentally  (Tejedor-Tejedor  and  Anderson,  1986).  Based  on  this 

increase  in  the  number  of  coordination  (Blesa  et  al„  1990),  which  can  explain  why  the 
enthalpy  of  protonation  is  similar  on  cither  side  of  pHo.  Further,  Blesa  et  al.  (1990)  have 
pointed  out  that  there  is  a difference  in  Gibbs  Dee  energy  between  forming  aqueous  metal 
ions  and  surface  metal  complexes,  and  inferred  that  the  difference  may  have  resulted 
from  the  dehydration  of  the  adsorbing  ions.  They  argued  that  to  favor  adsorption. 


molecules  front  the  first  coordination  sphere  must  be  removed.  They  concluded  that  the 
easy  dehydration  of  anionic  species  should  account,  at  least  partially,  for  the  observed 
difference  in  the  hydrolysis  in  solution  and  the  surface  complex.  These  results  suggest 
that  the  solvation  of  metal  ions  must  be  included.  Following  this  line  and  combining 
crystal  chemistry  (Parks.  1 965.  Yoon  et  al„  1 979)  and  solvation  theory  (Sverjensky, 


1993).  Sverjensky  (1994)  has  shown  that  surface  protonalion  on  a wide  variety  of 
minerals  can  be  accurately  calculated  from  the  dielectric  constant  of  the  solid  and  the 
ratio  of  the  Pauling  bond  strength  to  the  cation-hydroxyl  bond  length  of  the  solid 
panicles.  This  suggests  that  much  more  emphasis  should  be  placed  on  the  analogy 
between  the  bonding  of  the  surface  protonated  species  and  the  bonding  in  the  underlying 
crystal  structure.  In  Sverjenskys  approach  ( 1 994),  the  properties  of  crystals  have  been 
numeralizcd  with  its  electric  constant  instead  of  the  qualitative  terms  such  as  kinks,  edge 
position,  or  adatom  positions.  The  details  on  modeling  or  calculating  pHo  can  be  found 
in  the  literature  (Porks,  1965;  Yoon  ct  al„  1979;  Sverjensky,  1994;  Bicam,  1993). 

Fokkink  et  al.  (1987,  1989)  have  demonstrated  that  pH-  and  temperature- 
congruencies  exist  in  the  surface  charge  developing  processes  of  various  metal  oxides.  In 
a systematic  series  of  experiments  in  which  the  surface  charge  ao  of  a number  of  oxides 
[TiOj  (rutile),  RuO;  and  o-FejOj  (haematite))  in  aqueous  solution  of  KNO3  was 
measured  as  a (unction  of  pH,  ionic  strength  and  some  other  variables.  When  00  is 
plotted  as  a function  of  pH-pHo,  where  pHo  is  the  PZC,  the  curves  coincide  for  the  three 
oxides  at  three  ionic  strengths  even  though  they  are  different  in  pHo-  In  other  words,  the 
individual  identity  of  metal  oxides  makes  no  difference  to  the  S-shapcd  plot  of  00  vs.  pH- 
pHo.  Similarly,  in  the  measurements  of  the  temperature  dependence  of  surface  charge 
development,  changes  in  temperature  only  affect  the  positions  of  PZC,  not  the  trend  of 
the  surface  charge  developing  process.  They  conclude  that  the  electrical  double  layers  in 
metal  oxides  can  be  functionally  divided  into  a specific  and  a generic  part,  with  respect  to 
the  natures  of  oxide  and  electrolyte.  The  specific  part,  which  is  determined  by  the 
specific  interactions  of  a surface  with  protons  in  solution,  determines  the  PZC.  This  is 


Df  Svcijcnsky  <1994)  in 


controlled  by  the  solution  side  of  the  double  layer,  determines  the  surface  charge 
development  once  the  surface  charge  deviates  from  pHo.  It  can  be  described  well  by  the 
Ncmst  law  (ip0  =2.303  RT/F  (pH-pHo).  The  significance  of  this  result  is,  first,  that  it  has 
provided  a thermodynamic  basis  for  double  layer  model  that  divides  the  charged  interface 
into  a charged  surface  and  a diffuse  layer  near  the  surface;  second,  that  the  Nemst  law 
can  be  used  to  describe  the  charge  developing  process  in  metal  oxides.  It  has  been 

interface  is  a constant  (Blesa,  1988).  In  a practical  sense,  the  potential  can  be  treated  as  a 
constant  if  the  amount  of  If  in  the  interfacial  layer  is  much  higher  than  a change  caused 
by  the  charge  developing  process  (Blesa,  1988).  However,  for  metal  oxides  only  3-10 
sites  per  100  A (1  nm!)  are  usually  found  (Sigg  et  al.,  1981)  and  in  the  experimental  pH 
range  up  to  60%  of  the  sites  may  be  charged  (Blesa.  1988).  However,  if  the  notation  of 
the  binding  site  S(OH  WOH2)*>  we  discussed  above  is  more  representative  of  the  real 

neutral  species  (HjO)  on  surfaces,  and  thus  the  adsorbed  IT  would  not  change  much  with 
surface  charge.  This  may  result  in  the  Nemstian  behavior. 

In  the  presence  of  specific  adsorbing  ions  (SAls)  such  as  Ca1'  and  Mg2'  in 
solution,  the  charge  developing  process  becomes  more  complicated  (Lyklema,  1984; 
Ardixzonc,  1982).  Generally,  the  PZC  will  be  shilled  compared  with  the  one  in  the 
absence  of  those  SAls.  Lyklema  (1984)  defined  the  PZC.  even  in  the  presence  of  SAls, 

surface  charge  o0  = F(TH+-r0H-).  The  major  basis  for  the 


; definition  is  the 


sensitivity  of  oxides  to  H'  and  OH\  Consequently,  SAIs  belong  to  the  Stem  charge. 
Therefore,  the  PZC  in  the  presence  of  SAIs  remains  fully  defined.  In  their  definition 
PZC  is  sdil  the  pH  where  o0  =0.  If  the  cation  adsorbs  on  to  the  Stem  layer  PZC  is  lower 
than  the  PZC  in  the  absence  of  SAIs  because  cations  in  the  stem  layer  favor  the 
adsorption  of  OH"  over  H\  so  that  a lower  pH  is  required  to  restore  the  H*-OH' 
adsorption  balance.  Similarly,  specific  adsorption  of  anions  leads  to  PZCs  above  the 
pristine  value. 

In  fact,  in  the  presence  of  SAIs,  the  common  intersection  point  (CIP)  of  titration 
curves  with  various  ionic  strengths  has  also  been  observed.  Similar  to  the  situation  in  the 
absence  of  SAIs,  the  CIP  is  the  equal  compensation  point  (Lyklema,  1 984),  at  which 
charge-compensating  cations  and  anions  have  an  equal  affinity  to  the  surfaces.  From  this 
point  of  view,  they  explained  CIP  *PZC  in  the  case  of  specific  adsorption  as  follows:  at 
the  PZC  adsorption  of  metal  cations  are  favored  over  that  of  anions  if  cations  have  a 
higher  affinity  to  tile  surface.  In  order  to  reach  a situation  where  the  intrinsic 
adsorbability  of  cations  and  anions  is  identical  (equal  compensation  point),  a positive 
charge  on  the  surface  has  to  be  developed.  This  charge  should  be  more  positive  with 
higher  specific  affinity  of  cation  to  the  surface.  Once  the  equal  compensation  point  has 
been  reached,  further  increase  in  concentrations  of  cations  and  anions  do  not  lead  to 
further  shift,  i.c.,  all  successive  curves  at  increasing  concentration  must  pass  through  the 
point.  Similar  reasoning  applies  to  a situation  at  which  anions  have  a higher  affinity  over 
cations.  This  principle  is  important  to  understand  the  impact  of  soil  solution  chemistry 
on  soil  colloid  charge  development.  Following  this  point,  it  is  easy  to  understand  the 
effects  of  changes  in  species  of  cation  and  anion  in  solution  on  PZC  and  CIP. 


In  the  discussion  above,  we  have  been  dealing  wilh  well-defined  crystal  surfaces. 


in  which  the  dissolution  can  be  neglected  using  surface  complcxalion  approaches. 
However,  solubility  of  metal  oxides  has  to  be  taken  into  account,  especially  in  the  case  of 
amorphous  (hydroxides.  Under  such  conditions,  the  adsorption  and  desorption  of  the 
hydrolyzed  metal  complex  ions  may  become  an  alternative  to  determine  surface  charge 
other  than  the  surface  complexation  of  H*  and  OH”.  Based  on  the  minimum  solubility 
theory  (Parks  ct  al„  1962),  PZC  should  be  identical  to  the  isoelectric  point  (IEP)  of  an 
aqueous  solution  suspending  the  particles.  The  IEP  is  defined  as  the  pH  resulting  from  an 


1965),  which  is  often  found  at  the  pH  of  minimum  solubility  of  a solid.  Blesa  et  al. 
(1997)  have  found  that  the  PZC  and  IEP  for  metal  oxides  do  not  always  coincide; 
however,  the  deviation  from  IEP  (PZC  - IEP)  can  be  explained  well  if  considering  the 
fact  that  the  complex  cations  and  anions  require  different  dehydration  energy  when  they 
are  transferred  from  solution  to  surface.  They  explained  that  monomeric  cations  hydrate 
more  strongly  than  anions  in  aqueous  solution,  so  the  shift  should  be  small  and  negative 
when  the  charge  is  determined  by  monomeric  surface  complexes  and  the  minimum 
solubility  theory  is  generally  valid.  On  the  other  hand,  for  polymeric  species,  more 
specific  behaviors  are  expected.  Relative  large  positive  or  negative  shifts  may  result  if 
the  charge  surface  complexes  are  polymeric.  Obviously,  the  underlying  assumption  of 
these  approaches  is  that  the  surface  active  site  is  more  analogous  to  their  solution  species. 
'This  may  be  more  suitable  to  soil  environment  where  various  amotphous  minerals  are 
abundant  and  aqueous  species  are  more  understood  than  surface  species. 


In  general,  surface  charge  of  soil  panicles  can  be  classified  into  two  types: 
permanent  and  pH-depcndenl  charges.  Detailed  characterization  of  soil  panicle  charges 
has  been  developed  by  Charlct  and  Sposilo  (1987).  Anderson  et  al.  (1991)  and  Chorover 
and  Sposito  ( 1 995).  Charge  developing  processes  on  soil  panicles  are  more  complicated 
than  those  on  metal  oxides  because  the  surface  composition  and  structure  of  soil  particles 
varies  greatly  and  PDIsare  not  limited  to  just  H*  and  OH‘.  However,  similar  phenomena 
such  as  CIP  have  been  observed,  which  is  referred  to  as  the  point  of  zero  salt  effect 
(PZSE).  It  is  often  found  that  PZSE  deviates  from  CIP  (Chorover  and  Sposilo,  1 995). 

Soil  is  a mixture  of  various  minerals  and  organic  substances.  Funher,  there  con 
be  significant  amounts  of  amorphous  materials,  such  as  amorphous  Fe  and  AI.  It  is  well 
documented  that  the  existence  of  hydrolyzed  species  of  Al  and  Fe  have  significant  impact 
on  the  pH-dcpendcnt  charge  of  soil  particles  (Chorover  and  Sposilo.  1995,  Parker  et  al.. 

1 979),  and  their  complexation  with  organic  matter  further  complicates  the  charge 
developing  process  (Chorover  and  Sposito,  1 995).  In  addition,  mineral  dissolution  may 
be  significant  for  some  soils  and  may  affect  the  charge  developing  process  greatly  (Blesa 
el  al„  1997).  Combining  all  these  factors  with  the  dynamic  nature  of  solution  chemistry, 
one  may  conclude  at  the  present  time  it  is  impossible  to  accurately  predict  proton  surface 
charge  developing  process. 

In  principle,  however.  Blesa  et  al'sapproach  (1990;  1997)  is  more  promising. 
They  emphasize  the  role  of  partitioning  of  charged  species  between  solid  surface  and 
solution  during  surface  charge  developing  processes,  which  is  driven  by  the  difference  in 
Gibbs  free  energy  between  a charged  species  in  solution  and  on  surfaces.  The  underlying 
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assumption  of  this  approach  is  that  the  charged  species  in  solution  and  on  surface  are 
similar.  This  approach  has  relative  sound  theoretical  basis  and  is  consistent  with 
experimental  evidence. 

Charge  development  in  mobile  colloids 

Considerable  amount  of  research  has  been  conducted  to  understand  charge 
development  of  water  dispersive  and  mobile  colloids  (Heil  et  ah.  1993a,b;  Kretzschmar, 
1993;  1996).  It  is  commonly  observed  that  the  stability  of  colloidal  clays  in  soils  is  many 
times  greater  than  that  of  comparable  reference  clays.  Also,  clays  isolated  from  surface 
soils  were  found  to  be  more  dispersive  than  clays  from  the  subsurface  horizons  of 
kaolinitic  soils  in  the  Southeastern  USA.  Factors  contributing  to  the  high  dispersibility  of 
soil  clays  include  the  presence  of  adsorbed  organic  matter,  minor  quantity  of  smectite  in 
kaolinitic  clay,  and  larger  surface  roughness  of  soil  clays  compared  with  unweathered 
reference  clays. 

Organic  matter  is  an  important  source  of  negative  charges  in  soils.  Well 
decomposed  humus  may  have  a CEC  > 300  emol  kg'1  humus,  which  is  considerably 
greater  than  that  of  clays  such  as  kaolinitc  (3-1 5).  illite  (30-40)  and  montmorillite  (80- 
1 50).  It  has  been  estimated  that  20-70  % of  the  CEC  of  many  soils  attributes  directly  to 
the  soil  organic  matter  alone  (Vaughan  et  al„  1984;  Stevenson.  1982).  In  soils,  dissolved 
organic  carbon  (DOC)  tends  to  adsorb  onto  solid  particles  such  as  clays  and  metal  oxides, 
which  is  driven  by  ligand  exchange,  multivalent  ion  bridging.  Van  der  Waals  force  and 
hydrophobic  interactions  between  DOC  and  the  solid  minerals  (Murphy  et  al„  1995).  In 
particular,  for  Fe  and  Al  oxides,  when  the  pH  of  a system  is  lower  than  their  PZC, 
adsorption  of  DOC  onto  the  minerals  increases  significantly  because  of  the  electrostatic 
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attraction  between  the  two.  resulting  in  a significant  impact  on  their  surface  charge. 

There  are  numerous  studies  of  the  effects  of  organic  coating  onto  metal  oxides  (Tipping 
et  al..  1 982)  and  clay  minerals  (Kretzschmar.  1 996)  on  surface  charge.  Heil  and  Sposito 
(1993a,  b)  found  flocculation  of  illitic  soil  colloids  with  organic  coatings  increased  with 
pH  in  Ca  solution,  which  is  contrary  to  their  flocculation  behavior  when  organic  matter  is 
removed  by  HiO;.  They  concluded  that  the  competition  between  proton  and  Ca  onto  the 
acidic  functional  groups  of  organic  molecules  is  essential  to  the  charge  development  on 
coated  colloids.  A similar  result  was  also  found  by  Kretzschmar  (1 993; 1 997)  in  humic- 
coated  kaolinic  soil  clay.  On  the  other  hand,  it  has  been  reported  that  the  electrical 
properties  of  the  quartz  surface  dispersed  in  river  water  (Whitray  Bech,  UK)  is 
determined  essentially  by  their  interactions  with  inorganic  cations  such  as  Al  and  Ca 
instead  of  organic  matter  (Findlay  et  al.,  1 996).  In  concept,  however,  these  results  are 
consistent  with  others,  i.e.,  surface  electric  properties  of  minerals  are  a result  of 
interactions  both  in  the  interface  of  particle-solution  and  in  solution,  including  particle- 
organic  matter,  organic-inorganic  ions,  particle-inorganic  ions,  etc.  Particularly,  types  of 
cations,  such  as  proton  and  multivalent  cations,  and  their  concentrations  are  important  for 
the  charge  development  on  particles  in  addition  to  organic  matter. 

It  is  unrealistic  to  model  surface  charge  development  on  soil  particles  with 
reasonable  certainty  at  the  present  time.  However,  it  is  evident  soil  particles  can  be 
treated  as  assemblages  of  crystalline  and  amorphous  minerals  and  organic  residuals; 
Sophisticated  models  such  as  the  triple  layer  model  do  not  make  much  sense  in 
describing  the  electrified  interfaces  of  such  particles  considering  that  the  particle  surfaces 
are  irregular  and  the  Stem  layer  often  moves  inward  inside  of  the  physical  boundary  of 


(he  particles.  As  a first  approximation,  on  the  other  hand,  the  interface  of  soil  particles  is 
more  reasonable  described  by  diffuse  double  layer  with  a specific  part,  i.e.,  a charged 
surface,  and  a generic  part,  i.e.,  a diffuse  layer.  The  surface  charge  developing  process 
may  be  described  by  the  partitioning  of  charged  species  between  solution  and  surfaces 
assuming  the  charged  species  are  similar  between  the  two  phases. 

Hydration 

Water  molecules  may  be  orientated  surrounding  the  immersed  particles,  forming 
"Hydration  shells”  (Clifford.  1 975).  Once  the  particles  with  such  shells  approach  each 
other,  additional  forces  between  the  particles  arise.  The  origin  of  these  forces  is  the 
interactions  between  polar  electron  acceptors  (Lewis  acids)  in  water  and  polar  electron 
donors  (Lewis  bases)  from  colloidal  particles,  and  termed  AB  forces  as  discussed 
previously  (Oss  et  al,  1987;  1988).  Surface  electron  donicity  of  colloid  particles  is 
crucial  for  water  molecules  to  be  orientated  along  the  surface.  Strong  surface  electron 
donicity  causes  more  ordering  of  water  molecules  orientating  along  the  surfaces,  resulting 
in  a repulsive  force  between  the  particles.  Weak  surface  electron  donicity  causes  less 
ordering  of  water  orientating,  resulting  in  an  attractive  force.  Ossctal.  (1990)  have 
demonstrated  between  particles  dispersed  in  water,  the  interaction  free  energy  from  AB, 
whether  repulsive  or  attractive,  is  commonly  as  much  as  1 00  times  greater  than  LW 
energy,  and  1 0 or  more  times  greater  than  EL  energy  at  close  range  ( I -5  nm).  It  can 
account  for  up  to  85%  of  the  interactions  between  mineral  particles  in  aqueous  solution, 
such  that  the  DLVO  theory  combined  with  the  AB  force  can  adequately  describe  some 
anomalous  colloid  stability  (Oss  et  al,  1990;  Wu  et  al,  I994a,b).  Wu  et  al.  (1994a,b)  have 
experimentally  evaluated  the  interactions  between  the  particles  of  montmonttonite  and 
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calcilc  by  examining  flocculation  of  ihcir  suspensions  while  increasing  Ca2' 
concenlralion.  They  concluded  that  Ca!*  not  only  causes  a decrease  in  panicle  £ - 
potential,  but  it  also  drastically  lowers  the  electron  donicity  of  the  polar  surface  of  the 
panicles,  resulting  in  an  AB  attractive  force,  which  is  more  responsible  for  the 
flocculation  of  the  suspensions.  Obviously.  SAIs  such  as  Co2*  alter  not  only  £ - potential 
when  it  adsorbs  onto  a panicle,  but  also  decreases  the  surface  electron  donicity.  This 
result  means  that  surface  electron  donicity  and  £ - potential  of  panicles  are  related. 

Size  development  of  mobile  colloids 

The  behaviors  of  colloids  are  size-dependent.  Individual  colloidal  panicles  may 
aggregate  depending  on  the  magnitude  of  the  energy'  barrier  of  their  interaction  (Figure  2- 
1 ),  which  is  similar  to  colloid  deposition  that  will  be  discussed  later.  Briefly,  if  the 
energy  barrier  of  interactions  is  low,  fast  coagulation  occurs,  in  which  case  the  diffusion 
rate  of  colloids  is  the  limiting  factor,  the  aggregation  rale  is  otherwise  controlled  by  the 
interaction  of  colloidal  particles.  These  two  mechanisms  result  in  lower  and  higher 
fractal  dimensions  of  aggregates,  respectively  (Riscovic  et  al„  1 996).  However,  direct 
measurement  and  theoretical  simulation  of  colloid  size  distribution  in  natural  system  have 
shown  that  the  concentration  of  colloids  < 0. 1 pm  should  be  negligible  because  of  their 
almost  instant  coagulation  (Filella  et  al,  1993;  Buflle  et  al.,  1995).  Kaplan  et  al.  (1997) 
examined  the  possibility  of  aggregation  of  mobile  colloids  from  a reconstructed  soil 
profile.  They  analyzed  the  size  distribution  of  panicles  in  the  suspensions  with  and 
without  sonificalion,  which  is  supposed  to  break  down  aggregation  in  the  suspension. 
They  found  that  the  suspension  without  sonificalion  exhibits  almost  the  same  bimodal 
particle-size  distribution  as  the  one  with  sonificalion.  This  holds  for  soil  suspensions 
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collected  from  two  different  soils:  a loamy  sand  and  a sandy  soil  in  Southeastern  USA. 
Their  results  suggest  that  colloid  aggregation  in  soil  pore  space  is  very  limited  especially 
when  soil  colloid  particles  are  highly  charged  (Kaplan  et  al.,  1995;  1997).  Interestingly, 
similar  bimodal  size  distribution  of  mobile  colloids  has  been  reported  by  Ronen  et  al. 
(1992).  who  characterized  the  suspended  particles  collected  Dorn  groundwater  in  a 
coastal  plain  phreatic  aquifer  of  Israel.  Kaplan  cl  al.  ( 1 996)  explained  that  smaller 
particles  are  less  likely  to  be  strained  during  transport  in  porous  media  based  on  the 
Duration  theory  (Yao  et  al.,  1971).  This  suggests  that,  in  the  bimodal  distribution, 
colloids  of  smaller  size  (<  0.4  pm)  may  have  resulted  from  long-distant  translocation  to 
sampling  point,  while  those  of  larger  size  (<1  pm)  may  come  from  local  dispersion  at  the 
sampling  point. 


The  Characteristics  of  Porous  Media 

Mobile  colloids  interact  with  the  surfaces  of  porous  media,  which  is  more 
important  than  those  between  colloids  in  limited  pore  space.  Therefore,  pore  structures 
and  their  physical-chemical  surface  characteristics  have  tremendous  influence  on  colloid 
mobility.  Extensive  studies  have  been  done  in  the  systems  consisting  of  clean  surfaces  of 
porous  media  (Adamczyk  et  al.,  1983;  Kallay  et  al..  1987.  Elimelech  et  al,  1995). 
However,  large  discrepancies  in  colloid  deposition  rate  under  unfavorable  deposition 
conditions  between  theoretical  predictions  and  experimental  observations  have  been 
found.  This  has  been  attributed  to  the  hydrodynamic  effect  and  surface  heterogeneity  of 
the  porous  medio,  such  as  surface  roughness  and  local  charge  heterogeneity.  This  has 
been  considered  to  be  the  most  promising  approach  to  understanding  the  interactions 
between  colloids  and  porous  media  in  clean  systems  (Song  et  al.,  1 994).  There  are  some 


exhaustive  reviews  on  this  topic  (Swanton,  1 995;  Ryan  and  Elimelech,  1 996;  Elimclech, 
ctal.,  1995). 

In  a soil  system,  the  porous  media  is  an  assemblage  of  a wide  range  of  particles, 
which  vary  greatly  in  mineralogy,  surface  composition  and  dimensions.  Different  from 
pure  systems,  where  the  surface  properties  of  porous  media  may  be  far  different  tfom 
those  of  colloids,  the  surfaces  of  soil  porous  media  may  at  least  in  part  consist  of 
potentially  mobile  colloids.  In  sandy  soil,  a small  amount  of  clay-size  particles  tends  to 
pack  closely  along  large  sand  grains  driven  by  energy  minimum.  This  arrangement 
among  a grain  and  colloidal  panicles  has  also  been  confirmed  by  fractal  analysis  (Bartoli 

be  treated  as  those  between  colloidal  panicles,  as  a first  approximation,  which  has  been  a 
long-established  concept.  Swanton  (1995)  has  suggested  that  in  an  undisturbed  soil 
profile,  the  active  site  of  the  porous  media  may  be  completely  occupied  by  mobile 
colloids.  Ryan  and  Gschwend  (1994a)  have  successfully  related  the  observed  colloid 
release  rales  from  an  Fe  oxide-coated  sand  columns  to  EXP  (4mu-tmin)fkBT.  Their 
estimation  of  the  energy  barrier  is  based  on  the  interactions  between  colloids 

rather  than  those  between  colloids  and  porous  medio  surface.  The  underlying  assumption 

Deposition  of  Colloids 

Theoretical  background  in  well-defined  porous  media 

Colloid  mobility  through  porous  media  is  mainly  determined  by  the  net  rate  of 
colloid  deposition  and  release.  Generally  speaking,  deposition  of  colloids  onto  porous 
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media  can  be  viewed  as  a two-slep  process:  Ihc  transport  of  colloids  from  bulk  solution  to 
the  proximity  of  the  surfaces  and  then  their  attachment  to  the  surfaces,  which  depends 
upon  the  nature  of  particle-surface  interactions.  If  the  energy  barrier  Ibnis.)  is  < 0,  the 
Qux  of  colloid  transport  is  equal  to  that  of  deposition.  In  oilier  words,  the  collection 
efficiency  of  the  collector  (porous  media)  is  1 00%.  The  energy  barrier  <)„„  must 
otherwise  be  overcome  in  order  for  a colloidal  particle  to  attach  to  the  surfaces  (Figure  2- 
1).  Therefore,  deposition  rate  constant  (k*p)  is  related  to  the  energy  barrier  with  k*p 
« EXP  (-jw/koT),  where  ka  is  the  Boltzmann  constant  and  T is  absolute  temperature 
(Ruckenstein  et  al„  1 976).  The  major  mechanisms  of  colloid  transport  to  a collector  ate 
inertial  impaction,  interception,  sedimentation,  electrostatic  forces.  Brownian  diffusion, 
and  straining.  Collection  efficiency  contributed  from  each  of  them  has  been  desoribed  by 
Tien  (1989),  and  the  details  will  not  be  repeated  here.  In  porous  media,  however,  many 
factors  mentioned  above  may  be  operative  simultaneously.  The  overall  collection 
efficiency  can  be  theoretically  calculated  based  on  the  filtration  theory  (Yao  et  al..  1971). 
However,  the  great  discrepancy  between  prediction  and  experimental  observation  under 
unfavorable  deposition  conditions  has  been  observed.  This  has  been  attributed  to  the 
distribution  of  surface  and  physical  properties,  surface  charge  heterogeneity  of  solids, 
surface  roughness,  interfacial  electrodynamics  and  colloid  deposition  in  secondary 
minima  (Elimelech  et  al.,  1995).  In  order  to  take  this  discrepancy  into  account,  collision 
efficiency  (a=q/qo)  is  used,  where  q and  arc  observed  and  calculated  collection 
efficiencies,  respectively.  As  particles  deposit  on  a collector,  affer  an  initial  stage,  the 
deposition  rate  will  change  depending  on  the  nature  of  particle-particle  interaction.  If  the 
net  interaction  is  repulsive,  die  collector  surfaces  become  progressively  occluded  as 
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panicles  accumulate  and  colloid  deposition  tale  declines  accordingly.  This  surface 
exclusion  phenomenon  is  termed  blocking.  There  are  some  excellent  recent  reviews  on 
this  topic  (Ryan  and  Elimclcoh.  1996;  Elimelech  et  al.,  1995;  Swanton.  1995). 

Colloid  deposition  in  soil 

Transient  phenomena  in  colloid  deposition.  There  are  no  clean  or  pure  surfaces  in 
soils.  Therefore,  the  results  obtained  in  well-defined  porous  media  cannot  be  directly 
extrapolated  to  describe  colloid  deposition  in  soils.  In  fact,  in  any  given  time  in  most 
soils,  the  most  favorable  deposition  sites  are  likely  to  readily  be  occupied.  Alternatively, 
if  there  is  influx  of  mobile  colloids  as  a result  of  some  chemical  or  mechanical 
disturbance,  there  may  be  limited  deposition  during  the  early  stages  of  exposure.  This 
occurs  because  favorable  deposition  sites  arc  already  filled,  which  is  different  from  clean 
porous  media.  Over  a prolonged  period  of  exposure  for  a colloid  flux,  on  the  other  hand, 
soil  porous  media  may  develop  surfaces  that  are  more  similar  to  those  of  mobile  colloids 
(Swanton,  1995),  which  has  been  discussed  before.  In  soil  porous  media,  therefore, 
colloid  deposition  may  be  a transient  phenomenon. 

Dispcrsable  fraction  of  soil  clay  resulted  from  vigorously  shaking  in  water  can  be 
a relative  measure  of  the  amount  of  potentially  mobile  colloids.  Miler  ct  al.  (1986) 
observed  the  amount  remaining  in  suspension  after  36  h of  shaking  was  highly  correlated 
to  surface  soil  loss  in  the  Southeaster  US  under  high  intensity  rainfall.  Kaplan  et  al. 
(1997)  have  found  mobile  colloids  from  two  soils  are  similar  in  mineralogy  to  the  water 
dispersible  clays  and  that  they  arc  many  orders  of  magnitude  lower  than  the  amount  of 
water  dispersible  clay  in  the  soils.  In  other  words,  there  is  significant  amount  of 
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potentially  mobile  colloids  loosely  attached  to  the  soil  porous  media,  however,  during 
this  period  they  function  as  a part  of  stationary  porous  media.  They  may  be  detached 
from  the  porous  media  (mobilized)  once  they  are  exposed  to  a chemical  or  hydrological 
perturbation.  I have  referred  to  this  stage  a transient  phenomenon  in  colloid  deposition. 

In  fact,  this  is  not  a new  concept  at  all;  however,  it  has  to  be  emphasized  to  understand 
the  distinctively  different  features  of  colloid  mobility  in  soil  porous  media.  Take  surface 
soil  for  an  example.  A perturbation  such  as  rainfall  occurs  on  the  soil  surface  and  moves 
downward.  Such  a perturbation  may  result  in  a release  of  deposited  colloids,  and  then  the 
mobilized  colloids  move  down  over  the  surfaces  of  porous  media  consisting  of  the 
deposited  colloids.  Therefore,  at  the  front  of  the  downward-moving  perturbation,  there 
ate  two  opposite  processes  operating  simultaneously:  deposition  and  release  of  soil 
colloids.  The  former  is  mainly  influenced  by  the  blocking  effect,  and  the  latter  will  be 
discussed  in  the  following  section. 

Size  Straining  in  deposition.  If  the  effective  size  of  a colloid  particle  is  larger  than 
the  smallest  pore  through  which  fluid  is  flowing,  the  particle  is  retained  by  porous  media, 
which  is  referred  to  as  size  straining.  Straining  occurs  in  granular  bed  filtration  if  the 
ratio  of  the  suspended  particle  diameter  to  the  grain  diameter  is  greater  than  about  0.2 
(Hetzig  el  al.,  1970).  Straining  is  one  of  the  most  important  mechanisms  of  colloid 
deposition  in  soils  (Seta  et  al.,  1 997;  Jacobsen  et  al..  1 997;  Kaplan  et  al.,  1 997).  Seta  et 
al,  (1997)  have  examined  the  transportability  of  water  dispersible  clay  through  intact  soil 

concentrations  of  colloids  in  leachates  along  with  pH,  total  exchangeable  bases,  cation 
e capacity,  organic  carbon  content,  etc.  They  concluded  the  relative  low  colloid 
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transportability  is  attributed  to  larger  colloid  size,  i.c.  size  straining.  The  straining  may 
manifest  itself  in  another  way.  Jacobsen  et  al.  (1997)  examined  the  transportability  of 
illite  though  soil  columns  containing  mactoporcs,  and  found  illite  concentrations  moving 
through  soil  columns  are  significantly  related  to  the  sizes  of  active  macropores  in  soil 
columns.  However,  no  significant  difference  in  mobility  through  the  columns  between 
the  illite  and  that  coated  with  humic  acid  was  found.  This  suggests  size  straining  is  the 
predominant  process  in  determining  the  deposition  over  the  zeta  potential  and  steric 
effect,  resulting  from  the  coatings  of  humic  acid.  In  their  study,  this  is  further  confirmed 
by  the  fact  that  colloid  concentrations  in  leachates  increase  with  the  intensity  of 
infiltration,  as  they  argued,  a high  intensity  of  infiltration  depresses  the  effect  of  straining. 
However,  this  is  contrary  to  the  results  of  Kretzschmar  ct  al.  (1995).  Kretzschmar  ct  al. 
(1995)  found  colloids  coated  with  humic  acid  result  in  a much  lower  collection  efficiency 
in  columns  packed  with  "clean"  saprolite  particles.  In  fact,  it  is  this  contradiction  that 
demonstrates  the  difference  of  colloid  deposition  between  soil  porous  media  and  clean 
porous  media:  size  straining  may  be  a dominant  mechanism  if  blocking  effect  is  operative 


Release  of  Colloids  jn  jpil 

Different  from  colloid  deposition,  which  is  essentially  determined  by  the 
interactions  at  larger  separation  distance,  colloid  release  depends  on  interactions  at 
separation  between  the  primary  minimum  (<hvnm)  and  surface  (Figure  2-1)  and  is  greatly 
influenced  by  additional  repulsion  at  short  separation.  Colloid  release  is  generally 
expected  when  the  repulsion  is  generated  between  porous  media  and  colloids.  The 


efficiency  of  colloid  release  depends  on  whether  the  process  is  caused  by  diffusion  alone 
or  by  an  applied  external  force  (Kallay  et  al„  1987).  Without  external  forces,  the  rate  of 
particle  release  is  proportional  to  the  diffusional  escape  probability  of  the  deposited 
colloids  through  the  energy  hairier,  i.c.,  the  energy  well  Mw.-rfnm).  It  is  otherwise 
determined  by  the  shape  of  interaction  energy  profile  between  surface  and  the  primary 
energy  minimum  (Figure  2-1 ) as  well:  the  smaller  the  slope  (di^dh)  of  the  curve  the  less 
force  needed  for  colloid  release,  where  h is  the  separation  distance.  Up  to  this  point,  it  is 
clear  that  the  two  opposite  processes,  colloid  deposition  and  release,  are  controlled  by 
different  mechanisms  and  factors.  For  a given  system,  mobile  colloid  concentration  is 
determined  by  the  relative  magnitude  of  colloid  deposition  and  release  rates.  The 
residence  time  of  colloid-carrying  solution  in  the  system  also  plays  an  important  role.  A 
longer  residence  time  represents  a condition  closer  to  the  equilibrium  of  colloid 
deposition  and  release,  while  a short  one  may  discriminate  deposition  over  release  and 
vice  versa.  For  example,  colloid  dispeision  in  batch  experiments,  in  general,  is  a result  of 
balancing  between  colloid  deposition  and  release  because  of  the  longer  residence  time  in 
general;  whereas  colloid  release  from  a short  column  may  be  free  of  redeposition  effects. 
Therefore,  one  must  take  extra  caution  when  extrapolating  colloid  stability  obtained  from 
botch  experiments  to  colloid  transportability  in  column  experiments  where  kinetic  effects 
are  more  pronounced. 

In  a soil,  there  are  numerous  energy  profiles  for  various  colloids  depicted  in 
Figure  2-1  due  to  the  heterogeneity  of  colloids  and  soil  porous  media.  Each  energy 
profile  represents  the  interaction  between  a individual  particle  at  a specific  location  on 
porous  media.  For  simplicity,  we  imagine  that  those  curves  can  be  lumped  into  several 


groups  based  on  similarity,  each  of  those  lumped  profiles  representing  the  interaction  of  a 


energy  maximum  and  minimum.  Therefore,  the  deposition  and  release  of  soil  colloids 
may  manifest  themselves  in  a strong  chromatographic  manner  (Chapter  5). 

The  processes  of  colloid  release 

Similar  to  colloid  deposition,  colloid  release  rate  is  generally  determined  by  both 
colloid  detachment  from  porous  media  and  transport  to  bulk  solution.  Ryan  et  al  (1994a) 
examined  the  release  rates  of  colloidal  hematite  from  hematite-coated  sand  columns 
under  different  ionic  strengths  and  flow  rates.  They  found  that  colloid  release  rate 
decreases  as  ionic  strength  increases,  which  was  supposed  to  reduce  the  size  of  the 
energy  barriers.  In  addition,  colloid  release  rate  also  decreases  with  increasing  flow  rate, 
which  is  contrary  to  the  expectation  that  a greater  mobilization  would  occur  at  a greater 
flow  rate  due  to  greater  hydraulic  stress  on  deposited  colloids.  They  suggested  the  rate- 
limiting  step  in  colloid  release  is  the  transport  of  detached  colloids  to  the  bulk  fluid  when 
rapid  colloid  release  corresponds  to  conditions  where  the  energy  barrier  has  vanished 
from  the  potential  energy  profile.  Coincidentally,  a similar  phenomenon  has  been 
observed  by  Jacobsen  et  al.  (1997)  in  intact  soil-column  experiments  where  natural 
particles  are  leached  with  tap  water.  Their  results  showed  particle  mobilization  is  not 
influenced  by  increase  in  flow  rale.  Further,  the  plot  of  accumulated  amount  of 
mobilized  particles  versus  square  root  of  time  shows  a fairly  linear  relation,  implying 
diffusion  limited  kinetics.  One  necessity  for  diffusion  kinetics  being  dominant  is  the 
barrier  maximum  of  mobile  colloids  vanishes,  which  may  be  very  common  in  soil 
because  of  the  transient  phenomenon  discussed  before.  The  lack  of  flow  rate  effect  in  the 


latter  study  may  be  caused  by  the  smaller  slope  of  energy  profile  (Figure  2-1)  for 
potential  mobile  colloids  at  the  given  condition,  which  may  be  overcome  by  increases  in 
hydraulic  stress,  subsequently,  these  potential  mobile  colloids  become  truly  mobile. 
Therefore,  colloid  dilution  in  effluent  may  be  compensated  by  its  enrichment  from 
mobilizing  more  potential  mobile  colloids  by  increased  flow  rates.  This  is  based  on  the 
heterogeneity  of  the  interactions  between  colloids  and  media  in  soil. 

Mobile  colloids  and  water-dispersible  clay 

Kaplan  et  al  (1993;  1994;  1997)  examined  extensively  the  release  of  colloids  &om 
an  Ultisol  in  the  southeastern  coast  of  the  US.  Mobile  colloids  were  collected  during  and 
after  mild  rain  events  from  the  reconstructed  pedons.  They  found  the  mobile  colloids  are 
similar  to  the  water-dispersible  clay  in  mineralogy,  which  is  consistent  with  the  results 
reported  by  Seta  eta).  (1997).  However,  there  were  some  differences.  Kaplan  etal 
(1997)  found  that  the  sizes  and  compositions  of  the  mobile  colloids  differed  from  those  of 
the  water  dispersive  clay,  with  tile  former  being  smaller  and  generally  enriched  with 
kaolinite,  Fc  oxides,  gibbsite,  and  organic  carbon.  In  addition,  compared  to  the  total  clay 
fractions  in  the  reconstructed  pedon,  the  mobile  colloids  are  more  dilute  in  quartz  and 
HIV  (hydroxy-imerlaycred-vermiculite).  Based  on  the  results  from  scanning  electron 
microscope  (SEM)  and  photon-correlation  spectroscopy,  essentially  all  the  mobile 
colloids  (>90  %)  have  diameters  of  about  0.23  pm  and  moved  through  the  soil  as 
discrete,  non-aggregated  particles.  They  concluded  the  particles  enriched  with  mobile 
colloids  are  not  only  readily  dispeisible  but  also  smaller  in  size  than  water  dispersive 
clay.  Further,  they  proposed  colloid  mobilization  in  the  pedons  is  a result  of  two 
consecutive  processes:  dispersion  of  highly  charged  particles  due  to  changes  in  soil 
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chemislry  and  induced  water  flow,  which  is  determined  by  particle  composition  and  can 
be  evaluated  by  water  dispersibility  of  a soil,  and  their  transport  through  a pedon,  which 
is  size-dependent.  However,  it  has  to  be  realized  that  water  dispersibility  of  soil  particles 
obtained  from  a batch  experiment  is  not  always  equivalent  to  colloid's  ability  to  stay 
detached  (mobilized)  in  a soil.  Colloid  residence  time  has  to  be  taken  into  account, 
which  has  been  discussed  above  and  is  supported  by  Kaplan  ( 1 996). 

Influences  of  exchangeable  sodium  percentage  (ESP)  on  stability  of  water-dispersible 
clay  and  mobilization  of  soil  colloids 

Kaplan  (1996)  summarized  if  particle  straining  is  not  a limiting  factor  for  colloid 
release,  the  dispensability  of  soil  water-dispersible  clay  may  correlate  to  colloid 
transportability  in  soils.  The  dispersability  of  water-dispcrsable  clay  has  been  studied 
extensively  (Miller  et  al,  1990;  Sumner  et  al„  1993).  Even  though  we  cannot  completely 
understand  the  dispersability  of  soil  colloids,  it  is  undoubtedly  affected  by  factors  such  as 
pH,  ionic  strength  and  exchangeable  sodium  percentage. 

Kaplan  et  al.  (1997)  reported  that  colloid  concentrations  in  the  effluents  from  soil 
pedons  arc  highly  correlated  to  the  ESP  of  those  soils  (R!=0.9 1 -0.93,  p<  0.01 ),  but  not  so 
to  soil  pH  or  total  electrolyte  concentrations.  They  explained  the  dispersive  property  of 
the  highly  hydrated  Na'  ion  attributes  to  the  positive  correlation  as  the  dispersion  of  the 
soil  colloids  appears  to  be  a more  important  process  in  colloid  release  in  the  sandy  soils. 
In  fact,  the  positive  correlation  between  soil  colloid  dispersability  and  exchangeable 
sodium  percentage  has  been  found  by  other  researchers  (Frenkle  et  al.,  1978;  Suarez, 

1 984)  who  investigated  the  effects  of  irrigation  or  rainfall  on  hydraulic  conductivity  in  a 


sodic  soil. 


Seta  and  Karathanasis  (1997)  examined  the  relation  between  the  stability  and 
tmnsportability  of  water-dispersible  soil  colloids  by  pumping  water-dispersible  clay 
through  an  intact  soil  column.  They  found  colloid  recovery,  after  five  pore  volumes, 
depends  on  not  only  the  type  of  soil  column  but  also  the  characteristics  of  colloids. 

Among  the  colloid  properties,  pH  and  total  exchangeable  bases  are  significantly 
correlated  with  colloid  recovery.  The  effect  of  pH  on  colloid  migration  can  be  explained 
by  its  effect  on  colloid  stability:  the  higher  the  pH  above  the  pHo  (<  4.0),  the  greater  the 
colloid  stability.  The  strong  correlation  between  total  exchangeable  bases  and  colloid 
recovery  is  anticipated  because  of  its  high  correlation  with  pH.  In  addition,  they  stated 
"contribution  of  total  exchangeable  bases  to  colloid  transport  may  derive  from  cation- 
exchange  reaction  between  the  colloid-saturating  cations  and  the  column  matrix,  which 
reduces  colloid  interaction  with  matrix  surfaces  and  thus  enhances  colloid  transport' '. 
Finally,  they  concluded  the  best  single  independent  variable  predicting  colloid  recovery  is 
total  exchangeable  bases  (RJ=0-97). 

Even  though  there  is  some  inconsistency  among  the  results  in  different  studies 
described  above,  there  is  sufficient  information  to  show  ESP  is  consistently  correlated  to 
the  dispersability  of  soil  colloids  in  a soil.  Rcngasamy  (1982),  however,  demonstrated 
that  even  Ca-saturated  clay,  i.e„  ESP=0,  could  be  dispersed  provided  the  soil  is  free  from 
electrolyte  by  dialysis,  suggesting  exchangeable  Na  is  not  essential  for  dispersing  a soil. 
Further,  in  a soil  system,  solution  chemistry'  is  dynamic  and.  therefore,  concentrations  and 
types  of  cations  in  soil  solution  and  exchange  sites  may  vary  greatly.  In  this  sense, 
parameters  such  as  total  exchangeable  bases,  which  takes  cations  besides  sodium  into 
account,  may  be  more  reasonable  to  describe  colloid  transportability.  The  effects  of 
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cation  exchange  reactions  on  colloid  release  have  been  demonstrated  by  Roy  el  al. 
(1996b),  which  will  be-discusscd  in  the  following  section. 

Ion  transfer  processes  durine  colloid  release 

Observation  of  ion  transfer  during  colloid  release  in  wcILdefincd  systems.  Kallay 
et  al.  (1986)  examined  the  effects  of  neutral  electrolytes  on  the  detachment  of  spherical 
colloidal  particles  of  goethite  from  glass  surface  in  basic  media.  They  unexpectedly 
found  colloid  removal  was  enhanced  with  increasing  concentration  of  NaNOj.  Such 
unusual  behavior  could  be  predicted  if  constant  potential  is  assumed.  They  explained  this 
behavior  corresponds  to  a double  layer  that  undergoes  relaxation  during  particle 
detachment.  In  order  to  keep  the  potential  constant  as  the  distance  between  the  particle 
and  the  surface  increases,  the  adsorption  of  the  potential  determining  ion  (OH')  must 
occur  with  consequent  equilibration  of  all  ionic  species  in  the  Stem  and  diffuse  layers. 
During  particle  detachment,  ion  transport  from  solution  phase  into  the  intcrfacial  layer  is 
accelerated  with  increasing  distance  between  the  surfaces  due  to  solution  influx  from 
bulk.  Kallay  clal.  (1987)  summarized  this  unusual  behavior  can  be  observed  under  the 
following  conditions:  colloid  redeposition  is  minimal  (short  column);  particles  and  media 
have  alike  charges;  surface  potentials  are  constant  during  colloid  detachment.  They 
argued  this  apparently  unusual  behavior  can  be  understood  if  one  takes  into  consideration 
the  energy  profile  as  a function  of  distance  near  the  surfaces  (Figure  2-1).  The 
probability  of  colloid  detachment  depends  on  the  depth  of  energy  well  which 

is  a function  of  several  parameters.  On  the  other  hand,  the  probability  of  colloid . 
deposition  depends  on  the  height  of  the  energy  maximum  W™).  As  the  concentration  of 


NaNOj  in  the  system  increases,  the  depth  of  energy  well  reduces,  resulting  in  an 
increased  rate  of  colloid  detachment;  at  the  same  time,  the  energy  maximum  decreases, 
resulting  in  an  increased  rate  of  deposition.  The  amount  of  colloidal  particles  observed  in 
effluents  is  determined  by  the  net  rate  of  these  two  opposite  processes. 

The  essence  ofKallays  ( 1 986, 1987)  rationale  is  that  a change  in  the  electrical 
environment  of  the  interface,  resulting  from  colloid  detachment,  may  induce  ion  flow 
between  solution  and  interface  (EDL),  which  in  turn  may  facilitate  colloid  detachment. 
This  has  been  known  as  so-called  surface  charge-  or  potential  regulation,  which  has  been 
discussed  by  Chan  et  al.  (1974).  They  have  shown  there  may  be  significant  ion  flow 
between  solution  and  interface  during  EDL  interaction  if  the  particles  carry  pH  dependent 

since  pH-dcpcndcnl  charge  is  generally  abundant  on  soil  colloid  particles. 

Ion  transfer  during  soil  colloid  mobilization  in  literature 

Shainberg  ct  al.  (1981a)  examined  the  effects  of  electrolyte  concentration  on 
hydraulic  conductivity  (HC)  of  a sodic  soil.  One  of  the  major  mechanisms  of  the 
influence  of  electrolyte  concentrations  on  hydraulic  conductivity  is  by  mobilizing  soil 
colloids  and  then  plugging  them  into  conducting  pores.  They  found  both  clay 
disperaibility  and  hydraulic  conductivity  of  the  soil  were  very  sensitive  to  the  level  of 
exchangeable  Na  in  the  soil  and  to  the  salt  concentration  of  the  percolating  solution. 

and  HC  decreased  only  if  ESP>12%.  Conversely,  when  salt  concentration  was 
maintained  at  - 0.5  meq/liter,  clay  dispersion  increased  and  HC  decreased  when  ESP  >1- 
2%.  These  results  indicate  that  ESP  itself  can  not  determine  soil  dispersibility,  and  salt 
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concentrations  in  solution  have  to  be  taken  into  account.  A simple  principle  behind  this 
is  colloids  flocculate  when  the  critical  flocculation  concentration  of  u salt  is  reached.  The 
higher  the  ESP  of  the  clay  and  the  lower  the  salt  concentration  in  solution,  the  higher  the 
tendency  for  a soil  to  disperse,  and  consequently  the  HC  decreases.  Further,  they 
concluded  die  response  of  soils  to  low  ESP  and  leaching  with  low  electrolyte  water 
depends  on  the  concentration  of  electrolytes  in  the  soil  solution  that  the  solid  phases  of 
each  soil  maintains  (Shainberg  et  al„  1981b).  They  demonstrated  salt  concentration  in 
solution  was  determined  by  tire  dissolution  rate  of  soil  minerals,  and  higher  mineral 
dissolution  rate  resulted  in  less  impact  of  exchangeable  Na  on  soil  dispersion.  Obviously, 
soil  dispersion  results  from  equilibrium  of  colloid  deposition  and  release.  However,  it  is 
clear  that  the  interaction  between  cations  on  exchange  sites  and  bulk  solution 
significantly  impacts  colloid  mobilization. 

There  is  sufficient  evidence  exchangeable  Na  has  significant  influence  on  colloid 
release  from  soil  when  it  is  exposed  to  an  influent  with  low  ionic  strength.  Cummins  and 
Kelley  (1923)  first  demonstrated  the  hydrolysis  of  exchangeable  Na’  in  soil,  or  the 
replacement  ol  exchangeable  Na'  by  H*  from  the  dissociation  of  water.  They  found,  in 
the  absence  of  CaCOi  and  COj.  a Na-saturated  soil  leached  with  distilled  water  yields  a 
NaOH  solution.  In  an  experiment  with  Na-montmorillonite,  where  the  hydrolysis 
products  are  removed  continuously.  Bar-on  and  Shainberg  (1970)  found  a Na* 
concentration  of  0. 1 W in  the  effluent.  Shainberg  (1973)  demonstrated  Na- 
montmorillonite  releases  Na'  even  when  the  reaction  products  are  not  removed.  He 
found  the  specific  conductance  of  the  clay  suspension  is  proportional  to  the  square  root  of 
time.  These  observations  are  consistent  with  a hydrolysis  mechanism  consisting  of  two 
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consecutive  reactions:  a rapid  exchange  between  exchangeable  Na*  and  H*.  which  results 
in  an  acidic  surface;  a slower,  first-order  transformation  of  H clay  to  Mg  or  A1  clay, 
which  increases  the  amount  of  exchangeable  Na*  release.  In  soils,  similar  phenomena 
have  been  observed.  Oster  and  Shainbcrg  (1979)  demonstrated  that,  washing  three  arid- 
zone  soils  with  distilled  water,  the  release  of  electrolytes  can  be  related  to  the  square  root 
of  time  exhibiting  two  linear  rates  in  time  sequence.  They  concluded  that  the  first  rate, 
the  more  rapid  of  the  two,  which  occurred  right  after  the  washing  (<1  -2  h),  depends  on 
exchangeable  Na.  Increasing  ESP  causes  increases  in  the  release  rate  in  electrolytes. 

In  general,  the  flux  of  each  ion  species  i across  the  boundary  of  the  EDL  around  a 
panicle  is  governed  by  the  convective-diffusion  equation  (Van  der  ven,  1989): 

Ni  = -Zi  ui  F Ci  Vip  -D,Vc,  +c,  v (2-7) 

Where  Ni , Vip,  Vc,  and  v are  the  flux  of  ion  species  i,  electric  field  acting  on  i, 
concentration  gradient  of  ion  species  i from  bulk  solution  to  interface  and  the  bulk 
velocity  of  the  fluid  motion,  respectively.  They  are  all  vector  quantities.  Zj  F.  u, , c,  and 
D,  arc  the  charge  per  mole,  mobility,  diffusion  coefficient,  and  bulk  concentration, 
respectively,  of  ion  species  i.  Once  a colloidal  particle  starts  to  detach  from  an  electrified 
surface  in  a given  system,  the  electric  field  acting  on  i in  the  surrounding  EDL  will 
change  (d  (Vqrydt)  in  the  course  of  colloid  detachment,  resulting  in  changes  in  N, 
because  Vc,  and  v do  not  change  much  during  colloid  detachment  under  stagnant  flow 
during  a very  short  period.  With  this  in  mind,  the  experimental  observations  we  just 
discussed  above  can  be  understood  conceptually.  Once  the  electric  field  acting  on  a 
colloid  decreases  during  colloid  detachment,  exchangeable  Na  is  more  readily  able  to 
diffuse  to  bulk  solution  compared  to  Co2',  Mg2*,  K*  and  etc.  This  has  been  confirmed  by 


the  fact  thai  electrolyte  concentration  in  solution  increases  much  faster  at  higher  ESP  for 
a soil  (Osier  et  al.,  1 979).  In  essence,  this  is  a process  of  EDL  expansion,  which  results 
in  increases  in  repulsive  forces  between  soil  colloids  or  colloids  and  porous  media,  and 
thus  mobilization  of  soil  colloids.  Therefore,  it  is  easy  to  understand  the  observation  that 
colloid  mobility  are  proportional  to  ESP  (Kaplan,  1 996)  as  well  os  total  exchangeable 
bases  (Seta  et  al.,  1997).  This  is  consistent  with  the  study  by  Roy  et  al.  (1996),  which  is 
based  on  cation  exchange  reactions. 

Ion  transfer  during  soil  colloid  release  in  our  research 

We  have  examined  how  colloid  mobility  in  a Pb-contaminated  soil,  collected 
from  Montreal,  Canada,  is  affected  by  water-flooding  incubation  (Chapter  3).  The  soil 
packed  in  a series  of  short  columns  was  incubated  for  approximately  3, 20,  and  80  d, 
respectively.  After  the  standing  water  on  the  top  of  soil  columns  was  removed,  0.01  M 
CaCI;  solution  was  pumped  through  the  soil  columns  until  the  soil  was  saturated  with 
CaCU-  This  is  designed  to  eliminate  the  artifacts  resulting  from  column  packing  and  to 
use  Ca  as  an  index  for  ion  transfer  during  colloid  release.  The  influent  was  then  switched 
from  CaClj  to  deionized  distilled  water  (DDW)  and  effluent  of -1 1.5  mL  was 
continuously  collected  using  a fraction  collector.  Total  and  dissolved  metal 
concentrations  were  analyzed.  The  relation  among  colloidal  Al,  Fe,  dissolved  Ca,  and  pH 
with  pore  volumes  of  the  effluent  is  presented  in  Figure  3-2  to  demonstrate  the  proposed 
mechanisms  for  Ca  release. 

After  saturating  the  soil  columns  with  0.01  CaCb  solution,  Ca  was  the  dominant 
cation  in  both  bulk  solution  and  electrical  double  layers  surrounding  colloidal  particles  or 
stationary  solid  phases  in  the  soil.  When  the  effluent  is  switched  fiom  0.01  M CaCb 
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solution  to  DDW,  the  release  of  Ca  from  the  soil  is  determined  by  various  mechanisms, 
which  most  likely  comes  to  play  sequentially.  Firstly,  at  the  beginning  stage  of  switching 
to  DDW,  Ca  release  is  from  the  bulk  solution,  which  is  confirmed  by  the  fact  that  Ca 


concentrations  and  pH  in  the  effluent  at  pore  volume  =1  arc  approximately  those  of  the 
influent  ([Ca]=  0.01  M,  pH=6.95)  (Figure  3-3).  Secondly,  Ca  diffusion  from  EDL  to  the 
bulk  solution  driven  by  a concentration  gradient  occurred  from  two  to  17  pore  volumes. 
When  Ca  in  the  bulk  solution  was  depleted,  Ca  cation  moved  against  electrostatic 
attraction  away  from  the  surfaces,  extending  EDL  surrounding  colloids  and  solid  phases 
to  mobilize  colloids.  Obviously,  any  significant  expending  of  EDL  and  resulting 
colloidal  mobilization  have  to  be  associated  with  significant  amount  ofCa  released  into 
bulk  solution.  In  this  case,  the  excess  Ca2*  in  the  diluted  bulk  solution  may  be  combined 
with  hydroxyl  to  release  proton  (Ca2*  + HiO  <->  CaOH'  + H*)  (Oster  et  al.,  1979), 
resulting  in  a decrease  in  pH,  which  is  consistent  with  our  data  (Figure  5-3)  and 
consistent  with  the  mechanism  we  proposed  above.  Generally,  Ca  hydrolysis  is  weak  in 
solution  and . thus,  not  likely  to  bring  the  pH  down  os  much  as  one  unit  as  shown  in 
Figure  5-3.  In  a heterogeneous  system  such  as  this  one,  however,  deficiency  of  negative 
charge  existed  in  the  bulk  solution,  hence  it  is  possible  that  the  charge  deficiency  may 
further  drive  Ca  hydrolysis.  Thirdly,  Ca  desorption  from  colloid  surfaces  and  stationary 
solid  phases  occurred  from  1 8 to  23  pore  volumes  (Figure  5-3).  This  is  a slow  process 


between  colloids  and  stationary  solid  phases.  Increases  in  Ca  concentration  and  pH  were 
observed  in  the  effluent  in  pore  volumes  17-23  (Figure  5-3).  In  essence,  this  process  is 
equivalent  to  sorbed  Ca  displaced  by  H*.  However,  a direct  validation  is  needed. 


Nevertheless,  this  could  be  a mechanistic  illustration  for  ion  transfer  between  EDLs  and 
bulk  solution  as  described  by  Eq  (2-7).  Significant  ion  transfer  (Ca  and  OH")  during 
colloid  release  provides  a more  promising  approach  to  evaluate  colloid  mobility  in  a soil 
since  modern  colloid  theory  has  found  serious  difficulty  in  soils  because  of  its 
heterogeneity  in  both  structure  and  composition  of  soil  particle  surfaces. 

A^gdaliiMJ2LC9.1kld.^t.lLHeavv  Metals 

There  are  two  major  issues  involved  here:  sorption  capacity  of  heavy  metals  to 
soil  colloidal  particles,  which  has  been  well  documented  (Hayes  et  al.,  1990;  Dzombak, 
et  al..  1990)  and  will  be  briefly  described  in  the  following  section,  and  influences  of 
metal  sorption  on  colloid  mobility.  The  association  of  colloids  with  heavy  metals  can  be 
classified  into  adsorption,  precipitation  (surface  precipitation),  and  ion  exchange.  Their 
influences  on  surface  charge  are  possibly  by  virtue  of  heavy  metal  ions  functioning  as 
specific  adsorbing,  potential  determining,  and  indifferent  ions,  respectively.  Since  most 
heavy  metals  are  low  in  solubility  in  soils  and  their  concentrations  are  much  lower  than 
those  of  electrolytes  in  soil  solution,  it  is  rare  for  them  to  function  as  indifferent  ions. 
The  same  reasoning  may  also  apply  to  the  role  heavy  metal  ions  play  as  potential 
determining  ions.  However,  it  has  been  realized  that  adsorption  is  one  of  the  most 
important  mechanism  of  heavy  metal  association  with  colloids  (Hayes  et  al..  1990)  and 
colloid-facilitated  heavy  metal  transport  (Kretzschmar  et  al.,  1997),  which  will  be 
discussed  next. 

In  soils,  colloids  are  heterogeneous,  including  weathered  minerals  (clay,  metal 
oxides,  etc.),  CaCOj,  silica,  large  organic  molecules  and  cellular  debris.  However,  they 


may  be  classified  based  on  surface  fimclional  groups:  surface  hydroxyl,  carbonyl,  organic 
complexation  group,  etc.  As  an  example,  we  will  discuss  hydroxyl  surfaces  in  derail. 

Adsormion  of  heavy  meials  lo  colloids  with  surface  hydroxyl 

Surface-hydroxyl-bcaring  minerals  arc  abundant  in  soil,  such  as  metal  oxides, 
quartz,  and  the  edge  of  knolinite.  Similar  with  the  expression  for  the  de-  and  protonation 
processes  described  previously,  surface  complexation  reactions  of  surface  hydroxyl  with 


heavy  metals  can  be  written  as  follows  (Sturnm  et  al.,  1987): 

SOH  + M'*  = SOMu'"  + H*  K|.w  (2-8) 

2 SOH  + M"=  (SOhM"'11  + 2 H*  Kj.w  (2-9) 

Where  SOH  denotes  a surface  site  and  Mz*  represents  a metal  cation.  Ki.  *,p  and  K,  w arc 
the  apparent  surface  equilibrium  constants  and  defined  as: 

K,.app=  [SOM,'*)  *]  [H*J/[SOH]  [M'*|  (2-10) 

Kj.  ((SO)jM,'-3,](H1!/[SOHJ![M^  (2-1 1) 


where  |]  indicates  concentrations.  Because  of  the  charge  characteristic  of  solid  surfaces, 
ion  activity  at  the  surface  needs  to  be  corrected  to  obtain  the  intrinsic  equilibrium 
constants.  Accordingly,  the  intrinsic  equilibrium  constants  can  be  expressed  as: 

K,.i„=  K i . app  exp((z- 1 )F  iji/RT]  (2-12) 

K;.  K2.  gpp  exp((z-2)Fi|t/RT|  (2-13) 

where  v is  the  potential  difference  between  the  binding  site  and  bulk  solution.  The 
exponential  term  accounts  for  the  coulombic  contribution  to  the  intrinsic  equilibrium 
constants.  Because  the  surface  potential  cannot  be  determined  experimentally,  it  is 
generally  formulated  based  on  a variety  of  models,  such  as  constant  capacitance,  diffuse 
layer  and  triple  layer  models.  It  has  been  found  that  these  models  are  equivalent  in 


predicting  i 


Applying  mass  and  charge  i 


intrinsic  equilibrium  constant  for  metal  ions  binding  to  surface  can  be  determined; 
however,  it  is  somewhat  model-dependent. 

It  has  been  found  stability  constants  of  surface  complexes  correlate  with  those  of 
the  hydroxo  complexes  in  aqueous  phases  (Schindler  et  ah,  1987).  This  could  be 
supporting  evidence  for  Eqs  (2-8)  to  (2-12),  which  assume  a surface  complex  is  an 

intrinsic  and  coulombic  terms.  Recently,  Blesa  ( 1 990, 1 995)  pointed  out  the  importance 
of  ion  hydrolysis  during  adsorption,  and  Sveijensky  ( 1 993)  has  divided  the  intrinsic  terms 

AGM,  = dGii  + AGs  + AGsci  (2-14) 

In  this  equation,  the  coulombic  term  (AGaai)  represents  its  contribution  to  the  overall 
free  energy  of  adsorption  owing  to  the  interaction  between  the  ion  and  surface  charge,  the 
solvation  term  (AGs)  accounts  for  the  role  of  solvation  during  sorption  of  the  ion,  and  the 
remaining  term,  named  an  ion-intrinsic  term  (AG,,),  is  assumed  to  be  a property  of  the  ion 
alone.  With  this  approach,  the  experimental  results  for  Pb2’,  Cd2\  Cu2’,  Mg2'  and  Ca2* 
can  be  closely  reproduced  (Sverjensky,  1993). 

If  ligands  other  than  H2O  and  OH'  exist  in  solution,  such  as  Lewis  acids  or  bases  ( 
often  the  case  in  soil),  competition  may  occur  among  dissolved  ligands,  heavy  metals  and 
surface  adsorption  sites.  The  competing  anions  (e.g.,  sulfate,  phosphate  or  dissolved 
organic  carbon,  DOC)  can  alter  surface  sorption,  and  formation  of  aqueous  metal 
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the  other  hand,  metal  sorption  may  be  enhanced  by  the  formation  of  ternary  surface 
complexes  via  sulfate  and  phosphate  and  DOC.  At  the  present  time,  there  is  no  theory 
that  can  describe  these  complicated  processes  well.  However,  in  principle,  the  partition 

of  any  aqueous  species  may  be  described  under  the  framework  ofEqt4W4)  (Blesa,  1990; 

4* 

1995;  Svetjensky,  1993). 

In  metal  oxide  systems,  heavy  metal  ions  interact  with  surfaces  as  specific 
adsorbed  ions  (Lyklema,  1984;  Ardizzone,  1983)  to  alter  charge  development,  which  is 

Effeptsofheqv^jisyit)risMiMKiB.ffl|.lllP..t;llW|jW..9fj9j|jSll9Wi 

A number  of  studies  have  indicated  the  importance  of  particle-solution 
interacuons  in  determining  surface  properties  of  colloids  in  natural  waters.  Elcctrokinclic 
measurements  made  on  natural  particles  dispersed  in  sea  water  and  natural  fresh  waters 
have  shown  that  their  elecirophoreuc  mobilities  tend  to  fall  within  vety  limited  ranges  of 
negative  values.  The  apparent  uniformity  of  this  surface  electrical  property  has  been 
attributed  to  the  adsorption  of  organic  materials,  particularly  humic  compounds,  at  the 
particle  surfaces,  which  has  been  supported  by  a number  of  studies  using  synthetic 
systems.  It  has  been  shown  that  the  effects  of  adsorbed  organic  compounds  on  surface 
properties  of  minerals  present  in  natural  waters  are  modified,  to  some  extent,  by 


, the  interaction  of  heavy  metals  with  I 


and  pure  minerals,  rather  by  those  between  heavy  metal  and  the  complexcd  surfaces  of 
"real"  soil  colloids.  This  complication  includes  the  irregularity  of  colloid  shapes  and 
three-dimensional  interfaces  between  solution  and  colloids.  One  of  most  important 
examples  is  organic-coaled  oxides.  It  has  been  proposed  as  a model:  individual  colloid 
particles  mostly  consist  of  inner  mineral  cotes  with  sorbed  Fe,  A1  & Mn  (hydroxy )oxides 
and/or  organic  matter  (OM)  (Day  et  al.  1994).  which  generally  possesses  similar 
composition  with  soil  colloids  as  well  as  three  dimensional  interfaces. 

There  has  been  a fair  amount  of  research  done  dealing  with  how  ionic  strength  or 
pH  affects  the  mobility  of  the  colloids  with  the  constituent  components  above  (Day  et  al, 
1994).  However,  little  is  known  about  the  effects  of  heavy  metals  on  colloidal  charge 
development  and  stability  of  the  complex  colloids.  Only  recently,  Kretzschmar  et  al. 
(1997),  have  explicitly  taken  into  account  the  effects  of  adsorption  of  Pc'  and  Cu2*onto 
humic-coated  oxide  colloids  on  colloid  charge  development  and  transportability.  They 
investigated  the  influence  of  adsorbed  heavy  metals  (Cu2*  and  Pb2*)  on  the  transport  and 
deposition  kinetics  of  submicron  size  hematite,  coated  with  humic  acid  suspended  in  Ca2* 
solution,  through  a natural  soil.  They  found  that  replacement  ofCa2*  by  Pb2*  while 
holding  the  total  concentration  (Ca2**  Pb2')  fixed  resulted  in  a slight  decrease  in 
electrophoretic  mobility  of  humic-coated  hematite  colloids.  When  Ca2*  was  completely 
replaced  by  Pb2*,  the  suspensions  were  destabilized  and  aggregated  within  20  h.  In 
contrast,  replacement  of  Ca2'  by  Cu2*  hod  very  little  effect  on  electrophoretic  mobility 
and  colloidal  stability.  Both  Pb2*  and  Cu2*  are  known  to  bind  much  more  strongly  to 
humic  substances  and  Fe  oxide  surfaces  than  Ca2*.  Al  the  pH  of  the  suspensions 
investigated  (pH  5.7),  Cu2*  has  a higher  affinity  for  humic  substances  than  Pb2*,  but  Pb2* 
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adsorbs  more  sirongly  to  hematite  than  Cu2'.  In  fact,  the  interface  is  three-dimensional 
since  the  humic  coating  can  be  considered  an  extension  of  hematite  panicles  to  the 
solution  phase.  The  bonding  difference  may  direct  Pb2'  and  Cu2'  to  locate  preferentially 
at  Stem  and  diffitse  layers,  respectively.  Therefore,  it  is  expected  that  the  effect  is  more 
pronounced  by  replacing  Ca2*  with  Cu2*  than  with  Pb2’  in  the  colloid  suspensions  on 
electrophoretic  mobility  and  colloid  deposition.  Their  results  suggested  that  humic- 
coated  oxide  colloids  can  be  stable  and  mobile  in  the  presence  of  strongly  adsorbing  trace 

Partitioning  of  heavy  metals  in  soil  colloids 

The  ability  of  colloids  to  facilitate  heavy  metal  transport  is  largely  determined  by 
the  distribution  of  heavy  metals  among  solution,  mobile  colloids,  and  stationary  grains. 
Mills  et  al,  (1991)  first  incorporated  partition  coefficients  into  their  model  in  describing 
colloid-facilitated  metal  transport  in  porous  media.  However,  a mechanistic 
understanding  of  metal  distribution  among  different  phases  is  still  missing.  One  of  the 
major  reasons  is  the  dynamic  feature  of  the  distribution  of  heavy  metals,  and  phase 
transformations  among  them,  such  as  from  mobile  to  immobile  colloids,  once  heavy 
metals  are  introduced  to  a system.  Kalplan  et  al.  ( 1 995)  have  found  that  the  mineralogy 
of  mobile  colloids  from  contaminated  wells  is  significantly  different  from  those  from  the 
contaminated  soil  on  one  site,  suggesting  the  association  with  heavy  metals  may  alter 
greatly  the  mobility  of  original  mobile  colloids. 

Amrhein  et  al.  (1993)  examined  the  potential  of  colloid-facilitated  heavy  metals  in 
roadside  soils  receiving  deicing  salts.  There  are  two  deicing  agents  considered,  NaCl  and 
calcium  magnesium  acetate  (CMA).  In  their  experiments,  a series  of  60-mL  syringes 


were  packed  with  30  g of  soil.  The  columns  were  leached  on  a centurion  vacuum 
extractor  with  three  consecutive  30-mL  aliquots  of  either  0.1  mol  L'1  NaCI  orCMA. 
After  initial  leaching  with  either  of  the  salts,  the  columns  were  leached  with  three 
consecutive  30-mL  aliquots  of  deionized  water.  This  is  to  simulate  the  input  of  salty 
ntnolf  water  to  the  roadside  soil,  followed  by  snowmelt  or  rainfall.  A portion  of  each 
leachate  was  saved  without  any  further  nitration  and  the  remainder  of  each  solution  was 
tillered  through  a 0.45  pm  membrane  filter.  Then  tile  filtrates  were  immediately  placed 
into  a stirred  ultrafiltration  cell  with  a membrane  of  1000  MWCO.  This  procedure 
separates  the  particles  into  three  fractions:  >0.45  pm.  between  1 .0  nm  and  0.45  pm  and  < 
1 .0  nm.  The  results  have  shown  that  the  cumulative  leached  metals,  through  the  whole 
process,  vary  with  the  initial  salt  input.  In  general.  NaCI  tends  to  mobilize  soil  colloids 
more  than  CMA  ( Amrhein  et  al.,  1 993).  On  the  other  hand,  because  of  the  possibility  of 
artifacts  resulting  from  column  packing  in  their  study,  it  may  be  more  reasonable  to 
examine  the  leached  cumulative  metals  excluding  the  first  several  pore  volumes.  Figure 
2-2  is  a plot  of  leached  cumulative  metals  after  switching  the  influent  from  salt  solution 
to  deionized  water.  It  is  well  established  that  reduction  of  ionic  strength  may  mobilize 
soil  colloids.  As  shown  in  Figure  2-2,  the  concentrations  of  each  metal  (Cu,  Pb,  Fe,  Ni 
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Figure  2-2.  Cumulative  concentrations  of  metals  leached  from  the  soil  in  the  three  size 
fractions  after  switching  salt  solution  to  deionized  water.  Size  separation  was  by  0.45  pm 
membrane  filter  and  l-nm  cellulose  ultrafillrution  membrane  (adapted  from  Amrhein  et 
al„  1993) 
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and  Cr)  are  consistently  higher  in  the  leachates  from  soil  columns  pre-leachcd  with  NaCI 
than  that  with  CMA  in  all  three  size  fractions.  The  fraction  <1  nm  will  hereafter  not  be 
emphasized  since  it  is  generally  considered  as  dissolved.  Among  those  size  fractions, 
colloids  ranged  from  I -450  nm  is  most  important  for  colloid  transport  (Kaplan  et  al., 

1 995),  and  their  concentrations  of  each  metal  are  much  higher  for  the  columns  pre- 
leached  with  NaCI  than  those  with  CMA.  This  is  consistent  with  normal  intuition,  that  is, 
Na  promotes  colloid  mobility,  and  so  does  colloid-facilitated  metal  transport  if  one 
assumes  association  of  metals  with  mobile  colloids  are  not  significantly  affected  by 
changes  in  solution  chemistry.  Furthermore,  metals  in  fraction  1.0  -450  nm  are  much 
higher  than  those  in  fractions  > 450  nm  for  columns  pre-leached  with  CMA,  whereas 
there  are  no  significant  differences  for  the  column  with  NaCI  except  for  Fc  and  Cu.  This 
suggests  CMA  is  more  selective  in  mobilizing  smaller  colloidal  particle  than  NaCI. 

Along  with  this  fact,  the  authors  pointed  out  that  colloidal  Cu  did  not  correlate  well  with 
DOC  between  the  two  fractions  for  the  columns  pre-leached  with  NaCI,  nor  docs 
colloidal  Fe  with  other  colloidal  metals  pre-leached  with  CMA  (Figure  2-2).  This 
suggests  that  heavy-metal-bearing  mobile  colloids  vary  with  both  size  and  solution 
chemistry.  In  other  words,  phase  transformation  among  heavy-metal-bearing  colloids, 
such  os  mobile  to  immobile  colloids  or  among  mobile  colloids  with  various  sizes,  has  to 
be  taken  into  consideration. 

We  have  examined  the  effects  of  water  flooding  incubation  on  colloidal  Pb 
mobility  in  soil  columns.  Detailed  descriptions  about  the  experiment  procedure  can  be 
found  in  Chapter  5.  As  shown  in  Figure  5-2.  colloidal  Pb  concentrations  varied  with 
incubation.  Interestingly,  there  was  no  significant  difference  in  concentrations  of 


colloidal  Fe  and  A1  in  the  firs!  several  pore  volumes  when  the  columns  were  incubated 
from  3 to  20  d.  while  concentrations  of  colloidal  Pb  decreased  approximately  4 times 
after  20-d  than  3-d  incubation.  This  suggested  that  Pb-bearing  colloids  tended  to  be 
immobilized  or  the  tendency  of  association  of  Pb  with  mobile  colloids  decreased  with 
water-fiooded  incubation. 


Concluding  Remarks 

Soil  has  distinctively  different  features  in  colloid  mobility  compared  with  well- 
defined  porous  media  or  subsurface  systems.  At  the  present  time,  it  is  not  realistic  to 
predict  colloid  mobility  and  colloid-frtcilitatcd  heavy  metal  mobility  in  soil  with  certainty. 

Soil  colloids  arc  exposed  to  solution  where  they  take  part  in  chemical  processes. 
Therefore  it  is  difficult  to  predict  charge  development  of  soil  colloids  since  it  is 
determined  by  the  interaction  of  surfaces  and  solution.  Point  of  zero  charge  (PZC)  of 
mineral  particles  suspended  in  simple  solution  can  be  theoretically  predicted,  while  little 
success  has  been  achieved  in  soil.  In  general,  soil  particles  are  assemblages  of  crystalline 
and  amorphous  minerals,  and  organic  residuals;  for  such  particles  sophisticated  models 
such  os  triple  layer  model  do  not  make  much  sense  in  describing  the  electrified  interfaces 
since  particle  surfaces  arc  irregular  and  the  Stem  layer  is  often  inside  of  the  physical 
boundary  of  the  panicles.  As  a first  approximation,  however,  it  may  be  more  pratical  to 
describ  soil  particle  interfaces  by  diffuse  double  layer  with  a specific  pan  (a  charged 
surface),  and  a generic  pan  (a  diffuse  layer).  The  surface  charge  development  process 
may  be  described  by  the  panitioning  of  charged  species  between  solution  and  surfaces 
assuming  the  charged  species  are  similar  between  the  two  phases.  Similarly,  the 
association  of  heavy  metals  with  colloids  may  have  significant  influences  on  colloid 
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mobility,  which  is  subject  to  changes  in  dynamic  solution  chemistry'-  However,  very 
limited  information  is  available  at  the  present  time. 

Colloids  arc  ubiquitous  in  soils.  Different  from  clean  porous  media,  the  surfaces 
of  soil  stationary  phases  consist  of  deposited  colloids.  Consequently,  the  blocking  effect 
of  colloid  deposition  in  soil  is  obvious,  which  may  greatly  enhance  colloid  mobility  in 
soil  compared  to  that  in  a comparable  clean  porous  medium.  Under  such  a condition, 
size  straining  is  a dominant  mechanism  in  soil  colloid  deposition.  The  deposited  colloids 
may  release  depending  on  the  nature  of  interactions  with  the  stationary  phase,  that  is,  the 
interaction  energy  profiles  (Figure  2-1),  which  vary  because  of  heterogeneity  of  the 
surfaces.  Phenomenally,  this  process  is  time-dependent  or  condition-dependent,  and 
.therefore,  is  referred  to  as  transient  phenomena  in  colloid  deposition.  At  the  present 
time,  water  dispersible  clay  extracted  from  a soil  is  widely  used  to  indicate  colloid 
mobility.  However,  it  has  to  be  realized  that  in  both  coses  there  are  two  simultaneous, 
opposite  processes:  colloid  deposition  and  release.  In  a given  system,  the  relative  effect 
of  each  on  either  soil  clay  dispersion  in  batch  tests  or  colloid  mobility  in  a column 
experiments  depends  on  mobile  colloid  residence  time  in  the  systems. 

Despite  the  shortcomings  of  the  DLVO  (Derjaguin-Landau-Verwey-Overbeek) 
theory,  we  believe  it  provides  a basic  theoretical  framework  in  describing  soil  colloid 
dispersion  and  colloid  deposition  and  release,  which  is  simple  and  applicable  in  concept. 
It  has  to  be  pointed  out  that  all  modem  colloid  theories  (not  only  the  DLVO  theory)  are 
difficult  to  apply  to  soil  because  of  its  heterogeneity  in  both  structure  and  composition  of 
soil  particle  surfaces.  In  order  to  determine  the  potential  mobility  of  soil  colloids,  a 
phenomenological  approach  may  be  more  practical.  Soil  colloids  possess  high 


concentrations  of  exchangeable  ions  in  their  EDLs  (high  CEC),  which  does  not  always 
equilibrate  with  those  in  bulk  solution.  It  has  long  been  recognized  that  colloid  release  in 
soil  is  coupled  by  significant  ion  transfer  between  EDLs  and  bulk  solution.  By 
evaluating  the  possibility  of  the  ion  transfer,  potential  colloid  mobility  may  be  estimated 
qualitatively. 


CHAPTER 3 

RELATION  OF  PB  SOLUBILITY  TO  FE  PARTITIONING  IN  SOILS 

Introduction 

Lead  comoroinalion  in  soils  is  of  cnvironmemal  significance  due  to  its  toxicity  to 
both  humans  and  animals  (Ma  ct  a!.,  1 995).  In  general,  Pb  mobility  is  low  because  of  its 
low  solubility.  Lead  solubility  may  be  further  reduced  as  a result  of  its  interactions  with 
soil  solid  phase  via  sorption  and  ion  exchange.  However,  enhanced  Pb  solubility  has 
been  found  under  both  laboratory  and  field  conditions  (Amrhein,  et  al.,  1994; 
Glazovskaya,  1 994).  As  such,  Pb  may  migrate  through  a soil  profile  to  contaminate 
groundwater. 

Although  much  effort  has  been  spent  to  model  heavy  metal  solubility  (Cederberg 
et  a!.,  1985;  Sposito.  1984),  such  prediction  under  field  conditions  involves  large 
uncertainty.  It  is  partially  because  of  the  difficulty  in  assessing  the  effects  of  dynamic 
soil  solution  chemistry  on  heavy  metal  speciation.  However,  changes  in  solution 
chemistry,  such  as  pH,  redox  potential  and  ionic  strength,  may  shift  Pb  retention 
processes  significantly.  These  impacts  may  be  further  complicated  by  aqueous  Pb 
competition  with  other  cations  for  ligands,  which  may  enhance  Pb  mobility  under  certain 
conditions  (Amrhein.  ct  al.,  1994). 

in  natural  aquatic  environments,  Pb  (II)  and  Fe  (II)  have  similar  affinity  to 
complex  with  ligands,  and  thus  they  show  similar  patterns  of  species  distribution  (Turner 
et  al..  1981).  Theoretically,  both  Pb  (II)  and  Fe  (II)  are  in " transient  class"  in  teims  of 
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metal  classification  based  on  electronegativity  and  covalent  index  (Nicboer  and 
Richardson,  1980).  This  evidence  implies  that  it  is  important  to  examine  Fe  chemistry 
when  studying  Pb  solubility.  In  fact,  good  correlation  (r2  = 0.71-0.91,  p<0.05)  between 
concentrations  of  aqueous  Pb  and  Fc  in  sediments  lias  been  reported  (Lee  et  al..  1997; 
Routh  and  Ikramuddin,  1996).  In  soil  environment,  data  published  by  Karczewska 
( 1 996)  showed  that  concentrations  of  mobile  Fe  and  Pb  were  related  (r2  - 0.54,  p<0.05)  in 
several  polluted  soils  near  a copper  smelter.  The  theoretical  evidence  along  with  the 
published  experimental  data  suggests  that  Pb  and  Fc  solubility  may  be  related  in  a soil. 

In  addition,  there  are  two  other  possible  factors  for  correlation  between  Pb  and  Fe 
in  solution.  First  Fe  (oxy)  hydroxides  are  good  sorbents  for  aqueous  Pb  (Ainsworth  et 
al.,  1994;  Laxen.  1985;  Benjamin  and  Leckie,  1981),  and  their  dissolution  or  precipitation 
may  release  or  sorb  Pb.  In  soils,  dissolution  of  Fe  (oxy)  hydroxides  is  generally 
promoted  by  reducing  Fe  (III)  to  Fe  (II),  which  is  sensitive  to  soil  redox  status  (Lindsay, 
1979;  Gotohand  Patrick,  1974).  Under  oxidizing  condition,  Fe  (oxy)  hydroxides  tend  to 
immobilize  Pb  (Gambrell  et  al.,  1980);  whereas  under  reducing  condition  they  dissolve  to 
release  Pb  (Gambrell.  1994).  Secondly,  soil  organic  carbon,  especially  dissolved  organic 
C (DOC),  is  an  important  factor  controlling  Pb  and  Fe  solubility  in  soil  (Dorr  and 
Munnich,  1991;  Davis.  1 984;  Laxen,  1 98S).  DOC  functions  as  a ligand  to  complex  with 
Pb  and  Fe  to  increase  their  solubility  (Davis  and  Leckie,  1978).  On  the  other  hand,  DOC 
may  also  be  sorbed  by  Fe  (oxy)  hydroxides  under  certain  conditions,  via  ligand  exchange 
or  hydrophobic  interaction  in  soils  (Murphy  and  Zachara.  1 995),  thus  decreasing  Pb  and 
Fc  solubility. 


Despile  all  ihe  circumstantial  evidence  that  Fe  solubility  impacts  Pb  solubility  in 
soil,  there  has  been  little  effort  to  quantify  the  relationship  between  the  two.  However, 
understanding  the  relationship  between  heavy  metal  solubility  and  the  behavior  of  Fe 
may  provide  important  information  for  assessing  potential  mobility  of  heavy  metals  in 
soil  environments,  which  cannot  be  modeled  successfully  at  present  time. 

Soil  redox  status  varies  temporally.  The  intensity  of  redox  status  can  be  described 
by  redox  potential  in  general.  In  surface  soil  it  is  influenced  by  rainfall,  bioactivity,  and 
changes  in  land  use,  whereas  in  vadose  zone  by  fluctuation  in  water  table  (Buol  etal., 

1 997).  The  redox  potentials  of  sediments  that  remain  saturated  can  be  measured  readily, 
whereas  characterization  of  the  redox  status  in  a soil  is  still  a challenge.  As  mentioned 
previously,  however,  soil  redox  status  affects  metal  solubility  greatly  under  certain 
conditions.  To  evaluate  the  effects  of  varying  redox  status,  saturated/unsaturated 
(aerobic/  anaerobic)  incubation  techniques  have  recently  been  developed  (Amrhcin  ct  al., 
1 994;  Karczewska,  1 996).  Water-flooded  incubation  of  soil  was  used  to  effectively 
decrease  soil  redox  potential  to  study  the  spccialion  and  fate  of  heavy  metals  in 
contaminated  soils  (Karczewska.  1996).  Even  though  saturation  rarely  happens  in 
surface  soils,  some  valuable  information  may  be  obtained  from  such  a study.  In  this 
paper,  two  redox  statuses  in  soils  will  be  obtained  under  water-flooding  (anaerobic)  and 
no  flooding  (aerobic)  incubations. 

The  major  objective  of  this  chapter  is  to  examine  the  relationship  between  Pb 
solubility  and  Fe  partitioning  in  soils  with  different  solution  chemistry.  Throughout  this 
paper,  Fe  partition  will  be  expressed  os  a concentration  ratio  of  aqueous  to  sorbed  Fe  (U) 
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Locanon  and  characteristics  of  soil  sample 

The  soil  samples  used  for  this  study  were  collected  in  March  1996  from  a research 
site  in  Hawthorne,  Florida.  The  soil  is  an  acidic  fine  sand  (typic  quartzipsamment)  with  a 
spodic  horizon  below  2 m.  The  samples  were  collected  from  2-20  cm  below  the  surface 
after  removing  organic  residue.  The  samples  were  air-dried,  sieved  through  2-mm 


Column  experiment 

A series  of  60-mL  syringe  (columns,  d=2.6  cm),  which  were  prepacked  with  a 
layer  (~  3mm)  of  acid-washed  Ottawa  sand  (20-30  mesh)  in  the  bottom,  were  packed 
with  40  g of  soil  sample.  The  columns  were  then  packed  with  another  layer  of  the  acid- 

syringes  were  set  on  a Centurion  Vacuum  Extractor  (Centurion  International,  Lincoln, 
NE)  for  incubation  and  leaching.  The  experiment  procedure  was  as  follows: 

Prewetting  soil  columns.  30  mL  of  deionized  distilled  water  (DDW)  was  slowly 


poured  in 


I for  48  h befo 
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Tabic  3-1 . Characteristics  of  the  Florida  soil  used 


PH 

“IS?1' 

Extractable  Fc* 

Texture*  •• 

A1  1 Fe  1 Ca  i Na  | Pb 

Sand  1 Silt  1 Clay 

% 

5.4 

2067  | 976  | 87  1 27  | 1.3 

94.6  | 4.5  | 0.9 

•Heron,  eta!.,  1994. 
••Maetal..  1996. 
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Addition  of  aqueous  Pb.  30  mL  of  0.48  mM  Pb(NOj)2  was  added  to  12  soil 
columns  and  30  mL  of  3.68  raM  PbfNOj);  to  16  soil  columns.  The  solutions  were 
extracted  allcr  24  h at  2.5  mL  h'1.  This  resulted  in  lead-loading  rates  of  0.36  (low  Pb 
loading)  and  2.90  mmoles  kg'1  (high  Pb  loading)  in  two  groups  of  soil  columns. 

Addition  of  electrolyte  solutions  30  mL  of  different  electrolyte  solutions  was 
added  to  different  columns  to  vary  solution  chemistry.  DDW  of  pH  5.5  and  2.17  and 
4.35  mMNaCI  solution  were  added  to  the  first  12  soil  columns,  and  0.52, 1.08,  and  2.17 
mM  NaCI  solutions,  0.25, 0.55,  and  1.10  mM  CaClj  solutions,  and  DDW  of  pH  4,5  and 

5.5  to  the  other  16  soil  columns.  The  solution  was  extracted  out  after  24  hat  the  rte  of 

2.5  mL  h1 . For  the  first  1 2 soil  columns  with  low  Pb  loading,  each  treatment  was 
replicated  four  times  (three  treatments)  and  for  the  16  soil  columns  with  high  Pb  loading. 

Incubation  of  soil  columns.  The  above  28  soil  columns  (with  the  same  Pb  loading 

each)  for  incubation.  One  subgroup  was  filled  with  30  mL  of  DDW  (2  cm  above  the  soil 
surface)  for  water-flooded  incubation.  The  other  subgroup  was  incubated  as  it  is  for  non- 
water-floodcd  incubation.  All  soil  columns  were  incubated  for  40  d at  room  temperature. 

Leaching.  DDW  was  added  to  each  soil  column  to  make  30  mL  of  standing  water 
above  the  soil  for  all  columns  before  leaching.  The  leaching  was  conducted  at  60  mL  h*1 
and  approximately  2S-50  mL  leachates  were  collected.  The  soil  samples  with  moisture 


-25%  in 


efrigerator  for  further  analysis. 


Sample  separation  and  analysis 


The  pH  was  measured  immediately  after  (he  leachates  were  collected.  The  leachates 
were  then  filtered  through  0.22-pm  membrane  filters.  The  filtrate  was  analyzed  for  DOC 
and  total  Pb  and  Fe  concentrations.  The  frozen  soil  columns  were  immersed  in  warm 
water  to  push  the  soil  out  of  the  syringes  with  minimum  disturbance.  The  intact  soils 
from  the  columns  were  then  divided  into  three  equal  sections  (top.  middle  and  bottom) 
for  fitrther  analysis.  The  data  presented  in  this  paper  for  soil,  however,  was  the  average 
over  the  three  sections  since  there  were  no  significant  differences  in  Fe  concentrations  in 

method  described  in  Chapter  3 


Result  and  Discussion 

Pb  solubility 

Aqueous  Pb  concentrations  in  leachates  showed  strong  pH-depcndence  with  r= 
0.92  for  soils  with  Pb  loading  of2.90mM  kg'1  (high  Pb  loading.  Figure  3- 1 A).  Aqueous 
Pb  concentrations  decreased  as  pH  increased  from  4 to  6.  It  has  long  been  recognized 
that  pH  is  an  important  factor  affecting  metal  solubility,  with  aqueous  metal 
concentration  increasing  as  pH  decreases  (Chuan  et  al.,  1996).  Leachate  Pb 
concentrations  in  soils  with  Pb  loading  of  0.36  mM  kg"1  (low  Pb  loading)  were  also 
correlated  with  pH  with  correlation  coefficient  of  0.51  QK0.05).  However,  they 
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Fc  concentration  in  leachates  (pM)  concentrations  (kgL'1)  * 


Figure  3-1  Relation  between  Pb  solubility  and  leachate  pH,  dissolved  organic  carbon 
(DOC),  leachate  Fe  concentration,  and  ratio  of  aqueous  Fe  to  sorbed  Fe(II)  concentrations 
in  a sandy  soil.  Data  from  2.90  mmole  kg*1  Pb  loading  rate  read  from  left  y-axis  and  data 
from  0.36  mmole  kg'1  Pb  loading  rale  read  from  right  y-axis. 
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increased  as  pH  increased  from  6.1  lo  6.8  (Figure  3-1  A).  Obviously,  mechanism's  lhat 
determined  aqueous  Pb  concentrations  in  soils  of  low  Pb  loading  was  different  from  that 
of  high  Pb  loading.  It  seemed  that  DOC  was  a primary  factor  in  determining  Pb 
concentrations  in  soils  with  low  Pb  loading  as  reflected  by  their  correlation  coefficient  of 
1*=  0.47,  p<0.05  (Figure  3-1 B),  i.e.  Pb  concentrations  increased  as  DOC  concentrations 
increased.  However,  such  correlation  did  not  exist  for  soils  with  high  Pb  loading  (r*= 
0.24)  possibly  due  to  the  fact  that  Pb  solubility  was  more  strongly  correlated  to  pH. 

Beside  pH,  redox  status  is  another  factor  affecting  metal  solubility  in  soil  (Chuan 
ct  al.,  1996).  Stumm  (1984)  suggested  that  it  is  reasonable  to  measure  important  redox 
species  instead  of  redox  potential  to  indicate  redox  status  in  natural  systems.  We  used 
aqueous  Fe  as  a measure  of  redox  status  for  soils.  As  expected,  aqueous  Fe 
concentrations  were  much  higher  in  soils  under  water-flooded  incubation  than  non-water- 
flooded  incubation  (Figure  3-1C).  In  soils  with  low  Pb  loading,  aqueous  Fe 
concentrations  was  over  three  times  greater  under  water-flooded  incubation  (averaged 
2.12  nM)  than  under  non-water-flooded  incubation  (averaged  0.57  nM).  Similarly,  the 
average  aqueous  Fe  concentrations  in  soils  with  high  Pb  loading  were  4.46  and  1.74  nM. 
respectively.  Greater  aqueous  Fe  concentrations  in  soils  with  high  Pb  loading  than  with 
low  Pb  loading  was  mainly  due  to  pH  difference  in  the  soils  (Figure  3-1  A).  However, 
there  was  no  relationship  shown  between  aqueous  concentrations  ofFeand  Pb  for  soils  of 
either  Pb  loading  rates  (Figure  3-1C),  It  was  reported  that  aqueous  Pb  and  Fe 
concentrations  are  linearly  related  due  to  the  fact  that  heavy  metals  arc  released  from  Fe 
(oxy)hydroxide  surface  when  it  dissolves  (Chuan  et  al.,  1996).  In  our  study,  however, 
this  was  not  the  case  (Figure  3- 1C),  suggesting  that  aqueous  Pb  may  not  be  controlled  by 


Fe  (oxy)hydroxide  surface  alone.  Similar  results  have  also  been  reported  by  other 
researchers  (Gambrell  et  al„  1991).  These  data  suggested  that,  in  general,  aqueous  Fe 
concentration  alone  could  not  describe  Pb  solubility  well.  Nevertheless,  as  expected, 
aqueous  concentrations  in  soils  under  reduced  conditions  were  much  higher  than  those 
under  less  reduced  conditions.  Aqueous  Pb  concentrations  in  soils  with  high  Pb  loading 
were  52  times  greater  under  water-flooded  incubation  (averaged  490  pM)  than  under 
non-water-flooded  incubation  (averaged  9.4  pM)  (Figure  3-1).  Similarly,  aqueous  Pb 
concentrations  were  91  and  22  nM.  respectively,  for  soils  with  low  Pb  loading. 

Pb  solubility  and  Fe  partitioning 

As  discussed  previously,  pH  was  related  to  Pb  solubility  in  our  study.  However, 
even  for  the  same  soils,  the  relation  between  pH  and  aqueous  Pb  differed  with  Pb  loading 
rates  (Figure  3-1  A).  In  other  wotds,  there  is  no  simple  relation  between  pH  and  Pb 
solubility  when  soil  solution  chemistry  varies  greatly.  There  were  no  simple 
relationships  between  Pb  and  DOC  concentrations  and  Pb  and  Fe  concentrations  either 
(Figure3-IA,C) 

To  better  describe  Pb  solubility,  we  defined  the  ratio  of  aqueous  Fe  to  sorbed  Fe 
(II)  concentrations  as  a Fe  partitioning  index.  In  this  study,  sorbed  Fe  (II)  was 
operationally  defined  as  Fe(II)  that  was  extracted  by  0.5  M HCI  (Lovely  and  Philip,  1 987; 
Heron  ct  al„  1 994).  By  definition.  Fe  partitioning  index  represents  the  relative  affinity  of 
aqueous  Fe  to  solid  phases  in  a soil.  A larger  number  indicates  a lower  Fe  affinity  to 
solid  phases.  There  was  no  simple  relationship  that  existed  between  aqueous  Pb 
concentrations  and  Fe  partitioning  index  (Figure  3-ID).  For  soils  with  low  Pb  loading 
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and  incubated  under  water-flooded  condition,  the  relation  between  aqueous  Fb  and  Fc 
partitioning  index  was  approximately  parabolic.  In  general,  aqueous  Pb  concentrations 
increased  as  the  index  decreased.  However;  as  the  index  decreased  to  below  2 kg  L' , 
increases  in  Pb  concentrations  were  substantial  (Figure  3-ID).  A similar  trend  was  also 
observed  for  soils  with  high  Pb  loading  and  incubated  under  non-water-flooded 
condition,  i.e.  substantial  increase  in  Pb  concentrations  were  observed  when  the  index 
was  < 2 kg  L'1.  However,  no  enhanced  Pb  solubility  was  found  in  soils  with  low  Pb 
loading  and  incubated  under  non-water-flooded  condition  since  the  Fe  index  was  > 2 kg 
L*1.  However,  for  soils  with  high  Pb  loading  rate  and  incubated  under  water-flooded 
condition,  extremely  high  Pb  concentrations  were  observed  since  the  Fe  index  was  <2  kg 
L'1.  Obviously,  compared  to  pH,  DOC  and  Fe  concentrations,  a more  consistent  relation 
between  aqueous  Pb  concentrations  and  the  Fc  partitioning  index  in  soils  was  observed. 

For  a pH-dependent  process,  Kurbatov  equation  (Kurbatov  et  al„  1 95 1 ) has  been 
successfully  used  to  describe  the  macroscopic  behaviors  of  metal  partitioning  in  natural 
aquatic  systems  (Fuller  et  al..  1996;  Tcssieret  al..  1985,  Balistrieri  el  al.,  1983).  In  our 
system,  the  Fc  index  in  soils  with  both  high  ([*=0.61,  p<0.05)  and  low  ([*=0.55,  p<0.05) 
Pb  loading  rates  was  correlated  to  pH.  Thus,  the  Kurbatov  equation  can  be  applied  to  our 
system  as  follow's; 

log(Fc  index)  * -xpH  - log  Kp  + log  N„  3-1 

where  x is  a macroscopic  proton  coefficient,  Kp  a Kurbatov  partition  coefficient  and  N, 
the  total  number  of  exchangeable  sites  (Fuller  et  al.,  1 996,  Balistrieri  et  al.,  1 983). 
Assume  that  the  above  equation  is  correct,  then  the  Fe  index  is  related  to  pH  in  addition 
to  x,  Kp  and  N,.  The  number  of  exchangeable  sites  N,  in  the  equation  for  a given  soil  is 


generally  a constant.  On  the  other  hand. varies  with  metal  species  in  solution  and  on 
surfaces  and  K,  does  not  vary  much  unless  the  adsorption  density  exceeds  its  maximum 
(Balistrieri  et  al..  1983).  Therefore,  for  a given  system,  the  Fe  index  changes  mainly  with 
pH  and  %.  The  relation  between  pH  and  the  Fe  index  in  the  Kurbatov  equation  is  as 
follows.  For  soils  with  high  Pb  loading,  pH  was  inversely  related  to  the  Fe  index  as 
expected  with  r2=0.55,  i.e.  the  index  decreased  as  pH  increased  from  6.1  to  6.8. 

However,  for  soils  with  low  Pb  loading,  pH  was  positively  related  to  the  Fe  index  with 
r2=0.6l,p<0.05  (Figure  3-IA&D).  The  impact  of  y on  the  index  is  unclear  at  the  present 

As  shown  in  the  Kurbatov  equation  (Eq.  3-1),  the  Fe  index  is  a lumped  parameter 
that  presents  the  properties  of  both  soil  solution  (pH  and  x)  and  solid  (Kp  and  N,)  that 
determine  Fe  partitioning.  To  evaluate  the  relationship  between  Fe  index  and  other  metal 
solubility,  it  is  expected  that  the  more  Fe-like  a metal  acts,  the  more  related  its  solubility 
is  to  the  Fe  index.  As  mentioned  previously,  Pb(ll)  shows  some  similarity  to  Fe  (11)  in 
both  metal  classiftcatron  and  speciation  distribution  in  aquatic  systems  (Turner  et  al., 
1981;  Nieboer  and  Richard,  1 980).  Therefore,  it  is  expected  that  aqueous  Pb 
concentration  in  soil  is  related  to  the  Fe  index  (Figure  3-1 D). 

Turner  et  al.  ( 1 98 1 ) reported  that  in  natural  waters  Pb  (II)  and  Fe  (II)  exhibit  some 
similarity,  meaning  that  Fe  (II)  may  be  a major  competitive  cation  to  Pb  (II)  for  ligands. 
Therefore,  it  is  important  to  evaluate  the  competition  between  Fe  (II)  and  Pb  (II)  for 
ligands.  It  should  be  pointed  out  that,  in  this  chapter,  we  did  not  emphasize  tire 
competition  between  Fe  (II)  and  Pb  (II)  for  adsorption  sites  on  solid  phases.  Tire  reason  is 
that  some  researches  have  shown  that  the  competition  may  be  insignificant  (Rose  and 
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Bianchi-Moaquera,  1993;  Coughlin  and  Slone,  1995).  In  our  study,  soil  solution 
chemistry  varied  in  different  soil  columns.  These  variations  may  shift  the  competition. 
At  the  present  time,  it  is  difficult  to  fully  evaluate  these  effects  on  metal  solubility  in  the 
soil,  especially  when  DOC  is  present.  The  Fe  index  presented  in  this  paper,  however, 
may  serve  as  a simple  measure  of  competing  ability  of  Fc  with  Pb  for  ligands.  At  this 
stage  it  is  more  important  to  look  at  the  conceptual  nature  of  the  index  than  the  detailed 
mechanistic  interpretation,  which  need  to  be  explored  further  in  future. 

Assuming  that  Fe  and  Pb  compete  for  ligands  in  a soil,  then  the  following  can  be 
inferred.  As  Fc  partition  index  increases,  there  is  more  aqueous  Fe  and/or  less  sorbed  Fc, 
which  means  there  is  more  competition  for  ligands  from  aqueous  Fe  in  the  soil  and  thus 
results  in  less  aqueous  Pb  in  the  soil.  On  the  other  hand,  as  Fe  partition  index  decreases, 


ligands  from  aqueous  Fe  and  thus  results  in  more  aqueous  Pb  in  the  soil  as  shown  in 
Figure  3-ID. 


In  theory,  there  should  be  a minimum  for  the  Fe  index,  i.e.  the  ratio  of  the 
minimum  of  aqueous  Fc  concentration  to  the  maximum  sorbed  Fc  concentrations.  If 
aqueous  Fe  concentrations  approach  zero,  i.e.  all  aqueous  Fe  was  transferred  onto  sorbed 
form,  the  Fe  index  is  zero.  At  this  situation,  Pb  dominates  the  ligand  competition  in  the 
solution,  so  that  enhanced  Pb  solubility  would  be  expected.  The  overall  trend  shown  in 
Figure  3-ID  confirmed  this  rationale  in  principle.  In  reality,  however,  aqueous  Fc 
concentrations  cannot  be  zero.  Therefore,  the  Fe  index  must  be  larger  than  zero. 
Theoretically,  in  a soil  the  minimum  of  aqueous  Fe  concentration  is  mainly  determined 
by  the  solubility  of  Fc  minerals  and  the  maximum  sorbed  Fe  concentration  is  soil- 


dependent.  Therefore,  the  Fe  index  is  generally  insensitive  to  solution  chemistry  in  a 
soil.  Since  it  is  impossible  to  accurately  measure  the  maximum  sorbed  Fe  concentrations 
in  a soil,  no  attempt  was  made  to  evaluate  the  real  value  of  the  minimum  Fe  index  in  the 
soil  we  used  from  the  definition.  However,  the  minimum  Fe  index  when  enhanced  Pb 
solubility  occurred  was  approximately  2 kg  L'1  in  the  soil  studied,  which  was  estimated 
from  the  relationship  between  Pb  solubility  and  Fe  index  (Figure  3-ID),  As  expected,  the 
number  did  not  vary  much  with  the  incubation  conditions  and  Pb  loadings  (Figure  3-ID). 
It  has  been  noticed  that  Pb  Concentrations  in  some  soil  columns  still  stayed  low  even 
when  the  Fe  index  was  < ~2  kg  L'1  (Figure  3-1 D),  which  is  unclear  at  the  present  time. 
Nevertheless,  the  significance  between  Pb  solubility  and  Fe  index  was  that  significantly 
enhanced  Pb  solubility  does  not  occur  unless  the  Fe  index  was  < -2  kg  L'1  (Figure  3-ID). 

Relation  of  Pb  solubility  and  Fe  partitioning  in  published  data 

In  this  research,  sorbed  Fe  was  defined  as  the  Fc  (U)  fraction  extracted  by  0.5  M 
HCI,  which  is  supposed  to  extract  exchangeable,  some  adsorbed  and  freshly-precipitated 
Fc(II)  (Heron  ct  al..  1994).  In  literature,  however,  such  data  are  not  always  available.  In 
order  to  conceptually  test  the  relation  between  Pb  solubility  and  Fe  partitioning  index  we 
defined  in  this  paper  using  available  data  in  the  literature,  other  extractants  that  may  be 
equivalent  to  0,5  M HCI  have  to  be  used.  Two  published  data  sets  were  used.  In  both 
data  set,  the  Fe  index  was  calculated  using  the  total  sorbed  Fc  instead  of  Fe  (II)  since  the 
quantity  of  the  later  is  unavailable.  However,  Fe  index  calculated  from  cither  total  Fe  or 
total  Fe  (II)  should  be  consistent  conceptually  since  Fc  (III)  is  even  more  competitive 
compared  with  Fe  (II). 


Fe  partition  index  (kg  L"') 
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Figure  3-3  Relationship  between  soluble  Pb  concentrations  and  the  ratio  of  soluble  to 
sorbed  Fe  in  contaminated  soils.  Data  arc  adapted  from  Karczewska  (1996). 


The  first  dale  set  was  from  Lee  et  al.  (1997)  who  characterized  a sediment 
contaminated  by  heavy  metals  in  a retention  pond  in  France.  In  the  study,  concentrations 
of  aqueous  Pb  and  Fe  in  pore  water  of  the  sediment  were  analyzed,  and  Fe  concentrations 
in  the  sediment  was  determined  by  sequentially  extracting  into  five  fractions: 
exchangeable  extracted  by  I M MgCI;.  bound  to  carbonate  by  I M sodium  acetate,  bound 
to  amorphous  Fe  and  Mn  hydroxides  by  0.04M  hydroxylamine  hydrochloric  acicl  in  25% 
acetic  acid,  bound  to  organic  matters  and  sulfides  by  30%  H2O2  and  0.02  M nitric  acid, 
and  residual  by  concentrated  HNOj  and  HClOj.  Assuming  that  the  sum  of  Fe 
concentrations  in  the  first  two  fractions  (exchangeable  and  bound  to  carbonate)  was 
equivalent  to  that  extracted  by  0.5  M HCI,  a plot  of  aqueous  Pb  concentrations  to  the  Fe 

observed:  enhanced  Pb  solubility  occurred  only  when  the  Fe  index  was  low.  However. 

This  again  suggests  that  high  Pb  solubility  will  not  occur  unless  the  Fe  index  was  below  a 

The  second  data  set  was  taken  from  Karczewska  (1996)  who  determined  concentrations 
of  heavy  metals  via  a sequential  extraction  in  soils  polluted  by  a cdfper  smelter.  In  the 
study,  heavy  metals  were  fractionated  into  7 fractions,  i.e.,  soluble  (1M  NH4NO3), 
exchangeable  (IMNH4OAC),  bound  in  MnO,  (IM  NH20H-HC1/1M  NHrOAc)  and  etc. 
We  assumed  that  Pb  extracted  by  1M  NHjNOj  is  equivalent  to  the  aqueous  Pb 
concentration  and  sum  of  exchangeable  Fe  and  bound  to  MnO,  was  equivalent  to  Fe(II) 
extracted  by  0.5  M HCI.  Again,  a similar  trend  shown  in  Figure  3-ID  was  observed 
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between  soluble  Pb  and  Fe  index  (Figure  3-3),  i.c.  highest  Pb  solubility  was  observed  at 
the  lowest  Fe  partition  index  and  in  this  particular  case  the  Fe  index  was  zero.  Like  the 
dataofLccctal.  (1997),  the  soil  samples  were  collected  from  various  locations  near  a 
smelter  site  in  Poland  and  there  were  large  variations  in  soil  properties  among  the 
samples  (CEC:  10-84  meq  kg'1;  organic  C:  0.05-1.42  %;  clay  content;  1-6  %) 
(Karczewska.  1996).  This  suggested  that  the  relationship  between  Pb  and  Fe  partition 
index  may  be  applicable  to  field  data  even  with  large  spatial  variability  in  soil  properties. 

I.nip).icati.gn.pfjhj.jjj:esearoh 

Partition  coefficients  (kj)  are  important  parameters  in  assessing  the  potential 
impacts  from  metal  contaminated  soils.  However,  kj  is  not  a constant  in  a dynamic  soil 
environment.  It  is  determined  not  only  by  the  characteristics  of  solid  phases  but  also 
solution  chemistry.  In  a given  soil  system,  variation  is  mainly  determined  by  metal 
solubility.  There  are  several  speciation  models  available  to  describe  the  interactions  of 
metals  and  soil  components  to  predict  their  solubility.  Unfortunately,  such  predictions 
generally  lack  certainty  primarily  because  of  soil  dynamic  nature. 

At  the  present  time,  assessing  changes  in  solution  chemistry  in  field  condition  is 
difficult  if  not  impossible.  This  information,  however,  is  critical  to  estimate  the  solubility 
of  heavy  metals  in  speciation  models.  In  contrast,  the  approach  presented  in  this  chapter 
provides  a simple  relation  between  the  probability  of  enhanced  Pb  solubility  (or  kj)  and 
Fe  partition  index,  which  is  generally  not  sensitive  to  changes  in  solution  chemistry.  In 
this  approach,  we  simply  determine  the  minimum  Fe  index  in  a field  condition  to  predict 
the  possibility  of  enhanced  Pb  mobility  (the  lower  ka).  Of  course,  more  work  is  needed  to 


further  test  this  concept  in  dilTercnt  soils  and  determine  the  ranges  of  Fe  index  for 
enhanced  Pb  mobility. 


Conclusion 

Our  data  using  a synthetically  contaminated  soil  demonstrated  that  enhanced  Pb 
solubility  was  only  observed  at  low  Fc  partition  index.  Such  a trend  was  also  observed 
using  two  published  data  sets  in  the  literature,  thus  further  validating  the  relationship 
conceptually.  Since  data  used  in  this  paper  were  from  various  soil  environments 
(roadside  sediment  in  France.  Pb  contaminated  soil  in  Poland,  a sandy  soil  from  Florida), 
it  is  reasonable  to  conclude  that  Pb  solubility  may  be  indeed  related  to  Fe  partition  index 
in  soils.  This  concept  may  also  be  applied  to  solubility  of  other  heavy  metals  such  os  Cd, 
Cu,  Ni  and  Zn. 


CHAPTER  4 

HEAVY  METAL  MOBILITY  IN  CONTAMINATED  SOILS:  PART  1.  ROLE  OF 
EXCHANGE  SITES  IN  CONTROLLING  SOLUBILITY  AND  MOBILITY  OF 
HEAVY  METALS 

Introduction 

Heavy  metal  mobility  in  soils  is  of  environmental  significance  due  to  its  toxicity 
to  both  humans  and  animals  (Ma  et  al„  1 995).  At  a constant  water  flow  condition  in  a 

that  the  solubility  of  heavy  metals  in  soil  is  mainly  regulated  by  adsorption,  precipitation 
and  ion  exchange  reactions.  Although  much  effort  has  been  spent  to  model  heavy  metal 
solubility  (Cederberget  ol..  19S5;  Sposito.  I9S4).  such  a prediction  under  field  conditions 
contains  large  uncertainty.  It  is  partially  because  of  the  difficulty  in  assessing  the  effects 

solution  chemistry,  such  as  pH.  redox  potential  and  ionic  strength,  may  shift  the  retention 
processes  of  heavy  metals  significantly.  These  impacts  may  be  liulher  complicated  by 
the  competition  of  cations  for  ligands,  which  may  enhance  heavy  metal  mobility  under 
certain  conditions  (Amrhein.  et  al.,  1994). 

Soil  redox  status  varies  temporally  and  spatially.  In  a surface  soil  it  is  influenced 
by  rainfall,  bioactivity,  and  changes  in  land  use,  whereas  in  vadose  zone  mostly  by 
fluctuation  in  water  table  (Buol  et  al..  1997).  A reduction  in  redox  potential  may  cause 
changes  in  metal  oxidation  state,  formation  of  new  low  soluble  precipitates,  and  Fe 


: Chuan  cl  al„  1996; 
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ihc  role  of  exchange  phase  in  comrolling  the  solubility  and  mobility  of  heavy  metals  in 
incubated  soils  has  not  been  reported.  However,  this  knowledge  may  improve  our 
understanding  of  the  controversy  discussed  above. 

In  Chapter  3.  Pb  solubility  was  examined  in  a sandy  soil  spiked  with  Pb  and 
incubated  for  40  d under  water-flooded  or  non-water-flooded  conditions.  Solution 
chemistry  in  soil  columns  was  adjusted  using  different  concentrations  of  NaCl  and  CaClj 
and  deionized  water  of  varying  pH  before  incubation.  The  results  show  that  Pb  solubility 
in  the  incubated  soil  can  be  related  to  the  ratio  of  aqueous  Fc  to  0.5  M HCI  extracted  Fe 
(II)  from  solid  phases  instead  of  soluble  Fe  alone.  If  we  assume  0.5  M HCI  extractable 
Fe(II)  from  the  solid  phases  is  proportional  to  exchangeable  Fc  (II),  it  is  suggested  that 
exchange  sites  have  significant  effect  on  heavy  metal  solubility  in  soils. 

two  naturally  contaminated  soils.  In  addition  to  Pb.  concentrations  of  Cu  and  As  will  be 
presented  for  comparison.  The  major  objective  of  this  paper  was  to  examine  the  role  of 
exchange  sites  on  the  solubility  and  mobility  of  Pb.  Cu,  and  As  in  soils  during  incubation. 

Materials  and  Methods 
Location  and  characteristics  of  soil  samples 

The  two  Pb-contaminated  soils  used  in  this  study  were  from  Montreal,  Canada 
and  Tampa,  Florida,  which  were  exposed  to  Pb-battery  recycling  and  concurrent  smelting 
operation  in  the  past.  Selected  characteristics  of  the  soils  are  1 


listed  in  Table  4-1, 
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Column  Experiment 


mesh)  was  packed  in  ihc  bottom  of  60-mL  columns  (13  cm  in  length  x 2.6  cm  in' 
diameter),  then  30  g air-dried  soil  (moisture  content  = 9.0%)  was  packed  on  top  of  the 
sand  layer.  Preliminary  data  showed  that  the  sand  layer  helps  to  prevent  the  outlet 
clogging  in  a column  during  incubation  and  leaching  ad  has  no  detectable  impact  on  soil 
colloid  mobility.  During  each  step,  the  columns  were  gently  shaken  horizontally  for  a 
few  minutes  to  minimize  the  packing  effects.  This  resulted  in  2.6x5.2  and  x7.8-cm  soil 
columns  for  dte  Montreal  and  Tampa  soils,  respectively.  The  soil  columns  were  set 
vertically  and  prewetted  by  pumping  deionized  distilled  water  (DDW)  into  the  bottom  of 
the  columns  until  the  water  level  in  the  columns  was  above  the  soil.  These  columns  were 
then  incubated  under  water-flooded  condition  for  approximately  3, 20  and  80  d for  the 
Tampa  soil  and  3. 20, 30, 40. 50.  and  60  d for  the  Montreal  soil. 

Leaching.  Once  the  redox  potential  monitored  had  reached  to  a certain  level,  the 
water  on  top  of  the  soil  column  was  removed  and  a rubber  stopper  was  put  on  top  to  seal 
the  syringes.  The  syringes  were  then  turned  upside  down.  The  lower  side  of  the  column 
was  connected  to  a needle  as  an  inlet  for  the  influent,  0.01  M CaCL.  The  pumping  rate 
was  1.20  mL  min'1.  The  effluent  from  the  lop  of  the  column  was  collected  with  a 
fractionation  collector. 

Analysis  of  leachates.  Each  fraction  of  effluents  was  filtered  through  0.22  pm 
membrane  filter,  then  acidified  to  pH<2  to  yield  the  soluble  metal  concentrations.  The 
details  can  be  found  in  Chapter  3. 


Analysis  of  pore  water  of  .sfli! 


water  was  separated  from  the  soil  by  the  centrifugation  method  described  by  Nkedi-Kizza 
et  al.  (1982)  and  Dao  and  Levy  (1978).  Then  the  pore  water  was  filtered  with  0.22-pm 
membrane  filter  for  metal  analysis.  For  detailed  procedure  of  Fo(U)  analysis,  please  see 
Chapter  3. 

Result  and  Discussion 

Concentrations  of  heavy  metals  and  Fe  (11)  in  pore  water  varied  with  incubation 
lime  (Figures  4-1  and  4-2).  Fe  behaved  differently  in  the  two  soils  with  incubation:  in  the 
Montreal  soil.  Fc(ll)  decreased  with  incubation  (Figure  4-1 ) whereas  in  the  Tampa  soil  it 
increased(Figure  4-2).  Similarly,  Pb  increased  in  the  Montreal  soil  and  decreased  in  the 
Tampa  soil  with  incubation.  However,  in  both  soils  the  concentrations  of  Fe  (II)  and  Pb 
were  inversely  correlated  with  incubation  (Figure  4-1  and  4-2).  Concentrations  of  As  and 
Cu  increased  first  then  decreased  with  incubation  in  the  Tampa  soil. 
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Figure  4- 1 . Pb  and  Fc  (II)  concentrations  in  pore  water  of  the  Montreal  soil. 
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Figure  4-2  Concentrations  of  Pb,  As,  Cu  and  Fc(II)  in  pore  water  of  the  Tampa  soil 
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dissolved  Fe(II)  increased  due  10  Fe(III)  reduction  to  Fe(ir)  (Ararhein  et  al. , 1994). 
However  the  unexpected  behavior  of  Fe  (II)  in  the  Montreal  soil  may  be  caused  by 
formation  of  FeCOj,  which  will  be  discussed  in  Chapter  5 in  details. 

Chuan  et  al.  (1996)  examined  the  release  of  heavy  metals  from  a contaminated 
soil,  which  was  suspended  in  water  (soihwater  ratio  = 1:7),  and  found  that  aqueous 
concentrations  of  Pb,  Cd  and  Zn  ore  positively  correlated  with  those  of  Fe  during 
incubation  because  heavy  metals  are  released  from  Fe  (oxy)hydroxide  surface  os  it 
dissolves  (Chuan  et  al.,  1996).  Obviously,  it  is  contrary  to  what  we  observed  for  Pb 
concentrations  in  both  the  Tampa  and  Montreal  soils.  On  the  other  hand,  Amrhein  et  al, 
( 1 994)  reported  that,  under  water-saturated  (soihwater  ratio  - 1 :0.2)  incubation,  Cd  and 
Cu  concentrations  in  pore  water  decrease  whereas  Fe(II)  increases  with  time.  This  is 
consistent  with  our  results.  It  is  possible  that  the  different  relations  between  Fe(U)  and 
heavy  metals  in  solution  result  from  the  different  soil-water  ratios  used  in  different 
studies.  In  the  studies  discussed  above,  contributions  of  exchangeable  metals  to  aqueous 
metal  concentrations  during  incubation  varied  with  the  soihwater  ratios.  When  the  ratio 
is  low,  he.,  small  quantity  of  soil  is  suspended  in  large  amount  of  water,  ions  on  the 
exchange  sues  tend  to  difluse  into  the  water  to  maintain  their  chemical  potentials 
between  solution  and  exchange  sites.  Consequently,  the  importance  of  exchange  sites  in 

increases,  the  importance  of  exchangeable  metals  in  contributing  to  soluble  metals 
increases  since  the  amount  of  exchange  sites  increases  relative  to  solution  volume. 
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consideration,  which  is  mainly  controlled  by  metal  competition  between  exchange  sites 
and  aqueous  ligands  in  solution.  Therefore,  the  observation  of  Amrhein  et  al.  (1 994)  and 
the  one  presented  here  may  be  due  to  the  fact  that,  with  incubation,  more  heavy  metals  go 
to  exchange  sites  as  Fc  (II)  concentrations  increase  in  solution.  We  have  demonstrated 
that  Pb  solubility  in  soil  columns  is  inversely  related  to  aqueous  Fe  concentration  in 
Chapter  3.  Therefore,  in  this  sense,  all  the  studies  discussed  above  along  with  the  one 
presented  here  are  consistent.  i.e.,  with  Fe  dissolution  upon  reduction,  Pb  is  released  into 

Similarly,  arsenic  concentrations  in  pore  water  and  suspension  varied  differently 
with  incubation.  Masschclcyn  ct  al.,  (1991)  examined  Ihc-orscnic  release  from  a 

soluble  arsenic  increases  as  redox  potential  decreases  (+500  -200  mv).  It  was  attributed 
to  As (V)  reduction  to  As(Ill)  and  Fe  reduction  dissolution  to  release  adsorbed  As. 

different  trend  (Onken  and  Hossner,  1996).  They  found  that  soluble  arsenic 

occurred  at  20-30  d after  flooding.  Further  analysis  of  the  solid  phase  showed  the  loss  in 
soluble  atsenic  could  be  accounted  by  surface  bound  when  the  soil  was  incubated  longer 
than  20-30  d.  The  result  of  the  later  study  is  consistent  with  our  data  that  showed  As 
concentration  increased  and  then  decreased  with  Fc  reduction  dissolution  (Figure  4-2), 


incubation  time  (d) 


Figure  4-3.  Cumulative  Pb  and  Fe  leached  after  3 1 .8  pore  volumes  of  0.01  M CaCI2 
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Figure 4-4.  Cumulative Pb,  As,  Cuand  Fe  leached  alter  31.8  pore  volumesofO.Ol  CaCl; 
in  Tompa  soil. 


Heavy  metal  mobililv  with  incubation 


Mela!  mobility  in  this  paper  is  examined  using  cumulative  metais  leached  with  32 
pore  volumes  of  0.0 1 M CaClj  solution.  It  is  well  understood  that  the  metals  that  can  be 
leached  out  with  CaClj  solution  include  both  the  soluble  and  the  exchangeable.  For  the 
Montreal  soil,  the  cumulative  leached  Pb  concentrations  increased  with  incubation 
(Figure  4-3).  which  is  consistent  with  the  overall  trend  shown  in  the  pore  watets  (Figure 
4-1 ).  However,  the  cumulative  leached  Fe  concentrations  showed  a different  trend  (torn 
Fe(II)  concentrations  in  pore  waters  (Figure  4-3),  suggesting  that  Fe(II)  in  exchange  site 
is  not  proportional  to  that  in  aqueous  phase  assuming  that  Fe(Il)  * soluble  Fe.  For  the 
Tampa  soil,  changes  in  overall  trends  for  metals  with  incubation  can  be  divided  into  two 
groups:  I .)  Pb  and  Cu  decreased  with  incubation;  2.)  As  and  Fe  increased  with 
incubalion.  which  is  not  consistent  with  the  case  in  pore  water. 

In  fact,  the  reduction  of  cation  exchange  capacity  (CEC)  caused  by  the  association 
of  hydroiron  polymers  or  iron  hydroxides  with  soil  mineral  has  been  long  recognized 
(Coleman  and  Thomas.  1964;  Corsica  et  al.,  1970).  Hcndershot  and  Lavkulich  (1983) 
reported  that  Fe  coatings  on  illite  lower  CEC  and  the  ones  on  kaolinite  increase  anion 
exchange  capacities  (AEC)  in  comparison  to  the  uncoated  ones  at  pH<  7,  They 
suggested  that  the  Fe  coatings  reduce  the  CEC  and  increase  AEC  measured  at  low  pH  by 
cither  physically  blocking  or  electrostatically  canceling  the  permanent  negative  charge 
carried  by  the  crystalline  minerals.  Similarly,  Arias  etal.,  ( 1995)  found  that  the  point  of 
zero  charge  (PZC)  of  Fe-coaled  kaolinite  increases  linearly  with  the  amount  of  the 
coating  Fe.  In  the  experiments  by  Shainberg.  et  al.  (1987).  they  demonstrated  that  the 
CEC  of  the  mixture  Fe  (FeClj)  and  soil  decreases  by  52%  in  the  1 0 mmol  Fe  kg'1 
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treatment  to  72%  in  the  40  mmol  Fe  kg'  treatment,  and  correspondingly  the 
exchangeable  sodium  percentage  from  23.8  to  15.4.  Rcngasamy  and  Oades  (1977) 
reported  that  with  addition  of  ferric  hydroxide  the  electrophoretic  mobility  of  iron-clay 
increase  from  negative  to  positive  values.  All  these  data  surpport  that  with  the  reduction 
dissolution  of  Fe,  Fe  coating  fraction  decreases,  and  then  CEC  increases  for  a soil.  Even 
though  it  is  supported  by  sufficient  experimental  and  theoretical  evidence,  it  has.  to  our 

with  solution  chemistry,  can  have  significant  effects  on  heavy  metal  retention. 

concentration  in  exchangeable  phases  may  increase  as  well.  It  was  evident  that,  for  the 

amount  with  incubation  time  were  similar:  they  both  increased  with  incubation  (Figure  4- 

found  at  20  day  incubation  compared  to  Fe(ll)  concentrations  in  pore  water,  suggesting 
that  exchangeable  Fe  (II)  became  dominant  in  cumulative  leached  Fe.  In  fact,  there  was 
5.55  pmoles  Fe  leached  during  7-16  pore  volumes  (which  can  be  considered  as 
exchangeable)  among  7. 1 1 pmoles  of  the  total  cumulative  leached  Fe  after  20-d 
incubation. 

In  the  Tampa  soil,  changes  in  Pb  concentrations  in  pore  water  and  the  cumulative 
leached  Pb  with  incubation  time  were  similar:  they  both  decreased  with  incubation.  This 

controlled  by  competition  between  aqueous  Fe  and  Pb  for  ligands  and  exchange  sites.  It 


nted  that  Fe(II)  and  Pb(II)  have  similar  affinity  for  ligands  (To 


1981 ; Niebocr  and  Richard,  1980).  As  aqueous  Fe  (II)  concentrations  increase  with 
exchange  sites.  If  this  is  simply  the  case,  however,  the  cumulative  leached  Pb  was  not 
that  there  was  so  much  Pb  present  in  the  exchange  sites  without  some  of  it  adsorbing 


chemically  onto  the  surfaces 

. However,  it  should  be  pointed  out  that  there  was 

significant  difference  betwet 

:n  the  relations  of  Pb  concentrations  in  pore  water  and 

cumulative  leached  Pb  with 

incubation:  while  the  Pb  concentration  in  pore  water 

decreased  exponentially  the. 

cumulative  leached  Pb  decreased  only  linearly  (incubation 

time  <30  day).  It  is  suggest 

ed  that  less  reduction  of  the  cumulative  leached  Pb  than  that 

of  Pb  concentrations  in  pore 

In  the  Tampa  soil,  cu 

mulativc  leached  Cu  decreased  with  incubation  (Figure  4-4), 

significantly  greater  cumulat 

ivc  leached  Cu  afier  3-d  incubation  can  be  attributed  to  more 

Cu  coming  from  exchange  si 

tes  when  leached  with  CaClj.  In  fact,  Cu  (II)  has  greater 

affinity  for  exchange  sites  Ih 

an  Pb  (II)  (McBride,  1994). 

Interestingly,  As  con, 

centrations  in  pore  water  increased  first  and  then  decreased 

with  incubation,  whereas  the 

amount  of  cumulative  leached  As  showed  no  significant 

drop  with  incubation  (Figure  4-4),  As  discussed  above,  as  Fe  “coatings"  was  removed 


with  incubation  via  reductior 

i dissolution,  soil  CEC  increased  as  well.  This  caused  more 

Pb  (II)  to  transfer  from  soluti 

ion  to  exchange  sites,  as  a result,  local  charge  reversal  might 

occur,  therefore  positively-charged  particles  held  more  HsAsOT  and  /or  H3ASO4*.  It  has 
been  realized  that  at  slightly  reducing  condition  and  pH  = 4 (for  the  Tampa  soil  in  this 
study)  both  H2ASO3'  and  H3ASO4’  arc  not  stable  thermodynamically  (Cullen  and  Reimer, 
1989).  However,  local  pH  near  positively  charged  surface  may  be  significantly  greater 
than  in  bulk  solution,  therefore  H-AsOj'  may  be  stable  near  highly  positively-charged 
surfaces.  On  the  other  hand,  there  has  been  evidence  showing  that  even  under  reducing 
condition,  only  HjAsOT  is  associated  with  natural  particle  surfaces  (Belzile  and  Tessier , 

1 989).  The  existence  of  HjAsO/  under  reducing  condition  has  been  attributed  to  the 
oxidation  of  arsenite  to  arsenate  by  Mn  and  Fe  oxyhydroxides  (Peterson  and  Carpenter. 
1983). 

As  discussed  above.  CEC  increases  with  Fc  reduction  dissolution,  especially  at 
pH  < PZC  of  Fe  minerals.  The  consequent  redistribution  of  metals  between  solution  and 
exchange  sites  is  probably  determined  not  only  by  the  magnitude  of  CEC,  but  also  the 
characteristics  of  metal  ions,  the  competing  co-ions  and  counter-ions  (ligands). 

Therefore,  changes  in  metal  concenlradons  in  a soil  solution  with  time  under  incubation 
can  be  complicated.  However,  the  following  interactions  have  to  be  taken  into 
consideration: 

1 . Fe  reduction  dissolution  releases  Fc(II)  and  adsorbed  heavy  metals  into  solution  at  the 
same  time  frees  exchange  sites, 

2.  Fe  (II)  and  heavy  metal  cations  compete  for  ligands  to  stay  in  solution, 

3.  Fe  (II)  and  heavy  metal  cations  compete  for  the  exchange  sites,  transferring  from 
solution  to  solid  phase. 


4.  Charge  reversal  on  exchange  sites  results  in  negatively  changed  heavy  metal  to 
transfer  from  solution  to  exchange  sites. 

In  a particular  soil  system,  different  processes  may  be  dominant  in  controlling  metal 
concentrations  in  soil  solution. 

There  is  an  alternative  explanation  for  the  As  behavior  with  incubation.  It  has 
been  well  established  that  P04J‘  may  precipitate  with  Fe!*  to  form  minerals  such  os 
vivianite.  Fej(P04)j  8H  jO  at  low  soil  redox  potential  (Holford  and  Patrick,  1 979;  Patrick 
and  Khalid,  1974).  Because  of  the  great  similarity  between  HjPO«'  and  HiAsOj'or 
HjAsOj",  it  can  be  inferred  that  formation  of  a new  solid  phase  by  Fe!’  and  HjAsCV  or 
HjAsOj  limited  the  solubility  of  As  in  pore  water  as  the  incubation  time  approached  to  60 
days  (Figure  4-4). 


implication  of  this  research 

In  literature,  there  has  been  a fair  amount  of  papers  focusing  on  changes  in 
solubility  and  mobility  of  heavy  metals  with  incubation.  However,  reported  data  have 
often  been  contradictory  partially  because  of  differences  in  soil/water  ratios  used.in 


separated  from  soil  suspension,  and  leachates  from  a soil  leached  with  DDW  or  an 
electrolyte.  It  has  been  noticed  that  concentration  changes  of  heavy  metals  with 
incubation  are  often  confusing.  For  example,  metal  concentration  in  pore  water  or  metal 
mobility  leached  by  DDW  suddenly  decreases  with  incubation  when  occurrence  of 


precipitation  can  not  be  proved.  As  discussed  in  this  chapter,  the  roles  of  exchange  sites 
in  controlling  metal  concentrations  in  these  systems  have  to  be  taken  into  consideration. 
With  sufficient  evidence,  we  showed  that  CEC  of  a soil  change  with  incubation,  and  thus 
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may  greatly  affect  metal  solubility  and  mobility.  The  role  of  exchange  sites  on  the 
solubility  and  mobility  of  heavy  metals  is  crucial  to  understand  various  behaviors  of 
heavy  metal  upon  incubation.  Furthermore,  incorporation  of  these  results  into  fitlure 
metal  transport  models  may  improve  their  accuracy. 


CHAPTER  5 

HEAVY  METAL  MOBILITY  IN  CONTAMINATED  SOILS:  PART  2.  COLLOID- 
FACILITATED  METAL  MOBILITY  IN  A PB-CONTAMINATED  SOIL 

Introduction 

Understanding  metal  mobility  in  soils  is  of  significant  environmental  interest.  In 
general,  it  is  widely  accepted  that  metal  mobility,  such  as  Pb.  is  practically  negligible  in 
soils  due  to  their  low  solubility.  However,  enhanced  Pb  mobility  occurs  in  soils  under 
certain  conditions  (Newman  ct  al.,  1993).  Colloid-facilitated  metal  transport  in  soils  has 
been  proposed  as  one  of  the  mechanisms  (McCarthy  and  Zachara,  1 9S9). 

Based  on  the  DLVO  (Derjaguin-Landau-Verwey-Overbeek)  theory  (Israelachvili,  1992), 

electrical  double  layer  (EDL)  and  van  dcr  Walls  forces,  which  varies  with  separation 

secondary  minimum  first,  then  a primary  minimum  alter  an  energy  barrier.  After  that 
point,  firnher  panicle  move  towards  each  other  causes  drastic  increases  in  the  interaction 
energy,  which  is  referred  to  as  Bom  repulsion.  In  principle,  panicles  are  attached 

minimum.  On  the  other  hand,  both  the  magnitude  of  the  energy  barrier  and  the 
separation  distances  for  the  primary  and  secondary  minima  are  affected  by  solution 
chemistry.  In  order  for  a deposited  panicle  to  be  released,  repulsive  forces  have  to  be 
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generated  between  surfaces  of  colloid  and  stationary  solid  phases  such  as  resulting  from 
changes  in  solution  chemistry.  Changes  in  solution  chemistry  may  reduce  or  eliminate 
the  primary  or  secondary  minimum  so  that  the  release  of  deposited  particles  become 
favored  (McDowcll-Boyer,  1992;  Ryan  and  Gschwend,  1994b).  Numerous  studies  have 
shown  that  solution  chemistry  is  one  of  the  primary  factors  affecting  particle  deposition 
and  release  in  natural  porous  media  (Elemeleeh  and  O'Mclia,  1990a;  Elimelech,  1992; 
McDowell-Boycr,  1992;  Ryan  and  Gschwend.  1994;  Seaman  et  al.,  1995).  However 
influence  of  solution  chemistry  on  the  interaction  between  EDLs  of  particles  is  a rather 
complicated  issue.  In  addition,  solution  chemistry  in  soil  varies  spatially  and  temporally, 
therefore,  it  is  still  a challenge  to  model  the  interaction  of  EDLs  between  colloids  and 
stationary  solid  phases  in  soils. 

Nevertheless,  effects  of  solution  chemistry  on  colloid  mobility  are  closely  related 
to  potential  determining  ions  (Parks,  1967).  specific  adsorbed  ions  (Parks.  1975)  and 
indifferent  ions  (Park.  1967)  present  in  soil  solution.  Potential  determining  ions  are 
constituent  ions  of  solid  particles  (c.g.  H*  and  OH*  of  Fe  oxides,  Ca2*  and  CO)2*  of 
CaCOr)  and  their  concentrations  primarily  determine  particle  surface  potentials.  Specific 
adsorbed  ions  arc  those  adsorbed  onto  solid  surfaces  and  may  change  the  magnitude  and 
even  sign  of  the  surface  charge  (e.g.  Ca2*  and  Pb2‘),  whereas  indifferent  ions  are  those 
adsorbed  physically  to  surfaces  and  change  only  the  magnitude  of  the  surface  charge  (c.g. 
Na*  and  Cl").  It  is  important  to  consider  all  three  types  of  ions  when  studying  the  impacts 
of  dynamic  solution  chemistry  on  colloid  mobility  in  soils.  In  surface  soils,  solution  pH. 
ionic  strength,  and  redox  potential  are  influenced  by  rainfall,  bioactivity,  or  change  of 
land  use;  whereas  in  vadose  zone  they  are  mostly  influenced  by  the  fluctuation  of  water 
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(able  (Boul  el  al,  1 997).  In  other  words,  in  soils  all  three  kinds  of  ions  (potential 
determining  ions,  specifically  adsorbed  ions  and  indifferent  ions)  vary  drastically  with 
time  and  space  in  terms  of  species  and  concentrations.  Inevitably,  those  changes  will 
significantly  affect  the  interactions  between  colloid  surfaces  and  stationeiy  solid  phases, 
and  thus  colloid  release  and  deposition  rate. 

In  particular,  solution  chemistry  has  great  effect  on  surface  charge  heterogeneity 
of  soil.  Oxides  of  Fe,  Al  and  Mn  are  the  most  common  sources  of  surface  charge 
heterogeneity  in  natural  aqueous  environments  (Ryan  and  Elimelech,  1996).  At  neutral 
pH,  these  compounds  may  carry  a positive  surface  charge,  which  are  called  positively- 
charged  patches  in  matrix  minerals  such  as  quartz  and  feldspar  that  carries  negative 
surface  charge  (Song  et  al.,  1 994).  In  natural  aquatic  systems,  the  range  and  magnitude 
of  the  existence  of  these  patches  may  affect  colloid  deposition  and  release  significantly. 
Recently,  Johnson  et  al  ( 1 996)  have  demonstrated  that  colloid  deposition  rate  is 
controlled  by  the  extent  of  positively-charged  patches  developed  in  stationary  solid 
phases  in  porous  media.  In  dynamic  soil  environments,  however,  these  patches  may  vary 
temporally  and  spatially  since  solution  chemistry  (e.g.  pH  and  redox  potential)  greatly 
affects  oxide  solubility.  However,  the  effects  of  the  dynamic  feature  of  surface  charge 
hcierogencuy  in  soil  on  colloid  mobility  are  not  considered  in  most  studies, 

Metals  of  low  solubility  tend  to  associate  with  colloids  (Mills  et  al.,  1991).  They 
may  adsorb  or  precipitate  onto  colloid  surfaces,  or  even  form  their  own  colloids.  Take 
Pb2'  for  example,  in  the  first  case,  Pb2"  may  affect  the  surface  potential  between  solution 
and  solid  drastically  by  acting  as  a specific  adsorbed  ion  (Kretzschmar  et  al.,  1997).  In 
the  second  case,  Pb‘*  is  a potential  determining  ion  that  may  change  the  surface  potential 
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greatly.  In  other  words,  the  mobility  of  these  colloids  may  be  different  from  the  original 
ones  during  charge  developments.  The  mobility  of  metal-bearing  colloids  thus  needs  to 
be  examined  (Temmingholfetal.,  1997;  Krctzschmar  ct  al„  1997). 

The  objective  of  this  study  was  to  evaluate  the  effects  of  soil  solution  chemistry 
on  colloid-facilitated  metal  mobility,  such  as  Fe,  Al,  and  Pb,  in  a Pb-contaminatcd.  Such 
effects  were  limited  to  those  rising  from  surface  properties  of  colloids  and  stationery 
solid  phases,  which  are  impacted  by  soil  solution  chemistry.  In  our  experiment,  soil 
solution  chemistry  was  altered  through  incubation  under  saturation  condition.  Alter 
prewashing  the  soil  columns  with  0.01  M CaClj  solution,  the  influent  was  switched  to 
deionized  water  (DOW)  to  mobilize  colloids.  Short  columns  were  chosen  to  avoid 
colloid  redeposition.  Our  results  showed  that  colloid  mobility  varied  significantly  with 
incubation,  which  may  be  caused  by  changes  in  surface  charge  heterogeneity  on  colloids 
and  stationery  solid  phases. 


Materials  and  Methods 

Characteristics  of  soil  samples 

The  lead-contaminated  soil  used  in  this  study  was  collected  from  Montreal, 
Canada.  The  site  was  contaminated  via  Pb-battery  recycling  and  concurrent  smelting 
operation  in  the  post.  The  soil  was  air-dried  and  passed  through  2-mm  sieve.  Total 
elemental  concentrations  of  the  soil  was  determined  with  a Thermo-Jarrell-Ash  61 E ICP 
after  digestion  in  a CEM  MDS-2000  microwave  using  EPA  method  3051.  The  soil  was 
also  analyzed  using  x-ray  diffraction  and  thermal  analytical  techniques,  i.e. 


Table  5-1  - Selected  properties  of  the  Pb  contaminated  soil  used  in  this  study 


. Mi':::: . 

pH  ( 1 : 1 with  distilled  water)  7.7 

Fraction  of  organic  carbon  2.24 

Total  Elemental  Analysis  (Ma.  1996)  g kg'1 

Pb  1.60 

Fe  133.30 


Particle  size  distribution  (Day.  1965) 
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thermogravimetry  and  differential  thermal  analysis.  The  major  minerals  present  in  this 
soil  include  Fe  oxides  (amorphous),  quartz,  calcite,  dolomite,  vermiculite.  mica,  and 
feldspar.  Selected  soil  properties  arc  shown  in  Table  5- 1 . 


This  experimentation  is  a continuation  of  Chapter  4 

Leaching.  After  displacing  the  original  soil  solution  as  described  In  Chapter  4. 
the  influent  was  switched  to  DDW.  The  effluent  from  the  top  of  the  column  was 


Analysis  of  effluent.  Each  fraction  of  effluents  was  split  inti 
acidified  with  HNOj  to  pH  < 2 without  filtration.  The  metal  conceit 


i parts.  One  was 

>ns  determined 

soluble  Pb).  The 
ine  filter  to 


Redox  potential  normally  decreases  with  incubation  under  water-flooded 
condition.  In  our  study,  the  electrode  installed  in  soil  columns  measured  soil  redox 
potential  only  in  a relative  sense  as  discussed  previously.  In  addition.  Stumm  ( 1 984) 
suggested  that  it  is  more  reasonable  to  measure  important  redox  species  than  redox 
potential  to  indicate  redox  status  in  natural  systems.  Fc(II)  concentration  in  pore  wate 
was  thus  used  as  a measure  of  soil  redox  status. 


Pore  water  was  separated  from  the  soil  by  centrifugation  (Nkeda-Kizza  ct  al, 

1 982;  Dao  and  Lavy,  1 978)  after  removing  standing  water  on  top  of  the  soil  columns. 
The  pore  water  was  then  filtered  with  0. 1 -pm  membrane  filter  for  metal  analysis. 

Aliquot  ofO.I  mLofthe  filtrate  was  transferred  to  I0mLoffcrrozine(lgL',)in50mM 
HEPES  (N-2-hydroxyethylpiperazine-N-2-ethanesulfonic  acid)  buffer  at  pH=7  (Lovely 
and  Philips.  l987;Hcionetal.,  1994).  Concentration  ofFe  (II)  was  determined  by 
measuring  the  absorbance  ofthe  filtrates  at  562  nm  (Stookey,  1970).  In  addition,  Ca 


Analytical  methods 

labware.  All  chemicals  used  in  this  study  were  of  analytical  grade  or  better.  Double 
deionized  water  from  a Bamstead  NANOpure  water  system  was  used.  Total  metal 

spectrophotometer  (Perkin-EImer  2380)  equipped  with  a graphite  furnace  atomizer. 
Flame  atomic  absorption  was  used  to  analyze  metal  concentrations  > I mg  L’1  and 
graphite  furnace  atomizer  was  used  to  measure  metal  concentrations  < 1 mg  L'1.  Multi- 
level standards  (Fisher  Scientific)  for  all  metals  were  prepared  in  the  same  matrix  as  the 
extracting  reagents  to  minimize  matrix  effects.  Blanks  were  used  for  background 

run  with  every  20  samples  to  verify  method  precision.  The  spike  recovery  and  precision 
were  found  to  be  within  100%  + 10%. 


Result  and  Discussion 


£9M^jnaM!te.m«l.mg"b«ii°B 

Effluent  turbidity,  a measure  uf  relative  colloid  concentrations,  generally 
decreased  with  incubation  time,  with  longer  incubation  resulting  in  lower  turbidity 
(Figure  5-1 ).  Maximum  turbidity  was  reached  at  3-4  pore  volumes  for  all  soil  columns, 
suggesting  that  colloids  were  mobilized  by  switching  the  influent  from  0.01  M CaCl; 
solution  to  DDW.  Effluent  turbidity  stayed  more  or  less  constant  as  pore  volume 
increased  from  3-4  to  24.  Chloride  concentration,  which  was  20  mM  in  the  influent 
before  switch  to  DDW.  decreased  significantly  from  20  to  0.2  mM  as  pore  volume 
increased  from  0 to  1.6.  Meanwhile,  the  turbidity  increased  drastically  from  0.2  to  0.8 
NTU.  It  is  suggested  that  colloids  were  mobilized  by  DDW  immediately  while  0.01  M 
CaClr  in  the  soil  columns  was  displacing  by  DDW.  This  fast  colloid  release  has  been 
reported  by  Nocito-Gobel  et  al.  ( 1 996)  who  observed  that  polystyrene  latex  was  released 
from  a silica-sand  column  due  to  reduction  of  ionic  strength  in  the  influent.  They 
explained  that  reduction  in  ionic  strength  may  eliminate  prior  net  attractive  force  minima 
(primary  or  secondary)  so  that  colloid  release  rate  is  fast  relative  to  the  advective  solution 
flow  through  the  column. 

As  indicated  by  turbidity,  colloid  mobility  decreased  with  incubation  or  with 
decrease  in  aqueous  Fe(ll)  concentrations  from  0.91  to  0.1 7 pM  in  pore  water  (Figure  5- 
2A).  The  cumulative  leached  colloidal  Fe  or  Al  in  23  pore  volumes  after  3-d  incubation 
(aqueous  Fe  (II)  = 0.91  pM ) were  significantly  greater  than  they  were  after  20-  or  80-d 
incubation  (Fe  (1I)=0,27  or  0.17  pM ) (Figure  5-2  A).  It  has  long  been  recognized  that 


Figure  5-1  Changes  of  effluent  turbidity  with  pore  volumes  under  different  incubation 
times.  The  insert  is  a typical  breakthrough  curve  of  CT  when  deionized  water  displaced 
CaClj. 


redox  status  may  affect  colloid  mobility  in  soils  (McCarthy  and  Zachara,  1989; 
Grolimund  et  al.,  1996).  However,  there  is  little  information  available  on  this  topic  in  the 
literature.  In  Fe  oxide-coated  sand,  Ryan  and  Gschwend  (1994a)  reduced  goethite  with 
ascorbic  acid  to  remove  its  cement  effect,  and  found  that  the  amount  of  colloids  in 


effluents  is  directly  related  to  the  concentrations  of  dissolved  Fe.  This  is  seemingly 
consistent  with  our  result  in  that  cumulative  colloidal  Al  and  Fe  were  proportional  to 
aqueous  Fe  (11)  concentrations  in  pore  water  (Figure  5-2  A,  r2  - 0.92-0.97). 

In  this  study,  soil  redox  status  decreased  with  incubation,  which  was  confirmed  by 
decrease  in  redox  potentials  monitored  (Figure  5-2).  In  general,  a more  reducing 
condition  produces  a greater  aqueous  Fe  concentration  (Chuan  et  al.,  1996).  In  our 
system,  however,  the  opposite  was  observed,  i.c.  aqueous  Fe  (II)  concentration  in  pore 
water  decreased  from  0.9 1 to  0.17  pM  as  incubation  increased  from  3 to  80  d,  suggesting 
that  decementing  by  reductive  dissolution  of  Fe-oxide  may  not  be  the  only  factor 
controlling  colloid  mobility  in  our  system. 

Though  Fe  (11)  concentration  in  pore  water  decreased  with  incubation,  it  did  not 
change  much  when  the  incubation  time  was  increased  from  20  to  80  d (Figure  3-2), 
suggesting  that  Fe  concentrations  may  be  controlled  by  mineral  solubility  in  the  period. 
From  a thermodynamic  point  of  view,  sideritc  (FeCOs)  can  control  Fe1*  solubility  in  soil 
depending  on  redox  conditions  and  COj1"  concentrations  (Lindsay,  1979).  In  this  study, 
soil  redox  status  became  more  reducing  with  incubation,  which  was  confirmed  by  redox 


electrode  measurement  (Figure  5-2).  Therefore,  aqueous  Fe  (II)  was  released  by  Fe 
reduction  dissolution.  The  soil  used  in  this  study  contained  -8%  of  calcite  (CaCOy,  log 
K “ 9.74),  which  is  much  more  soluble  than  siderite  (log  K = 7.92,  Lindsay,  1979).  In 


Fe(II)  concentration  in  [tore  water  (uMI 


Ca  concentration  in  pore  water  (mM) 


Figure  5-2  Relation  of  aqueous  Fe  (II)  and  Ca  concentrations  in  the  pore  water  of  soil 
columns  on  the  cumulative  colloidal  Fe  and  Al  concentrations  in  the  effluent  after  23 
pore  volumes  during  incubation. 


; favored  as  a stable  ir 


ron  mineral,  which  was  confirmed 
by  M1NTEQA2  calculation  (Allison  et  a!.,  1991).  This  may  explain  why  aqueous  Fe(ll) 

Most  minerals  present  in  the  soil  used  in  this  study  carried  mostly  negative 
charges  as  expected  (Table  5-2).  For  minerals  carrying  pH  dependent-charge,  their  point 
of  zero  charge  (PZC)  was  generally  lower  than  pH  =7.7  of  the  soil  (Table  5-2).  The 
mobile  colloids  in  the  effluents  has  ^-potential  = -13  mv  in  0.1  M NaCl.  These  evidences 
supported  that  the  soil  was  negatively-charged,  however,  the  positively-charged  patches 
may  also  exist.  It  is  commonly  believed  that  Fc  oxides  are  major  sources  of  positive 
charge  near  neutral  pH.  This  led  us  to  expect  that  as  more  Fe  oxides  dissolved  with 
incubation,  positively-charged  patches  decreased,  colloid  mobility  would  increase. 
However,  the  opposite  relation  between  incubation  and  Fefll)  concentrations  were 
observed  in  this  study  (Figure  5-2A). 

As  discussed  previously,  Fe  oxides  may  be  transformed  to  siderite.  This 
transformation  promoted  aqueous  Ca  concentrations  in  pore  water,  which  was  consistent 
with  the  fad  that  Ca  concentrations  increased  with  incubation  (Figure  S-2B).  On  the 
surfaces  of  carbonate  minerals,  there  are  sufficient  examples  (Charlet  et  al.,  1990; 
Cappetlen  et  al.,  1993.)  showing  that  the  positive-charge  development  can  be  mainly 
described  as 

^OjH"  + Me2'  <-»  >COjMe'  + H*  (5-1 ) 

where  >COsH°  and  >COjMe*  stand  for  bicarbonate  sites  and  its  complex  with  tree  metal 
cation  (Me)  at  the  surface  of  carbonates,  respectively.  In  our  system,  Ca2*  was  dominant 
among  all  the  Me2*  in  pore  water,  for  example,  Ca  concentrations  was  much  greater 
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Table  5-2  Minerals  in  the  soil  and  iheir  Point  of  zero  charge  (PZC) 


Minerals  PZC 

Vormiculite 

Mica 

Fe  oxides  (amorphous)  7.9-8.2 
Quartz  2-3 

Calcite  5.5-6 

Dolomite 8-9 


Dzombaketal.,  1990 
Parks.  1967 
Parks  1967 

Parks  1967 
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than  Fc(II)  (Figure  5-2).  In  addition,  no  changes  in  pH  among  soil  columns  were  found 
(pH=~7.7),  which  may  be  due  to  the  buffering  capacity  of  carbonates.  Nevertheless,  This 
gave  us  a basis  to  assume  Ca2‘  concentrations  were  proportional  to  that  of  aqueous  Ca  in 
the  pore  water.  Obviously,  Eq  (5-1)  suggested  that  the  positive-charged  patches  on 
calcite  increased  with  aqueous  Ca  concentrations  and  thus  that  colloid  mobility  decreased 
because  more  mobile  colloids  had  deposited  on  the  stationary  grains  with  a greater 
bonding  (primary  minimum  deposition).  This  is  consistent  with  the  observation  that 
colloidal  A1  and  Fc  were  invctscly  related  to  aqueous  Ca  (Figure  5-2,  R!=0.43  and  0.86, 
respectively). 

Colloid  elution  curves 

Concentrations  of  colloidal  Fe  and  A1  changed  with  pore  volumes  arc 
demonstrated  in  Figure  5-3.  Different  from  the  observations  made  in  well  defined 
systems  (Nocito-Gobel  et  al„  1996)  that  the  elution  curves  consists  of  one  major  narrow 
peaks,  the  colloidal  metal  concentrations  peaked  over  a wide  range  of  pore  volumes  or 
repeatedly  at  several  numbers  of  pore  volume.  It  has  confirmed  that  the  patterns  .of  the 
elution  curves  are  reproducable.  The  velocity  difference  between  fluid  (tracer)  and 
colloidal  particles  have  been  long  recognized  (Small,  1974).  Size  exclusion  of  colloidal 
particles  could  cause  them  to  move  faster  than  the  tracer.  On  the  other  hand,  colloid 
deposition  in  a secondary  minimum  may  retard  colloids  significantly  (Nocito-Gobel  et 
al.,  1996).  To  our  knowledge,  however,  no  observation  has  been  reported  that  soil 


Figure  5-3  Elution  curves  of  colloidal  Fe,  Al,  and  Pb  concentrations  in  effluents  with  pore 
volumes  under  various  incubation  times.  Each  point  represents  the  mean  of  two  replicates 
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Figure  5-4  Relationship  among  colloidal  Fe  and  A1  concentrations,  dissolved  Ca 
concentrations  and  pH  in  effluents  after  3 d of  incubation.  Each  point  represents  the  mean 
of  two  replicates.  The  trend  is  reproducable.  The  arrows  indicate  the  boundary  Beyond 
which  the  Ca  release  mechanism  from  the  soil  may  change. 


In  facl,  after  the  displa 


»f  the  origin 


ail  solution  ir 


0.01  M CaCI;  solution.  Ca  was  the  dominant  cation  in  both  bulk  solution  and  electrical 
double  layers  (EDL)  surrounding  colloidal  particles  and  stationary  solid  phases  in  the 
soil.  When  the  effluent  was  switched  from  0.01  M CaCI2  solution  to  DDW,  we  expect 

play  probably  in  a time  sequence.  At  the  beginning  stage  of  switching  to  DDW.  the 
release  of  Ca  was  font  the  bulk  solution,  which  was  confirmed  by  the  fact  that  Ca 
concentrations  and  pH  in  the  effluent  at  pore  volume  =1  were  approximately  equal  to 
those  of  the  influent  ([Ca]=  0.01  M,  pH=6.95)  (Figure  5-4).  Secondly,  Ca  diffusion  from 

number  of  pore  volumes  increased  (between  the  two  arrows  in  Figure  5-4)..  When  Ca  in 
the  bulk  solution  was  depleted.  Ca  cations  moved  against  electrostatic  attraction  away 
from  the  surfaces  driven  by  osmotic  pressure,  extending  EDLs  surrounding  colloids  and 
solid  phases  to  mobilize  colloids.  Obviously,  any  significant  expanding  of  EDL  and 

released  into  bulk  solution.  In  this  case,  the  excess  Ca!'  in  the  diluted  bulk  solution  may 
be  combined  with  hydroxyl  to  release  proton  (Ca1’  + HjO  «-»  CaOH'  + H*)  (Osier  and 
Shainberg,  1979),  resulting  in  a decrease  in  pH.  Between  pore  volumes  6-17,  significant 

colloid  concentrations  (Figure  5-4).  This  is  consistent  with  the  mechanism  we  proposed 
above.  Generally  Ca  hydrolysis  is  w eak  in  solution,  it  is  thus  not  likely  for  it  to  bring  the 


1 5-4.  In 
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this  study,  however,  the  deficiency  of  negative  charge  existed  in  the  buik  solution,  hence 
it  was  possible  that  the  charge  deficiency  may  further  drive  dissociation  of  H2O,  releasing 
protons.  Thirdly,  desorption  ofCa  from  colloid  surfaces  and  stationary  solid  phases 
occurred  from  1 8 to  23  pore  volumes  (afier  the  second  arrow  in  Figure  5-4).  This  was  a 
slow  process  and  therefore  occurred  at  high  pore  volumes,  resulting  in  an  increase  in 
repulsive  forces  between  colloids  and  stationary  solid  phases.  An  increase  in  Ca 
concentration  and  in  pH  were  observed  in  the  effluent  in  pore  volumes  between  17-23 
(Figure  5-4),  In  essence,  this  process  is  equivalent  to  sorbed  Ca  displaced  by  H*. 
However,  a direct  validation  is  needed. 

Obviously,  dissolution  of  Ca  minemls  in  the  soil  may  also  play  a role  in 
regulating  the  release  of  Ca  from  soil  columns,  especially  during  pore  volumes  of  18  to 
23.  During  pore  volumes  of  3 to  23,  however,  it  was  impossible  for  dissolution  to  play  a 
significant  role  in  Ca  releasing,  because  the  pulse-like  feature  of  Ca  releasing  cannot  be 
explained  by  dissolution  of  Ca  minerals.  Of  course,  overlapping  between  the  different 
stages  may  occur.  However,  each  individual  mechanism  may  dominate  a particular  range 
of  leaching  (Figure  5-4). 

As  discussed  above,  the  colloid  elution  curves  can  be  mainly  divided  into  two 
regions  between  pore  volumes  of  3-23,  which  were  regulated  by  Ca  diffusion  or 
desorption  (Figure  5-3).  This  suggested  that  a wide  range  of  surface  heterogeneity 
existed  in  soil  colloids  and  stationary  solid  phases.  In  theory,  different  from  the  colloids 
mobilized  by  Ca  release  via  desorption,  those  mobilized  by  diffusion  must  contain  less 
amount  of  reactive  groups  such  as  >COiH.  -OH  on  the  surfaces.  However,  there  was  no 
direct  evidence  to  support  this  conclusion  in  this  study.  In  the  range  of  Ca  diffusion. 
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colloid  release  showed  a clear  elution  pattern.  In  a well-controlled  system,  Nocito-Gobel 
ct  ol.  (1996)  have  demonstrated  that  polystyrene  colloids  are  released  from  the  column 
packed  with  silica  sand  in  a similar  manner,  i.e.  as  pore  volumes  increased  colloid 
concentrations  appeared  in  two  pulse-like  peaks.  They  attributed  these  to  the  inherent 
heterogeneity  of  the  silica  sand  surface,  and  suggested  the  two  separate  peaks  are  most 
likely  related  to  the  primary  and  secondary  deposition  of  colloids.  In  our  study,  the 
interactions  between  colloids  and  stationary  solid  phases  are  much  more  complicated 
because  of  the  heterogeneity  of  both  colloids  and  stationary  solid  phases.  Therefore,  one 

manner  since  each  group  of  colloids,  which  possess  similar  surface  properties,  would 
appear  in  two  peaks  in  effluents,  which  is  caused  by  colloids  releasing  from  the  primary 
and  secondary  deposition.  This  is  consistent  with  the  trend  at  pore  volumes  less  than  1 7 
shown  in  Figure  5-3. 

As  shown  in  Figure  5-3,  colloid  mobility  changed  with  incubation  time. 

However,  the  peak  position  of  colloid  concentrations  under  different  incubation  basically 
matched  each  other,  even  though  some  shifted  slightly  (Figure  5-3).  This  suggested  that 
incubation  did  not  significantly  change  the  major  types  of  mobile  colloids  and  their 
mobility  rather  it  affected  the  mobile  population  of  each  type  colloid  when  they  were 

Colloid-facilitated  Pb  mobility 

Similar  to  colloidal  Fc  and  Al.  colloidal  Pb  concentrations  varied  with  incubation 


(Figure  5-3).  Interestingly, 


significant  differ 


5-3),  while  the  concentrations  of  colloidal  Pb  decreased  approximately  4 times  after  20-d 
than  3-d  incubation  (Figure  5-3).  This  suggested  that  Pb-bearing  colloids  tended  to  be 
immobilized  or  the  tendency  of  association  of  Pb  with  mobile  colloids  decreased  with 
water-flooded  incubation.  However,  a recent  research  (Krelzschmar  et  al..  1997) 
suggests  that  the  association  of  Pb  with  colloids  does  not  affect  colloid  mobility 
significantly,  especially  at  low  surface  coverage.  This  may  exclude  the  first  possibility. 

It  is  well  known  that  Pb  is  most  likely  sorbed  onto  Fe  oxides  (Benjamin  and 
Leckie,  1981).  However,  when  the  environment  becomes  reducing.  Pb  is  released  as  Fe 
oxides  dissolve  (Gambrell  et  al..  1980).  This  suggests  that  Pb  carried  by  mobile  colloids 
may  redistribute  with  incubation.  One  possibility  is  that  the  released  Pb  resorbed  onto  the 
surface  sites  bicarbonate  (>COjH)  of  carbonate  minerals  (Zachara  et  al.,  1 991 ),  and  the 
particles  bearing  carbonate  mineral  may  carry  positive  charge  and  not  tend  to  be 
mobilized  as  discussed  before.  Therefore,  colloidal  Pb  decreased  with  incubation. 
However,  further  study  is  needed  to  verify  this  mechanism. 

Conclusion 

In  this  Pb-contaminntcd  soil,  colloid  mobility  decreased  with  incubation,  and  it 
was  positively  and  inversely  correlated  to  the  concentrations  of  dissolved  Fe  and  Ca  in 
the  pore  water,  respectively.  This  cannot  be  understood  by  either  the  decementing  effect 
from  reducing  dissolution  of  Fe  oxides  and  charge  heterogeneity  rising  Horn  them. 
Instead,  it  is  possible  that  the  decrease  of  colloid  mobility  with  incubation  was  caused  by 
an  increase  of  positive-charged  patches  of  carbonate  surfaces,  which  is  determined  by 
Ca2*  concentration  in  pore  water. 
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A elution  pattern  of  colloidal  Fc  and  Al  concentration  changes  with  pore  volumes 
was  observed.  In  general,  a peak  of  colloid  release  accompanied  a higher  dissolved  Ca 
concentration.  Three  possible  mechanisms,  i.e.,  bulk  solution  displacement,  diffusion  out 
of  EDL  and  desorption  from  surfaces,  were  proposed  and  played  a major  role  in  a time 
sequence  in  Ca  release.  These  mechanisms  may  explain  the  chromatographic  pattern  of 
colloidal  Al  and  Fe  concentration  changes  with  incubation  time.  Colloid  particles  can 
facilitate  Pb  transport,  however,  the  ability  decreased  with  water-flooded  incubation  for 
this  soil.  It  is  possible  that  Pb  redistribution  from  Fe  oxides  to  carbonates  occurred 
during  incubation. 


CHAPTER  6 

RELEASE  AND  DISPERSABIL1TY  OF  COLLOIDS  IN  TWO  CONTAMINATED 
SOILS 

Inlroducnon 

certain  conditions  they  may  be  mobilized  and  facilitate  transport  of  heavy  metals  in  soils 
(McCarthy  and  Zachara,  1989;  Newman  et  al„  1993),  possibly  imposing  a threat  on  the 

Net  release  of  colloids  from  a soil  is  mainly  determined  by  rates  of  two  opposite 
processes:  colloid  release  (Kallay  et  al„  1987)  and  deposition  (Adamczyk  et  al„  1983). 

colloids  and  stationary  solid  phases,  which  can  be  described  by  the  DLVO  theory 
(Dcrjaguin-Landau-Verwey-Overbeek ) (Israelachvili,  1992).  Based  on  the  DLVO 
theory,  the  net  surface  interaction  energy  is  equal  to  the  sum  of  the  interactions  between 
electrical  double  layer  (EDL)  and  van  der  Waals  forces,  which  vary  with  the  separation 
distance  between  colloids  and  stationary  solid  phases  (Figure  2-1).  As  colloids  move 
close  to  stationary  solid  phases,  the  net  surface  interaction  energy  experiences  a 
secondary  minimum  first,  and  then  a primary  minimum  (9mm)  after  a maximum  energy 
barrier  ($„»*)•  After  that  point  ($„,;„),  further  movement  of  colloids  towards  the  stationary 
solid  phase  causes  drastic  increases  in  the  interaction  energy,  which  is  referred  to  as  the 
Bom  repulsion.  Theoretically,  particles  are  deposited  when  the  net  attractive  forces  are  at 
either  the  primary  or  secondary  energy  minimum.  Colloid  deposition  is  considered  as  a 
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two-step  process:  transport  of  colloids  from  bulk  solution  to  the  stationary  solid  surfaces 
and  then  their  attachment  to  the  surfaces,  which  depends  upon  the  nature  of  particle- 
surface  interaction.  Colloid  attachment  rate  is  exponentially  related  to  the  maximum 
energy  barrier  ■>«»«  (Ruckenslein  and  Prievc,  1976).  Similar  to  colloid  deposition,  colloid 
release  also  consists  of  two  consecutive  steps:  detachment  from  the  surfaces  and  transport 
to  bulk  solution,  and  the  detachment  rate  is  exponentially  related  to  the  depth  of  the 
energy  well  (4mn-4min)  (Figure  2-1). 

To  validate  the  DLVO  theory'  in  describing  colloid  dispersion  and  coagulation  in  a 
suspension,  colloid  stability  test  using  batch  experiments  was  developed  (Reerink  and 
Overbeek,  1954a  and  b).  Such  a test  has  been  partially  adopted  to  evaluate  the 
interactions  among  soil  colloids.  i.e.  water  dispersability  of  soil  colloids  (Miller  et  al., 
1990;  Shainberg  et  al„  1981).  In  general,  water  dispersability  of soil  colloids  decreases 
with  ionic  strength  at  a given  pH.  There  is  a critical  flocculation  concentration  (CFC)  of 
electrolyte  for  a given  soil  suspension,  above  which  flocculation  occurs.  There  has  been 
effort  to  relate  this  CFC-based  dispersability  of  soil  colloids  to  their  mobility  in  soil 
(Kaplan  et  al.,  1996;  Seta  and  Karathanasis,  1997,  Shainberg  et  al.,  1981a  and  b),  i.e. 
applying  batch-experiment  based  results  to  predict  colloid  behaviors  in  column 
experiment  since  colloid  stability  tests  are  normally  conducted  in  batch  experiments 
whereas  colloid  mobility  tests  are  done  in  column  experiments.  Soils  with  higher  CFC 
may  have  greater  colloid  mobility  since  colloids  can  be  more  stable  in  a suspended  state 
at  a given  condition.  However,  no  consistent  result  was  achieved  when  the  CFC  of  soil 
dispersive  clay  was  used  to  predict  soil  colloid  mobility  (Kaplan  et  al.,  1 996).  They 
summarized  that  CFC  predicts  colloid  dispersability,  but  not  mobility  when  size  straining 


is  dominoni  in  controlling  colloid  deposition,  which  occurs  when  colloid  particles  are 
larger  than  pore  opening  in  porous  media.  Different  from  the  colloid  stability  test,  in 
which  colloid  deposition  (aggregation)  is  dominant,  colloid  mobility  in  a soil  is 
determined  by  the  net  release  rate  of  colloids  at  a given  flow  condition  as  shown  in  Eq. 
(6-1). 

Release  Transport 

Immobile  colloids  c ^ > Mobile  colloids  «c » Colloids  in  leachate  (6-1) 

Colloid  release  can  be  induced  by  a rapid  reduction  of  influent  ionic  strength  in 
well-defined  porous  media  (Ruckenstein  and  Prieve.  1 976;  Roy  and  Dzomk,  1 996a; 
Nocito-Gobel  and  Tobiason,  1996)  and  soil  (Seaman  et  al.,  1995;  Jacobsen  et  al„  1997). 
However,  this  may  not  always  be  the  cose  since  the  opposite  result  was  also  observed  by 
Kallay  et  al.  ( 1 987).  They  attributed  to  the  elimination  of  colloid  redeposition  in  short 
columns  as  a possible  cause  for  their  observation.  Hus  suggests  that  the  net  colloid 
released  often  is  a result  of  both  colloid  deposition  and  release.  The  effect  of  ionic 
strength  reduction  on  colloid  release  has  been  attributed  to  the  repulsion  between  colloids 
and  porous  media  arising  from  EDLs  expansion,  resulting  m a reduction  in  the  depth  of 
the  energy  well  (Kallay  ct  al.,  1987).  Different  from  colloid  deposition,  however,  colloid 
release  is  more  difficult  to  evaluate  both  theoretically  and  experimentally.  Using 
hematite  colloids  and  quartz  columns,  Ryan  and  Gschwend(  1 994a)  have  shown  that  the 
transport  of  detached  colloids  to  the  bulk  solution  is  the  rale-limiting  step  during  rapid 
colloid  release  when  the  maximum  energy  barrier  {(w,)  approaches  zero  (Figure  2-1).  In 
soil,  a similar  phenomenon  has  been  observed  by  Jacobsen  et  al.  (1997)  in  an  intact 
column  experiment  in  which  the  colloids  were  leached  with  tap  water.  By  combining  this 
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with  ihe  discussion  that  the  net  colloid  release  rate  is  determined  by  the  relative  rates  of 
colloid  telease  and  deposition  as  shown  in  Eq.  (6-1 ),  one  may  conclude  that,  at  a rapid 
release  condition,  the  net  colloid  release  rate  is  determined  by  the  rates  of  the  net  colloid 
diffusion  to  bulk  solution  and  attachment  onto  stationary  solid  surfaces. 

Based  on  the  above  discussion,  we  hypothesize  that,  under  a given  water  flow 
condition  in  a soil,  the  net  colloid  release  rate  may  be  controlled  primarily  by  colloid 
attachment  rate,  which  can  be  evaluated  using  the  colloid  stability  ratio  (W)  in 
suspensions.  Colloid  stability  ratio  is  the  ratio  of  colloid  flocculation  rale  of  a suspension 
at  normal  observation  condition  (slow  coagulation)  to  that  under  rapid  flocculation 
condition  when  the  energy  barrier  (jta,,)  = 0 (Reerink  and  Overbeek,  1 954).  Under  slow 
flocculation  condition,  only  a fraction  (W)  of  colloid  collisions  leads  to  coagulation.  In 
fact,  CFC  represents  only  a specific  case  of  flocculation  kinetic  where  W =1.  Therefore, 
W instead  of  CFC  will  be  used  to  evaluate  dispcrabililty  of  soil  colloids  in  this  paper. 
Assume  that  the  rate  of  flocculation  is  related  to  the  changes  in  light  absorbency  of 
solution,  W can  be  obtained  by  dividing  the  rate  of  absorbency  change  observed  at  rapid 
coagulation  condition  by  that  at  slow  coagulation  condition  (Reerink  and  Overbeek, 
1954).  In  soil,  measurement  of  light  absorbency  has  been  widely  used  to  evaluate  CFC 
of  colloids  (Heil  and  Sposito,  1993;  Kretzscltmar  et  al„  1997). 

We  have  demonstrated  that  colloid  release  from  soil  columns  is  influenced  by  soil 
redox  status  (Chapter  5),  which  cannot  be  simply  related  to  the  decementing  effect  of  Fe 
oxides.  Instead,  colloid  release  can  be  related  to  changes  in  major  potential  determining 
ions  with  changes  in  soil  redox  potential.  In  this  paper,  a different  approach  will  be  taken. 
We  will  1)  evaluate  colloid  dispersnbility  by  comparing  relative  colloid  stability  ratios; 
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and  2)  determine  the  relation  between  the  relative  colloid  stability  ratio  and  colloid 
release  in  two  selected  contaminated  soils. 

Materials  and  Methods 

Column  preparation 

The  two  Pb-contaminated  soils  used  in  this  study  were  from  Montreal.  Canada 
and  Tampa,  Florida,  which  were  exposed  to  Pb-battcry  recycling  and  concurrent  smelting 

Soil  columns  were  prepared  as  previously  described  (Chapter  4)  and  were  briefly 
outlined  here.  5.0  g of  acid-washed  sand  (20-30  mesh)  was  placed  in  the  bottom  of  60- 
mL  columns  (13  cm  in  length  x 2.6  cm  in  diameter),  then  30  g air-dried  soil  was  placed 
on  top  of  the  sand  layer.  This  resulted  in  2.6x5.2  and  x7.8-cm  soil  columns  for  the 
Montreal  and  Tampa  soils,  respectively.  The  soil  columns  were  set  vertically  and 
prewetted  by  pumping  deionized  distilled  water  (DDW)  into  the  bottom  of  the  columns 

incubated  under  water-flooded  condition  for  approximately  3, 20  and  80  d for  the  Tampa 
soil  and  3, 20, 30, 40, 50,  and  60  d for  the  Montreal  soil.  Then  the  standing  water  on  top 

with  rubber  stoppers.  The  rubber  stoppers  on  bottom  of  the  columns  were  then  connected 
to  a needle  as  inlet  for  influent.  Then,  0.01  M CaClj  solution  was  pumped  through  the 
columns  until  Ca  concentration  in  the  effluent  was  approximately  equal  to  that  in  the 
influent,  i.e.  when  CaCl;  completely  displaced  the  soil  solution  in  the  column.  The  soil  at 


prepared  (three  replicates  for  each  treatment),  18  of  which  were  used  for  the  column 
leaching  test  (two  replicates)  and  nine  of  which  for  the  water  dispersability  test. 

Column  leaching  lest 

After  the  soil  columns  were  saturated  with  Ca,  the  influent  in  1 8 of  the  soil 
columns  was  switched  from  0.01  M CaCl2  solution  to  DDW  at  a pumping  rate  of  1 .20 
mL  per  minute.  The  effluent  from  the  top  of  each  of  the  column  was  collected  with  a 
(radio nation  collector.  Concentrations  of  Fe,  Al,  and  organic  C were  determined  in  the 
unfiltered  and  filtered  (0.2  pm)  solutions.  The  differences  between  the  two  fractions 
were  considered  as  colloidal  Fe,  Al,  and  organic  C.  Based  on  preliminary  experiment, 
colloidal  Fe  and  Al  ore  the  major  components  of  colloids  for  the  Tampa  soil  and  organic 

and  colloidal  organic  C were  used  to  approximate  colloid  concentrations  in  these  two 
soils.  Cumulative  amount  of  colloids  leached  after  the  first  10  pore  volumes  was 
calculated. 


After  the  soil  columns  were  saturated  with  Ca,  9 soil  columns  (each  for  each 
treatment)  was  unpacked  for  the  water  dispersability  test  The  soil  core  inside  the  column 

into  16  polycarbonate  centrifuge  tubes  (40-mL)  with  half  containing  40-mLof0.0l  M 
CaClj  and  half  containing  40-mL  of  0.06  M NaCI  solution.  All  tubes  were  placed  in  an 
ultrasonic  system  (FS28H,  Fisher  Scientific)  for  3 min  at  room  temperature.  Preliminary 
research  showed  that  soil  dispersion  increased  with  time  initially  then  leveled  off  after 


reaching  a maximum.  After  3 minutes  of  sonification,  90  % of  dispersible  panicles  in 
both  soils  were  dispersed  (data  not  shown).  The  resulted  suspensions  were  centrifuged 
(Beckman.  Model  J2-2 1 ) at  1 000  rpm  for  2 minutes  at  20  “C.  This  process  supposedly 
removes  all  panicles  > 2 pm  with  a panicle  density  of  Z4  g cm'!  (Tanner  and  Jackson. 
1947).  The  top  portion  of  the  supernatant  (~20  mL)  in  each  tube  was  removed  carefully 
from  the  tube  using  a pipette.  Those  containing  the  same  electrolyte  solution  (CaClj  or 

disposability  test 

developed  by  modifying  those  of  Thcllier  and  Sposito  (1989),  Kretzschmar  et  al.  (1997) 
and  Seta  and  Karathanasis  (1997).  A total  of  14  polystyrene  test  tubes  (13-mL)  were 

capped.  A subsample  of  2.0  mL  was  removed  at  2 cm  below  the  surface  from  two  test 
tubes  (duplicate)  after  2, 4, 6. 8 . 10. 20  and  24  h.  The  light  absorbency  of  the  sub- 
samples  was  then  read  immediately  at  525  nm  on  a Shimadzu  UV-160 
Spectrophotometer.  All  procedures  were  conducted  under  N;  atmosphere. 

Estimation  of  relative  colloid  stability  ratio  (RW) 

Reerink  and  Overbeek  (1954)  determined  the  rate  of  coagulation  of  silver  iodide 
particles  from  light  absorbency.  In  a short  period  of  time  (—15  min),  light  scattering 
increases  as  particle  sizes  increases  due  to  aggregation,  and  so  does  the  reading  of 


absorbency.  Thispn 


I»-eV0JCo 

I,-e£V,(t)!Ci(t) 
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(6-2) 

(6-3) 

Where  lo,  Vo  and  Co  are  inilial  absorbency,  average  volume  per  particle  and  particle 
concentration,  respectively;  It  Vi(t)andCi  (t)  are  absorbency,  volume  and  concentration 
of  an  aggregated  particle  consisting  of  i particles  at  time  t,  respectively;  and  c is  an 
optical  constant  For  soil  colloids,  however,  their  aggregation  kinetics  is  much  more 
complicated  because  of  the  heterogeneity  in  colloid  composition,  mineralogy  and  size. 

To  our  knowledge,  using  colloid  flocculation  rate  to  estimate  stability  ratio  has  not  been 

to  minimize  the  complicated  details  of  soil  colloid  aggregation  with  time.  In  our 
experiment  longer  settling  time  (>  2 h)  resulted  in  large,  aggregated  particles  settling  out 
of  the  sampling  zone  in  the  test  tubes  via  gravity,  and  thus  the  absorbency  decreases  with 

of  light  absorbency)  with  time  in  our  method  is  different  from  the  one  established  by 
Reerink  and  Overbeek  (1 954).  However,  assuming  the  size  and  optical  constant  of  the 
suspending  particles  in  the  sampling  zone  of  the  test  tube  are  constant  (settling  time  > 2 
h),  the  basic  rationale  for  the  method  of  Reerink  and  Overbeek  ( 1 954)  is  valid  in  our 
system.  In  our  system.  V^t)  could  be  treated  as  a constant  V,  then  Eq  (6-3)  can  be  written 


I,=sV2C(t) 


(6-4) 


Eq  (6-4)  states  that  flocculation  rate  can  be 
absorbency  with  time  in  a linear  stage  of  lie 


! obtained  by  i 


However,  we  found  that  colloid  flocculation  rale  varied  with  time  in  a 
complicated  manner,  which  is  similar  to  those  reported  by  Fcrretti  et  al.  (1997).  To 
minimize  uncertainty  in  calculating  colloid  flocculation  rate,  we  choose  a linear  region  of 
absorbency.  Two  types  of  absorbency-time  curves  were  observed  in  our  study:  1) 
Absorbency  decreased  somewhat  exponentially  with  time  at  initial  measurement.  The 
rate  was  calculated  using  the  slope  of  the  linear  region  (6-10  h)  after  an  initial  sharp 
decrease  in  absorbency  (Figure  6-la);  and  2)  No  distinct  "initial  drop”  period  was 
detectable.  The  rate  was  determined  by  linear  regression  of  the  data  obtained  in  the  fust 
8 h of  measurement  (Figure  6-lb).  The  above  procedures  yield  an  apparent  flocculation 
rate  in  the  linear  stage.  Although  it  is  subject  to  bias,  it  is  useful  to  calculate  and  compare 
such  apparent  flocculation  rates  based  on  absorbency-time  curves  of  a soil  under  different 
conditions. 

Based  on  the  definition  of  colloid  stability  ratio  (W)  of  Rccrink  and  Overbeek 
(1954),  W can  be  expressed  as  follows: 

W « Co!  U(  dl,/dt)  (6-5) 

Where  dVdt  is  the  slope  of  the  absorbency-time  curve  in  a linear  stage  and  L is  the  length 
of  the  light  path  through  the  suspension.  In  addition.  W can  be  estimated  using  the  ratio 
of  the  rates  of  rapid  flocculation  to  slow  flocculation.  Rapid  colloid  flocculation  rate,  in 
this  paper,  was  measured  using  a Ca-saturated  soil  in  0.01  M CaClj  electrolyte 
background  solution.  Colloids  in  the  Ca-saturated  soils  were  mobilized  when  the  influent 
was  switched  from  CaClj  solution  to  DDW.  However,  in  the  absence  of  background 
electrolytes  the  measurement  of  the  light  absorbency  in  a soil  suspension  is  subject  to  a 
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Figure  6-1  Schematic  representation  of  typical  absorbency-time  curves  observed.  Linear 
regression  was  used  to  calculate  an  apparent  flocculation  rate  indicated  by  solid  line.  The 
curves  are  not  drawn  in  the  same  time  and  absorbancy  scales.  Curve  "a"  stands  for  those 
observed  for  the  Montreal  soil  in  0.06  NaCI  M solution,  and  Curve  “b”  for  those  for  the 
Tampa  soil  in  0.06  M NaCI  and  in  0.01  M CaCfe  solutions,  and  the  Montreal  soil  in  0.01 
MCaCb  solution. 
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great  uncertainty.  Therefore.  0.06  MNaCI  solution  was  used  as  background  electrolytes. 
Slow  colloid  flocculation  rate  was  thus  measured  using  Ca-saturated  soil  in  0.06  M NaCl 
background  electrolyte  solution.  Same  ionic  strength  in  both  background  solutions  was 
designed  to  minimize  the  effect  of  soil  mineral  dissolution  caused  by  electrostatic 
interaction  between  solution  and  mineral  on  colloid  dispersion.  This  effect  has  to  be 
considered  here  since  the  residence  time  of  solution  is  much  longer  than  that  in  column 
leaching.  Therefore,  for  a soil  incubated  after  a given  lime,  its  relative  W (RW)  is 
defined  as 

RW=  Wm/Wca'c  «dVdt)c.(Co)1N.l/{(dVdt)N.  <C„)!c.)  (6-6) 

Where  the  subscripts  Na  and  Ca  arc  the  measurement  taken  in  0.06  NaCl  and  0.01  CaCh 
M solutions,  respectively.  Combining  Eq  (4-2)  and  (4-6)  yields 

RW  «c  ((dI/dt)c,(Io)!N«)'l(dl|/dt)N,  (Ioj’ci  (6-7) 

Eq  (6-7)  was  used  in  this  paper  to  estimate  RW  of  a soil  incubated  after  a given  time  in 

0. 06  M NaCl  solution  relative  to  in  0.01  M CaCI;  solution.  Based  on  the  definition 
above,  RW's  relative  magnitude  represents  the  colloid's  tendency  to  disperse  in  a soil, 

1. e.  a larger  number  indicates  greater  tendency  to  be  dispersed. 


Result  and  Discussion 

Figure  6-2  is  a typical  example  of  absorbency-time  curves  observed  for  the  Ca-saturated 
Tampa  soil  in  0.06  M NaCl  background  electrolyte  solution.  The  overall  trend  of  these 


curves  is  similar  to  Figure  6-1  b.  Therefore,  the  slope,  (dVdt)u,  and  initial  absorbency  (lo) 
Na  were  determined  from  the  initial  linear  stage  (2-8  h).  Table  6-2  lists  alt  the  parameters 
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Figure  6-2  Absorbcncy-lin 
0.06  MNaCI  solution. 


■ves  observed  in  Ca-salurnled  Tampa  soil  suspended  in 
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s.  Clearly,  RW  varied  from  0.S  to  17.7 


needed  lo  calculate  RW  via  Eq  (6-7)  and  RW  values 
and  0.86  to  6.9  for  the  Tampa  and  Montreal  soils  with  incubation  time,  and  reached  a 
maximum  at  20  d and  3 d,  respectively  (Table  6-2).  This  is  consistent  with 
our  earlier  findings  (Chapter  5)  that  incubation  significantly  influences  colloid  release 

The  basic  mechanism  for  dispersion  of  the  Ca-saturated  soil  in  the  Na  solution 
was  that  Ca  in  the  Stem  layer  and  in  the  diffiisc  layer  was  displaced  by  Na  in  bulk 
solution,  leading  to  electrical  repulsion  as  a result  of  the  EDL  expansion  (Roy  and 
Dzombk,  1996b).  As  moreNa  replaced  Caon  the  exchangeable  sites,  the  soil  became 
more  dispersed.  i.e.  more  colloids  were  released.  As  a result,  RW,  a measure  of  relative 
colloid  stability,  will  increase  with  an  increase  in  Na  concentrations  on  the  exchangeable 
sites.  This  ion  exchange  process  can  be  described  by  the  Gaincs-Thomas  equation: 

Kn*c,  ={[Na']IEaxi)/{[CaJ*]E3N1x)  (6-8) 

Where  Kitec.is  Gaines-Thomas  selectivity  coefficient,  [Na*]  and  [Ca!‘]  represent  Na  and 
Ca  molar  concentrations,  and  Eox2  and  Exax  are  the  equivalent  fractions  of  Ca2*  and  Na* 
in  the  exchange  phase  of  soil  particles  which  is  denoted  here  as  X.  Larger  ■(■<»<*  value 
indicates  greater  selectivity  of  soil  exchangeable  sites  for  Ca.  Since  only  1 g of  Ca- 
saturated  soil  was  mixed  with  40  mL  of  0.06  M NaCl  solution,  it  is  expected  that  [Na*] 
»[Ca2*],  Eq  (4-8)  can  be  rewrinen  as  (McBride,  1994): 

Ktwc.  =(Eax!/  E2*.x}/|nWmc.)MT  . (6-9) 

Where  ms,  and  me.  are  moles  of  Ca2*  and  Na'  in  the  solution  phase  and  Mt  is  the  total 
molarity  of  the  solution.  In  our  experiment,  Mr  can  be  considered  a constant.  Thus, 
significant  increase  in  RW  (Table  6-2)  at  20  d and  3 d incubation  for  the  two  soils  was 
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most  likely  caused  by  the  decrease  in  Ec.xVE!n>x-  Therefore,  Eq.  (6-9)  suggests  that  Kn» 
on  cannot  be  a constant  under  difTerent  incubation  times  especially  when  there  is 
significant  change  in  the  ratio  of  Eoa/  E2n«x.  This  is  inferred  by  the  fact  that  the 
decrease  in  the  ratio  of  Eooq/  E2nix  cannot  be  compensated  by  the  increase  in  the  ratio  of 
nWmca  on  the  right  side  of  Eq.  (6-9).  It  is  possible  that  significant  decrease  in  the 
selectivity  coefficient  Ku»ca  led  to  increase  in  RW  for  the  two  soils  under  different 
incubation  times.  Thus,  it  seemed  that  RW  and  Kn*c.  were  inversely  related  in  this 
experiment.  In  fact,  Kn*cb  is  influenced  by  many  factors  such  as  pH,  ionic  strength,  and 
exchangeable  site  geometry  (McBride,  1994),  However,  to  our  knowledge,  there  is  no 
miormauon  available  about  how  incubation  time  affects  Kn*ch  in  a soil.  In  addition,  the 
heterogeneity  of  exchangeable  sites  in  soil  colloids  make  it  difficult  to  describe  the 
characteristics  of  soil  particles  when  Knur,  was  low.  In  the  stability  lest  we  actually 
examined  only  a small  portion  of  the  soil  colloids  in  a lest  tube,  the  heterogeneity  of  the 
soil  colloids  thus  must  be  taken  into  account.  However,  with  the  current  knowledge  it  is 
difficult  to  further  discuss  this  topic  with  confidence.  Nevertheless,  high  value  of  RW  for 
both  soils  may  be  due  to  the  low  value  of  Knm:„  which  resulted  from  the  water-flooding 
incubation.  In  Chapter  5.  we  suggested  that  the  surface  properties  of  colloid  particles  may 
be  altered  with  incubation,  and  thus  changing  colloid  mobility.  This  is  consistent  with  the 
result  we  present  here  conceptually. 

Cumulative  mobile  colloids  in  the  first  1 0 pore  volumes  in  the  two  contaminated 
soils  under  different  incubation  times  were  plotted  against  its  corresponding  relative 
colloid  stability  ratios  (Figure  6-3).  Since  organic  matter,  and  AI  and  Fe  are  the  major 
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Figure  4-3.  Relation  of  cumulative  colloids  and  relative  stability  ratio  (RW) 
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components  of  the  mobile  colloids  in  the  Tampa  and  Montreal  soil,  respectively,  they 
were  used  to  approximate  the  mobile  colloids  released  from  these  two  soils.  Generally, 
RW  was  positively  related  to  the  amount  of  cumulative  colloids  (~  10  pore  volumes)  in 
leachates  for  both  soils  (Figure  6-3).  Interestingly,  though  the  characteristics  of  the  two 
soils  were  significantly  different  (Table  6-1 ) and  colloid  composition  in  leachates  were 
different  (inorganic  Al  (Fe  vs.  organic  matter),  the  relation  between  cumulative  colloids 
and  RW  for  the  two  soils  was  similar.  It  is  possible  that  the  decreases  in  the  amount  of 
cumulative  colloid  with  WR  for  the  two  soils  in  Figure  6-3  was  caused  by  transition  from 
dispensability  control  to  size  straining  control  of  existing  mobile  colloid  as  RW  increased. 

Our  results  suggest  that  RW  may  be  used  to  predict  overall  colloid  release  in 
soils.  However,  theoretically,  colloid  release  is  determined  by  the  magnitude  of  the 
energy  barrier  (^mn),  whereas  colloid  deposition  is  determined  by  the  depth  of  the  energy 
well  (Figure  2-1).  Using  hematite  colloids  and  quartz  grains.  Ryan  and 

Gschwend  ( 1 994a)  found  that  both  colloid  release  rates  and  the  magnitude  of  the  energy 
barriers  decreased  os  ionic  strength  increased  under  rapid  colloid  release  condition.  They 
suggested  that  rapid  colloid  release  is  controlled  by  transport  of  detached  colloids  to  the 
bulk  solution  instead  of  detachment  itself  since  the  energy  barrier  has  vanished  from  the 
potential  energy  profile  due  to  changes  in  solution  chemistry.  A similar  phenomenon  has 
been  observed  by  Jacobsen  ct  al.  (1997)  using  intact  soil  column  experiments  in  which 
the  natural  particles  were  leached  with  tap  water.  Their  results  showed  that  panicle 
release  from  the  columns  is  not  affected  by  an  increase  in  flow  rate.  Further,  their  plot  of 
accumulated  amount  of  released  panicles  versus  square  root  of  time  show  a fairly  linear 
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relation,  implying  diffusion-limited  kinetics  is  a dominant  process.  One  condition  for  this 
to  occur  is  that  the  maximum  energy  barrier  of  mobile  colloids  must  vanish.  This  may  be 
very  common  in  soil,  which  is  consistent  with  the  transient  phenomenon  we  discussed  in 
Chapter  2.  Under  such  a condition,  the  existing  colloid  concentration  at  a given  residence 
time  and  water  flow  condition  in  a soil  is  controlled  by  colloid  attachment  rate,  i.e„  the 
size  of  the  energy  barrier  ((■«»)•  Therefore  from  this  point  of  view,  it  can  be  understood 
that  the  cumulative  amount  of  released  colloids  in  leachates  can  be  well  related  to  RW, 
even  though  RW  is  only  a measure  of  the  magnitude  of  energy  barrier  not  the  depth  of 
the  energy  well. 

In  fact,  there  has  been  considerable  effort  to  use  colloid  dispersability  to  predict 
colloid  release  in  soil.  Kaplan  et  al.  (1996)  has  shown  that  CFC  of  soil  colloids  predicts 
colloid  mobility  reasonably  well  unless  size  straining  is  dominant  in  controlling  colloid 
deposition,  which  arc  consistent  with  those  reported  by  Seta  and  Kantthanasis  (1997)  and 
in  the  result  presented  here.  However,  different  from  CFC  the  RW  we  defined  in  this 
paper  included  kinetic  information  of  soil  colloid  flocculation  when  the  system  is  subject 
to  a perturbation  in  solution  chemistry.  In  general,  higher  CFC  indicates  greater  colloid 
stability  and  thus  greater  absorbency  of  suspension  at  a given  condition.  However,  CFC 
tailed  to  explain  the  colloid  release  in  our  leaching  test.  For  example,  the  Tampa  soil 
incubated  after  60  d has  the  highest  readings  of  absorbency.  i.e.  higher  tendency  for 
dispersion  (Figure  4-2),  however,  its  cumulative  colloids  was  almost  the  lowest  (Table  6- 
2). 
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Conclusion 

The  release  of  soil  colloids  into  ground  walcr  has  attracted  considerable  attention 
because  of  the  possibility  of  their  association  with  insoluble  contaminants.  Much  effort 
was  spent  to  use  the  parameters  derived  from  colloid  deposition  to  predict  the  overall 
release  of  colloids  in  a soil.  However,  based  on  modern  colloid  theory,  colloid  deposition 
and  release  are  different  processes  and  thus  have  different  dependence  on  the  interaction 
between  surfaces.  In  this  chapter,  we  defined  and  estimated  relative  colloid  stability  ratio 
of  a Ca-saturated  soil  from  absorbency-time  curves  in  0.01  CaCI;  and  0.06  NaCI  M 
solutions.  The  rapid  and  slow  flocculation  were  defined  as  those  in  0.01  CaClr  and  0.06 
NaCI  M solutions,  respectively.  These  are  consistent  with  the  condition  we  used  in  the 
leaching  test  to  some  extent.  The  dispersion  of  a Ca-saturated  soil  in  0.06  NaCI  M 
solution  was  mainly  caused  by  cation  exchange  between  Ca  in  exchangeable  phase  and 
Na  in  bulk  solution.  We  inferred  that  the  significant  increase  in  stability  ratio  for  a soil 
with  incubation  was  probably  caused  by  a decrease  in  Kn*ci.  However,  the  mechanistic 
understanding  is  missing.  Interestingly,  a very  obvious  relation  was  observed  that  overall 
release  of  colloids  from  the  soils  increased  as  the  stability  ratio  increased,  and  the  overall 
trend  is  similar  for  Tampa  and  Montreal  soils.  This  result  suggests  that  the  overall 
release  of  colloids  from  these  soils  could  be  evaluated  using  the  relative  stability  ratio 
qualitatively,  even  though  they  are  different  conceptually. 


CHAPTER  7 
CONCLUSION 

Metal  solubility  and  mobility  in  sail 

Our  data  showed  that  enhanced  Pb  solubility  was  only  observed  at  low  Fe 
partition  index  for  a Pb-spiked  sandy  soil.  Such  a trend  was  also  observed  using  two 
published  data  sets  in  the  literature.  Since  data  used  in  this  paper  were  from  various  soil 
environments  (roadside  sediment  in  Fiance.  Pb  contaminated  soil  in  Poland,  a sandy  soil 
from  Florida),  it  may  be  reasonable  to  conclude  that  Pb  solubility  may  be  related  to  Fe 
partition  index  in  soils. 

For  the  two  Pb-contaminated  soils,  metal  solubility  was  examined  by  analyzing 
the  pore  waters  of  the  incubated  soil  columns  whereas  metal  mobility  was  examined  by 
leaching  the  columns  with  0.01  M CaCIs  solution.  The  data  showed  that  metal  solubility 
in  pore  water  and  mobility  with  CaCl2  solution  was  not  always  directly  related.  There  are 
sufficient  evidence  that  increases  in  CEC  occur  with  Fe  reduction  dissolution.  However 
metal  redistribution  between  solution  and  exchange  sites  with  incubation  is  determined 
not  only  by  the  magnitude  of  CEC,  but  also  the  characteristics  of  metal  ions,  competing 
co-ions,  and  counter-ions  (ligands). 

State  of  Colloid  Deposition  and  Release  in  Soil  and  Their  Association  with  Heavy  Metals 

Soil  has  distinctively  different  features  in  colloid  mobility'  compared  with  well- 
defined  porous  media  or  subsurface  systems.  Soil  colloids  are  also  exposed  to  more 
dynamic  solution  chemistry,  it  is  therefore  more  difficult  to  predict  the  charges  of  soil 
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colloids  since  it  is  determined  by  the  interactions  between  solid  surfaces  and  solution.  In 
general,  soil  particles  are  assemblages  of  crystalline  and  amorphous  minerals,  and 
organic  residuals.  For  such  particles,  the  interface  of  soil  particle-solution  may  be  more 
practical  to  be  described  by  diffuse  double  layer  with  a specific  part,  i.c.,  a charged 
surface,  and  a generic  part,  i.e.,  a diffuse  layer.  The  surface  charge  development  process 
may  be  described  by  partitioning  of  charged  species  between  solution  and  solid  surfaces 
assuming  the  charged  species  are  similar.  Similarly,  association  of  heavy  metals  with 
colloids  may  have  significant  influences  on  colloid  mobility,  which  is  subject  to  dynamic 
changes  in  solution  chemistry.  However,  very  limited  information  is  available  at  the 
present  time. 

Different  from  clean  porous  media,  the  surfaces  of  soil  stationary  phase  consist  of 
deposited  colloids.  Consequently,  the  blocking  effect  of  colloid  deposition  in  soil  is 
significant,  which  may  greatly  enhance  colloid  mobility  in  soil  compared  to  that  in  a 
comparable  clean  porous  medium.  Under  such  a condition,  size  straining  is  dominant 
mechanism  in  soil  colloid  deposition.  Phenomenologically,  this  process  is  time- 
dependent  or  condition-dependent,  which  is  referred  to  as  transient  phenomena  in  colloid 
deposition.  At  the  present  time,  water  dispersible  cloy  extracted  from  a soil  is  widely 
used  to  indicate  colloid  mobility.  However,  it  has  to  be  realized  that  there  are  two 
simultaneous,  opposite  processes:  colloid  deposition  and  release.  In  a given  system,  the 
relative  effect  of  each  on  either  soil  clay  dispersion  in  batch  test  or  colloid  mobility  in  a 
column  experiment  depends-on  colloid  residence  time. 

Despite  shortcomings  of  the  DLVO  (Derjaguin-Landau-Verwey-Ovcrbcek) 
theory,  it  provides  a basic  theoretical  framework  in  describing  soil  colloid  dispersion. 


deposition  and  release,  which  is  simple  and  applicable  conceptually.  It  has  to  be  admitted 
that  the  whole  modem  colloid  theory  (not  only  the  DLVO  theory)  has  found  serious 
difficulty  in  applying  it  to  soil  because  of  heterogeneity  in  both  structure  and  composition 
of  soil  particle  surfaces.  To  determine  the  potential  mobility  of  soil  colloid  in  a soil,  a 
phenomenological  approach  may  be  mote  practical  at  the  present  time.  Soil  colloids 
possess  high  concentrations  of  exchangeable  ions  in  their  EDLs  (high  CEC).  It  has  been 
long  recognized  that  colloid  release  in  soil  is  coupled  by  significant  ion  transfer  between 
EDLs  and  bulk  solution.  Potential  colloid  mobility  may  be  estimated  qualitatively  by 
evaluating  ion  transfer. 


Our  results  showed  that  colloid  mobility  varied  greatly  with  incubation.  For  the 
Montreal  soil,  it  decreased  with  incubation,  and  was  directly  correlated  to  dissolved  Fe 
concentrations  and  inversely  to  Ca  concentrations  in  pore  water.  This  cannot  be 
understood  by  either  the  decementing  effect  from  reducing  dissolution  of  Fe  oxides  and 
charge  heterogeneity  rising  from  Fe  oxides.  Instead,  it  is  possible  that  the  decrease  in 
colloid  mobility  with  incubation  was  probably  caused  by  an  increase  in  positively- 
charged  patches  of  carbonate  surfaces,  which  is  determined  by  Ca‘"  concentration. 

In  the  Montreal  soil,  a pattern  of  colloidal  Fe  and  A1  elution  was  observed.  In 
general,  a peak  colloid  release  was  accompanied  by  a greater  dissolved  Ca  concentration. 
Three  possible  mechanisms.  i.e.,  bulk  solution  displacement,  diffusion  out  of  EDL  and 
desorption  from  surfaces,  were  proposed  to  have  played  a major  role  in  a time  sequence 
in  Ca  release.  These  mechanisms  may  explain  the  pattern  in  changes  of  colloidal  A1  and 
Fe  concentration.  Colloids  panicles  can  facilitate  Pb  transport;  however,  this  more 
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reduced  with  water-flooded  incubation  for  this  soil.  It  is  possible  that  Pb  redistribution 
from  Fe  oxides  to  carbonates  occurred  during  incubation. 

In  addition,  we  deflncd  and  estimated  the  relative  colloid  stability  ratio  of  a Ca- 
saturated.  incubated  soil  from  the  absorbency-time  curves  in  0.01  M CaClj  and  0.06  M 
NaCI  solutions.  -Rapid  and  slow  flocculations  were  defined  as  those  in  0.01  M CaCI;  and 
0.06  M NaCI  solutions,  respectively.  These  arc  consistent  with  the  condition  used  in  the 
leaching  test.  The  dispersion  of  a Ca-saturated  soil  in  0.06  M NaCI  was  mainly  caused 
by  cation  exchange  between  Ca  in  exchange  sites  and  Na  in  bulk  solution.  -Significant 
increase  in  the  stability  ratio  with  incubation  in  a soil  was  probably  caused  by  a decrease 
in  the  selectivity  coefficient  Ks,.cj.  Interestingly,  the  overall  release  of  colloid  from  the 
soils  increased  as  the  stability  ratio  increased,  and  the  overall  trend  is  similar  for  both 
Tampa  and  Montreal  soils.  These  data  suggest  that  the  overall  release  of  colloids  from 
these  soils  may  be  evaluated  qualitatively  using  the  relative  stability  ratio,  even  though 
they  are  different  conceptually. 
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